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PREPARATION OF THIS DOCUMENT 


The author, Emygdio Cadima. now retired, was an FAO scientist in the Fisheries Department 
until 1974, when he returned to the Institute de Investigate das Pescas e do Mar (IPIMAR) in 
Portugal, having also been a Professor at the University of Algarve until 1997. At the end of 
1997 he was the lecturer of a course in Fish Stock Assessment in IPIMAR, which became the 
basis for the preparation of this manual, requested and supported by the Project FAO/DANIDA 
GCP/INT/575/DEN. This manual also incorporates notes from courses in Fish Stock Assessment 
held at several different venues in the world, mainly in Europe, Latin America and Africa. These 
courses had an active collaboration of fisheries scientists from all over the world, especially 
Portugal. These scientists are also co-responsible for the orientation, for the matters treated and 
particularly for the elaboration of the exercises. 

This manual aims to present the basic knowledge on the problems and methods of fish stock 
assessment to young scientists, post-graduate students, and PhD students. This is a scientific area 
in permanent development, where the knowledge of fisheries biology is applied in order to make 
a rational and sustained exploitation of the fishing resources. 

The “Manual of Fish Stock Assessment” is mainly concerned with the theoretical aspects of the 
most used models for fish stock assessment. The practical application (i.e. the exercises solved in 
a spreadsheet), is considered as a complementary part to help the understanding of the theoretical 
matters. 

The editing of the manuscript was made by Siebrcn Venema, manager of 
Project GCP/INT/575/DEN and Ana Maria Caramelo, Fishery Resources Officer in the FAO 
Fisheries Department. 
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Cadima, E.L. 

Fish stock assessment manual. 

FAO Fisheries Technical Paper. No. 393. Rome, F AO. 2003. 1 61 p. 

ABSTRACT 

The manual follows the same order of the lectures in the last course held in 1PIMAR 
(November/December 1997). It starts with an introduction to the mathematical models applied in 
Fish Stock Assessment and some considerations on the importance of fisheries. The need for a 
rational management of the fishing resources is then stressed, this being indispensable for an 
adequate exploitation, aiming at conservation, to occur. The basic assumptions about a model 
and the concepts of different variation rates of a characteristic in relation to time (or to other 
characteristics) are presented, highlighting the most important aspects of the simple and 
exponential linear models which are used in the chapters that follow. After some considerations 
on the concept of cohort, models for the evolution in time of the number and weight of the 
individuals that constitute the cohort are developed, including models for the individual growth 
of the cohort. In the chapter concerning the study of the stock, the fishing pattern and its 
components are defined, the most used models for the stock-recruitment relation are presented, 
as well as the short and long term projections of a stock. With regard to fishing resources 
management, the discussion is focused on the biological reference points (target points, limit 
points and precautionary points) and fisheries regulation measures. The last chapter, which 
presents and discusses theoretical models of fish stock assessment, deals with production models 
(also designated as general production models) and with the long and short-term projections of 
the catches and biomasses. Finally, the general methods of estimating parameters are described 
and some of the most important methods are presented, with special relevance to the cohort 
analysis by age and length. Then a solution of the exercises from the last course held in IPIMAR, 
is presented by the author and the scientist Manuela Azevedo. 
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TO MY FIRST MASTERS AND OLD-TIME FRIENDS 


Ray Beverton 
John Gulland 
Gunnar Saetcrsdal 


PREFACE 


This work is essentially orientated to present an introduction to the mathematical models applied 
to fisheries stock assessment. 

There are several types of courses about the methods used in fish stock assessment. 

One type considers practical application as the main aspect of the course, including the use of 
computer programs. The theoretical aspects are referred to and treated as complementary 
aspects. 

A second type is mainly concerned with the theoretical aspects of the most used models. The 
practical application, considered as the complementary part, facilitates the understanding of the 
theoretical subjects. 

In this work, the second type was adopted and exercises were prepared to be solved in a 
worksheet (Microsoft Excel). The table of contents indicates the exercises corresponding to each 
subject. 

This manual is the result of a series of courses on Fish Stock Assessment held in the following 
places. Portugal : Instituto de Investigate das Pescas e do Mar - IPIMAR (ex-INIP) in Lisbon. 
Faculdade de Ciencias de Lisboa, University of Algarve and Instituto de Ciencias Biomedicas de 
Abel Salazar in Oporto. Other courses were held at Instituto de Investigayao das Pescas in Cape 
Verde, at the Centro de Investigayao Pesqucira in Angola, at the Instituto de Investigayao das 
Pescas in Mozambique, at the Centro de Investigacion Pesquera - CIP in Cuba, at the Instituto 
del Mar del Peru - IMARPE in Peru, at the Instituto Espanol de Oceanografia - IEO (Vigo and 
Malaga - Spain). It is also a result of some lectures integrated into cooperation courses held in 
several countries and organized by FAO, by SIDA (Sweden), by NORAD (Norway) and by 
ICCAT. 

Other fisheries scientists cooperated in these courses and they are also co-responsible for the 
orientation of the subjects studied and very particularly for the elaboration of the exercises and 
the editorial work. With no particular criterium, these are some of the collaborators to whom 1 
express my appreciation: Ana Maria Caramelo, Manuel Afonso Dias, Pedro Conte de Barros, 
Manuela Azevedo Lebre, Raul Coyula, Renato Guevara. 

Lisbon, December 1997 
E. Cadima 
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GLOSSARY OF TECHNICAL TERMS USED IN THE MANUAL 


Abundance index (U) - A characteristic preferably proportional to the available biomass of the 
resource. The catch per unit effort, cpue (especially when the effort is expressed in 
appropriate units) is an important index.. 

Biological Limit Reference Point (LRP) - Biological reference point indicating limits of the 
fishery exploitation with regard to stock self-reproduction, aiming at conservation of the 
resource. 

Biological Precautionary Reference Point (PaRP) - biomass levels (Bpa) and fishing levels 
(Fpa), established under the precautionary principle, concerning the reproduction of the 
stock, aiming at conservation of the resources. The assumptions and methods used to 
determine the PaRPs should be mentioned. 

Biological Reference Point (BRP) - Values of F and B, taking into consideration the best 
possible catch and/or ensuring the conservation of the fishery resource. There are BRPs 
based on long term projections (LP), BRPs based on values observed during a certain 
period of years and BRPs based on the two previous criteria. The BRPs can be Target- 
Points (TRP), Limit-Points (LRP), and Precautionary Points (PaRP). In this manual the 
following biological reference points are referred to: F m „, Fo.i, Fhi g h, F raed , Fmsy, Fi„„, 
Fcra*h- B ma „ Ba.i, B ml ,d, B>i S y, B| 0 „, MBAL. Other biological reference points, used in 
management, like Fjoy.spr, are not mentioned in this manual. 

Biological Target Reference Point (TRP) - Biological reference point indicating long term 
objectives (or targets), for the management of a fishery, taking into consideration the best 
possible catch and ensuring the conservation of the stock. 

Biomass (B) - Weight of an individual or a group of individuals contemporaneous of a stock. 

Capturability Coefficient (q) - Fraction of the biomass that is caught by unit of fishing effort. 

Carrying capacity (k) - Capacity of the environment to maintain the stock living in it. It is, 
theoretically, the limit of the non exploited biomass (see intrinsic gross rate of the 
biomass, r). 

Catch in number (C) - Number of individuals caught. 

Catch in weight or Y'ield (Y) - Biomass of the stock taken by fishing. Yield does not 
necessarily correspond to landed weight. The difference between the two values, yield 
and landings, is mainly due to rejections to the sea of part of the catch which, for some 
reason (price, quality, space problems or even legal reasons), is not landed. 

Cohort - Set of individuals of a fishery resource bom from the same spawning. 

Exploitation pattern of a gear (s) - Fraction of the individuals of a given size, available to the 
gear, which is caught. Also designated by Selectivity or partial recruitment. 

Individual growth coefficient (K) - Instantaneous relative rate of change of a function of the 
individual weight, w, that is, H(w„)-H(w), where vv_ is the asymptotic individual weight 
and H(w) is a function of w (frequently a power function, including the logarithmic 
function). The adopted models for the function H(w) have two constants, w„ and K. 
Some models introduce one more parameter, b, which is used to obtain a general relation 
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to include the most common individual growth relations. The constant K has the physical 
dimension of time 1 . 

Individual Quota (IQ) - Quota attributed to a vessel. 

Individual Transferable Quotas (ITQ) - System of fisheries management characterized by the 
sale, at auctions, of the fishing annual vessel quotas. 

Minimum Biomass Acceptable Level (MBAL) - Biological reference limit point that indicates 
a spawning biomass level under which the observed biomasses during a period of years, 
are small and the associated recruitments are smaller than the mean or median 
recruitment. 

Number of individuals of a cohort or of a stock (N) - Number of survivals of a cohort (or a 
stock) at a certain instant or over an interval of time. 

Partial recruitment - (see exploitation pattern) 

Precautionary principle - This principle establishes that a lack of information does not justify 
the absence of management measures. On the contrary, management measures should be 
established in order to maintain the conservation of the resources. The assumptions and 
methods used for the determination of the scientific basis of the management should be 
presented. 

Production models - Models that consider the biomass of the stock as a whole, that is, they do 
not take into consideration the age or size structure of the stock. These models are only 
applied in analyses that consider fishing level changes, as they do not allow the analysis 
of the effects of changes in the exploitation pattern, on catches and biomasses. 

Quota (Q) - Each of the fractions in which the TAC was divided. 


RATES 

Absolute Instantaneous Rate of y, air(y) - Velocity of the variation of the function y(x), at the 
instant x. 

Absolute Mean Rate of y, amr(y) - Mean velocity of the variation of the function y(x), during 
a certain interval of x. 

Annual Survival Rate (S) - Mean rate of survivals of a cohort during one year, relative to the 
initial number. 

Exploitation Rate (E) - Ratio between the number of individuals caught and the total number of 
individuals dead, over a certain period of time, that is, E = C/D. 

Fishing mortality instantaneous rate (F) (Fishing mortality 1 coefficient) - Relative 
instantaneous rate of the mortality of the number of individuals that die due to fishing . 

Intrinsic rate of the biomass growth (r) - Constant of the Production models that represents 
the instantaneous rate of the decreasing of the function H(K)-H(B), where B is the 
biomass, H(B) is a function of the total biomass, usually a power-function, (including the 
logarithmic function that can be considered a limit power function) and k is the carrying 
capacity of the environment. Some models introduce one more parameter, p, which is 
used to obtain a more general relation. 
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Natural mortality instantaneous rate (M) (Natural mortality Coefficient) - Instantaneous 
relative rate of the mortality of the number of individuals that die due to all causes other 
than fishing. 

Relative instantaneous rate of y, rir(y) - Velocity of the variation of the function y(x), relative 
to the value of y, at the instant x. 

Relative mean rate of y, rmr(y) - Mean velocity of the variation of the function y(x) relative to 
a value of y, during a certain interval of x. 

Total mortality instantaneous rate (Z) (Total mortality coefficient) - Relative instantaneous 
rate of the mortality of the number of individuals that die due to all causes. Z, F and M 
are related by the following expression : Z=F+M . 

Recruitment to the exploitable phase (R) - Number of individuals of a stock that enter the 
fishery area for the first time each year. 

Selectivity - (see exploitation pattern) 

Spawning or adult hiomassa (SP) - Biomasss of the stock (or of a cohort) which has already 
spawned at least once. 

Stock - Set of survivals of the cohorts of a fishery resource, at a certain instant or period of time. 
It may concern the biomass or the number of individuals. 

Stock-Recruitment (S-R) relation Relation between the parental stock ( spawning biomass) 
and the resulting recruitment (usually the number of recruits to the exploitable phase). 
The models have two constants, a and k. The constant k has the physical dimension of 
weight and a has the dimension of weight' 1 . Some models introduce one more 
parameter, c, which is used to obtain a general relation that includes the most common 
relations. 

Structural models - Models that consider the structure of the stock by ages or sizes. These 
models allow one to analyse the effects on catches and biomasses, due to changes in the 
fishing level and exploitation pattern. 

Total Allowable Catch (TAC) - Management measure that limits the total annual catch of a 
fishery resource, aiming to indirectly limit the fishing mortality. The TAC can be divided 
into Quotas (Q) using different criteria, like countries, regions, fleets or vessels. 

Total number of deaths (D) - Total number of individuals that die during a certain period of 
time- 

Virgin biomass (VB) - Biomass of the stock not yet exploited. 
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SYMBOLS 


Symbols Indicating : 


A 

a 

amr(y) 

air(y) 

B 

B 

SP, SB 

C 

C 

D 

E 

F 

Cconst 

Cte 

H 

ITQ 

K 

k 

k 

L,L 

MBAL 

M 

N 

P 

Q 

R 

R 2 

R 

rmr(y) 

rir(y) 


Constant of the simple linear model (intercept of the straight line ) 

Constant of the Stock-Recruitment relations (limit value of R/S when S— >0) 

Absolute mean rate of variation of y 

Absolute instantaneous rate of variation of y 

Constant of the simple linear model (slope of the straight line) 

Biomass 

Spawning Biomass 
Catch, in number 

Constant of the Stock-Recruitment relations (generalizes the models) 

Total number of deaths 
Exploitation rate 
Fishing mortality coefficient 
Non defined constant 
Non defined constant 
General power function 
Individual Transferable Quotas 

Constant of the individual growth models (associated to growth rate) 
Constant of the Stock-Recruitment relations 
Constant of the production models (Carrying capacity) 

Total length of an individual 

Minimum Biomass Acceptable Level (biological reference limit point) 
Natural mortality coefficient 

Number of individuals of a cohort 

Constant of the Production models (generalizes the models) 

Capturability coefficient 

Constant of the Production models (intrinsic rate associated with the biomass 
growth) 

Determination coefficient 
Recruitment to the exploitable phase 
Relative mean rate of variation of y 
Relative instantaneous rate of variation of y 
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Symbols Indicating : 

S Annual Survival rate 

S Adult or total biomass (in the relations S-R) 

s Exploitation pattern (selectivity) 

SQ Sum of the squares of the deviations 

S-R Stock-Recruitment relation 

t Instant of time 

T Interval of time between 2 instants 

TAC Total Allowed Catch 

TRP Biological Target Reference Point 

l) Stock abundance index 

V Function to be maximized for the determination of Fo.i 

W Individual weight 

Y Catch in weight 

Z Total mortality coefficient (total mortality instantaneous rate) 


xv 


Copyrighted material 



SUBSCRIPTS 


The characteristics of this glossary are usually shown with indices; that is why it was considered 
necessary to present the meaning of those subscripts. 

Subscripts Indicating : 


§ Economical value of the respective characteristic of the cohort 

X Maximum age 

0.1 Value of F (and of other characteristics of the cohort) corresponding to the 

air of the biomass equal to 10 percent of the virgin biomass 

c Recruitment to exploitable phase 

crash Value of F which, at long term, corresponds to the collapse value of the 

spawning biomass 

E Value of the characteristics of the cohort corresponding to an equilibrium 

point 

i Age 

infl Value of the characteristic corresponding to an inflection point of any 

relation between that characteristic and other variable. 

I Length 

lim Value of B or of F corresponding to a biological reference limit point 

loss Value of B or of F corresponding to the minimum spawning biomass 

observed 

Max Value of F (and of other characteristics of the cohort) where the yield per 

recruit is maximum 

Med Value of F (and of other characteristics of the cohort) which, at long term, 

will produce a spawning biomass per recruit equal to the median value of the 
spawning biomasses per recruit observed during a certain period of years 

MSY Value of F (and of other characteristics of the stock) where the long term 

total yield is maximum 

R Recruitment to the exploitable phase 
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CHAPTER 1 - INTRODUCTION 


1.1 THE IMPORTANCE OF FISHERIES 

The importance of fisheries in a country cannot only be measured by the contribution to 
the GDP, but one must also take into consideration that fisheries resources and products 
are fundamental components of human feeding and employment. 

Another aspect that makes fisheries resources important is the self renewable character. 
Unlike mineral resources, if the fishery resources or any other biological resources are 
well managed, their duration is pratically unlimited. 

An important conclusion is that the fundamental basis for the conservation and 
management of fisheries resources stems from the biological characteristics. (This does 
not mean that social, economic or any other effects are not important for management). 

In Portugal, the fisheries contribution to the GDP is less than 1.5 percent. However, with 
regard to food, the annual consumption value of 60 kg of fish per person, is very high. 
Only countries like Iceland, Japan and some small insular nations reach a higher value. 
We still have to consider that of the total amount of protein necessary in our food 
consumption, 40 percent comes from fisheries. This corresponds to 15 percent of the total 
amount spent on food by the Portuguese population. 

From a social point of view, we estimate that there arc, at present, 34 000 fishermen in 
Portugal. Assuming that each job at sea generates 4 or 5 jobs on land (canning, freezing 
and fish meal industry, commercialization, administration, research and training, etc.) one 
can estimate that about 150 000 Portuguese work in the several sectors of fisheries. 
Consequently, taking a minimum of 3 people per family, it is not unreasonable to say that 
about half a million Portuguese people depend on fisheries activities for their livelihoods. 


1.2 FISHERIES RESOURCES MANAGEMENT 

Saetersdal ( 1 984) defined a general principle of fisheries management as : 

“to obtain the BEST POSSIBLE utilization of the resource 
for the benefit of the COMMUNITY” 

It will be necessary to define, in each particular case, what best, possible and community 
mean. 
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In fact, best can be taken as : 


• Bigger yield 

• Bigger value of the catch 

• Bigger profit (difference between the value of the landing and the costs of 
exploitation) 

• More foreign currency 

• More jobs, etc. 

Community may also be taken as : 

• The population of the world 

• The European Community 

• A country 

• A region 

• Groups of interests ( fishermen, shipowners, consumers, . . . ) 


Possible 

reminds us that we cannot forget the self renewable character of fisheries resources and 
consequently, that the conservation of the fishery resource must be guaranteed in order to 
allow the application of the general principle for a long period of time. This statement 
means that conservation of an ecosystem does not imply that one should attribute the 
same importance to all its components. 


1.3 FISHERIES RESOURCES RESEARCH 

Figure 1.1 shows that the research on fisheries resources covers several sectors of the 
fishing activity The assessment models are the main concern of this manual. Among the 
several works and books on fish stock assessment, the books and/or manuals by Beverton 
& Holt (1956), Ricker (1958, 1975) and Gulland (1969, 1983) are historical standing 
references, 
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FISHERIES RESEARCH 


Gears 

Transformation 
Legislation 
Economical aspects 


Etc... 



Figure I.I The several sectors of fishing activity 


1 .4 FISH STOCK ASSESSMENT 

The following are necessary to assess a fish stock: 

• The appropriate data bases 

• Analyses of the available data 
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• Short and long-term projections of the yield and biomass 

• To determine long-term biological reference points 

• To estimate the short and long-term effects on yield and biomass of different 
strategies of the fishery exploitation 

The different steps to assess a stock can be summarized as follows : 

a) To define the objectives of the assessment according to the development phase of 
the fisheries and the available information. 


b) To promote the collection of information : 

• Fisheries commercial statistics : total and by resource landings, catch per effort, 
fishing effort (number of trips, days, tows, time spent fishing, etc.), and 
characteristics of the gears used. 

• Types of operation of the fleets and of its fishing gears, etc. 

• Biological sampling in the landing ports. 

• Biological sampling (and information about the fishing operation) on board 
commercial vessels. 

• Biological sampling on board research vessels. 


c) To analyse the stocks 

The knowledge gained about the resource and the available basic data, determine the type 
of models that should be used and consequently the type of analyses that can be done. As 
an illustration, let us look at some general situations : 
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Fishery resource with little information 


Analyses using particular methods to estimate biomasses and potential yields. 


Fishery resource with data on catches and catch per effort (CPUE) or stock abundance 
indices during several years 


V 

Analyses using production models in order to make projections of yield and catch per 
effort. 


Fishery resource with information collected over several years on : 

• Biological distribution of the catches by species, by length, by ages, etc. 

• Commercial catches 

• Fishing effort or CPUE 

• Research cruises (distribution of the stock by areas, by length, by ages, etc.) 

V 

Historical analysis of the stock (VP A) 

Long and short-term projections with different conditions (scenarios) 


Comments 

1. The lack of information may prevent certain projections, but allows other types of 
analyses. 

2. The Precautionary Principle forces one to estimate and to project catches and 
biomasses, even if they are not very precise. This will be discussed later. 
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CHAPTER 2 - MODELS AND RATES 


2.1 MODELS 

Science builds models or theories to explain phenomena. One observes phenomena and then 
looks for relations, causes and effects. Observations are made about the evolution of a magnitude 
(characteristics) with time (or with other characteristics) and possible causes (factors) are 
explored. Examples: 

• Physics - phenomenum of the movement of the bodies (characteristics - distance related to 
time spent) 

• Biology - phenomenum of growth (characteristic - length or weight, related to time). 


2.1.1 STRUCTURE OF A MODEL 
Basic assumptions 

The assumptions to serve as a basis for a model should : 

• simplify reality 

• be simple and mathematically treatable (manageable) 

• not be contradictory 

• not be demonstrated 

• be established with the characteristics 

Usually basic assumptions are related to the evolution of the characteristics. So, they are 
established on the variation rate of those characteristics and they do not need to be proved. 

Relations (properties) 

• they are deduced from the basic assumptions by the laws of logic (mathematics). 
The properties are also designated by 


“results” or “conclusions” of the model. 


Verification 

• the results of the model must be coherent (to agree) with reality. 

This implies the application of statistical methods and sampling techniques 
to check the agreement of the results with the observations. 
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Improvement 

• if agreement is approximate, it is necessary to see if the approximation is enough or not. 

• if the results do not agree with reality, then the basic assumptions have to be changed 

• the changes can aim to the application of the model to other cases. 

Advantages 

• it is easier to analyse the properties of the model than the reality. 

• the models produce useful results. 

• they allow analysis of different situations or scenarios by changing values of the factors. 

• to point out the essential parts of the phenomenon and its causes. 

• they can be improved in order to adjust better to the reality. 


2.1.2 SOME TYPES OF MODELS USED IN STOCK ASSESSMENT 
Production Models 

The production models are also designated as General Production models. Global models. 
Synthetical models or Lotka-Volterra type models. These models consider the stock globally, in 
particular the total abundance (in weight or in number) and study its evolution, the relation with 
the fishing effort, etc.. They do not consider the structure of the stock by age or by size. 

Structural Models 

These models consider the structure of the stock by age and the evolution of the structure with 
time. They mainly recognize that the stock is composed of individuals of different cohorts, and, 
consequently, of different ages and sizes. So, they analyse and they project the stock and the 
catches for the coming years, by following the evolution of its different cohorts. 

This manual will not follow the chronological construction of the models. It was thought to be 
more convenient to deal firstly with the structural models and afterwards with the production 
models. 
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RATES 


■ basic assumptions of a model, for the evolution of a characteristic, require the concept of 
■ation rate of the characteristic related to time (or to other characteristics). 



Figure 2.1 Evolution of the length (L) of an individual with time (or age) (t) 

In order to generalize the study of the rates, the characteristic L in the example above will be 
substituted by y, and the associated variable will not be time, t, but the variable x. To study the 
stock assessment models and to make this study easier, it will be considered that the function y 
will only assume real and positive values. 


2.2.1 ABSOLUTE MEAN RATE - amr (y) 

Consider y a function of x and the interval i with the limits (x„ Xj+i) 



Figure 2.2 Function y= f(x) with variation in the interval i 

Let : 

Axj = x j+i - x, be the size of the interval 
yi = the value of y when x = xt 
yi-H = the value of y when x = Xt+i 
The variation of y in the interval Ax, will be Ay, = y hi - yi 

The absolute mean rate, amr (y), of the variation of y within the interval Ax,, will be : 
amr(y)= 

Ax, 
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Graphically : 


Figure 2.3 



Absolute mean rate of the variation of v within the interval Axi 


Ay 

Slope of the secant = — - = amr (y) during Axj 
Ax, 

Note: amr (y) is known in physics as the mean velocity of the variation of y with x, in the interval 
Ax,. 


2.2.2 ABSOLUTE INSTANTANEOUS RATE - air (y) 

Let y be a function of x 

The absolute instantaneous rate of y at the point x = Xj is the derivative of y in order to x at that 

point 


air (y),.„, 



Graphically : 

y 


■+- 

xi 


x 



air(y)„ = , equal to the slope of the 
tangent to the curve at the point x-xi 


Figure 2.4 Absolute instantaneous rate of v at point xi 

Note: air(y) is known as the instantaneous velocity of variation of y with x at the point x. 
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Properties 

1 . Given the value of air(y) the calculation of the function y is obtained, by integration, being 
y = f(x) + Constant, where f(x) = Primitive of air (y) and Constant is the constant of 
integration. 

If the initial condition x* , y* is adopted, where y* is the value of y corresponding to x = x*, 
eliminating the Constant, then one can write y = y* + f(x)-f(x*) 

2. The angle made by the tangent to the curve y with the xx's axis is designated by inclination. 

The trigonometric tangent of the inclination is the slope of the geometrical tangent, 
air (y) = derivative of y = slope = tg (inclination) 

3. If, at point x : 

air (y) > 0 then y is increasing at that point 
air (y) < 0 then y is decreasing at that point 
air (y) = 0 then y is stationary at that point (maximum or minimum) 

4. If air (y) is constant (= const) then y is a linear function. From property 1, it will be : 

y = Constant + const, x or 
y = y* + const.(x-x*) and vice versa 


5. If y(x) = u(x) + v(x) then air (y) = air(u) + air(v) 

6. If factors A and B cause variations in y, then factors A and B considered simultaneously 
cause a variation of y with: 

air (y) ltJ11 ii - air (y) cmka + air (y) cimb 

d 2 y 

air(air(y)) = — y-- acceleration of y at the point x 
d‘x 

7. If the acceleration at the point x is increasing, then air (y) is positive and if that acceleration 
is decreasing, then air (y) will be negative. 


2.2.3 RELATIVE MEAN RATE - rmr (y) 

Consider y a function of x and the interval (x*, x j+t) 

Let : 

Axj = x s-H - Xj = the size of the interval 
yi = value of y when x = xi 
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y, + i = value of y when x = xm 

Xi* = a certain point in the interval (x„ x ,»i) 

y,* = value of y when x = x,* 

Xj* can be xi , x i+ i , x t , etc. 


The mean rate ofy relative toy f will be : 


nnr(y)= — or 

y, Ax, 

rmr(y) = — .amr(y) 


Comments 

1 . rmr (y) is associated with the mean rate of the variation of the percentage of y in relation to 
y*, that is: 


A(4) 

y, 

Ax* 


„ , , Ax, 1 , „ _ 

2- Let x OTMl be x^ =x i +-j‘- = -.(x i +x w ) = x i 

3. It is convenient to designate by y^^, the value of y in the interval (Xj , Xj»i) when 

X = X central, • 

fy + y ) 

Notice that y can be different from the mean, 3 ' 1 

4. It is frequent to calculate rmr(y) in relation to y ccmnl of the interval. 


2.2.4 RELATIVE INSTANTANEOUS RATE - rir(y) 


Let y be a function of x. 

The relative instantaneous rate of y at the point x = x, is 


rir(y) = — 

1 



or 


rir(y) = 


air(y)„, 

y, 
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Properties 

1 . Given rirfy), the calculation of the function y is obtained by integration, being 


y = f(x) + Constant, where f(x) = Primitive of rir(y) and C is the constant of integration. 

If one accepts the initial condition x* , y* , where y* is the value of y corresponding to 
x = x*, one will get, eliminating the Constant, y = y* + f(x) - f(x*) 


2. If, at a point x : 

rir (y) > 0 then y is increasing at that point 

rir (y) < 0 then y is decreasing at that point 

rir (y) = 0 then y is stationary at that point (maximum or minimum) 

3. rir(y) = air (lny) as can be deduced from the derivation rules. 

4. If rir (y) = constant = (const) then y is an exponential function of x, that is, 

y = Constant . e “ nMJl or 

y = y* e and vice-versa 


5. If y(x) = u(x).v(x) then rir(y) = rir(u) + rir(v) 

6. If the factors A and B cause variations in y, then simultaneously, 
factors A and B cause a variation in y, with: 

rir (y),„ u i = rir(y) cause A + rir (y), ausc B 

2.3 SIMPLE LINEAR MODEL 

Let y = f(x) 

Basic assumption of the mode i 

air(y) = Constant = b in the interval (x„ Xj+i) with AXi = Xj+i-Xi 

Initial Condition 

x* = Xj => y* = yi 



Figure 2.5 Graphical representation of a simple linear model 
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Properties 


1 . General expression 

2. Value yi+i at the end of the 
interval. Ax, 

3. Variation, Ayj , in the interval, 
Axi 

4. Central value, y c=1M1 
of the interval. Ax, 

5. Cumulative value, y ctlm 
during the interval. Ax, 


6. Mean value, y, , in the 
interval. Ax, 

Other useful expressions 

7. Cumulative value, y^ 
during the interval, Axj 

8. Mean value, y, , during the 
interval. Ax* 

9. Mean value, y, , in the 
interval. Ax, 

10. Mean value, y, , during the 
interval. Ax, 

1 1 . Relation between amr(y) et 
air(y) 


y = y, + b-(x-Xj) ; y = a + bx 
y. + . =y, +b-Ax, 

Ay, = y,„ -y, = b.Ax, 

y^, = y. + b -( x c«-i -x,) = y,+bY 

y««,, = Jydx =Ax, (a + bx,) 
or from the Property 1 
y™™, = Ax l [y,+b(x,-x j )] 

y, _ >W _ a + k x where 

' Ax, 

y,=v^ = y, + b (x,-x i ) 

Ax, 

y.«, = Ax, -y, 

y, =y, +b.(x, -x,) where y,=a + bx. 



y i “ y «nml, 

amr (y , ) = ~ L = b = air(y ) 
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12. If Ay, < 0 then b < 0 et vice-versa 

13. In the linear model, the arithmetic mean of y, and y„, is equal to the mean value, y, , 
and equal to the central value y callrii; 


Important demonstrations 


General expression If tia(y) = b in the interval Ax, then y is linear with x and considering 

Property 1 the initial condition it will be: y = y;+ b.(x-Xj) 


Central value 
Property 4 


y«™i, -y, + b.(x OT „, - )= y, + b^(x, + ^ ) - X,j = y, + b.^- 


Cumulative value 
Property 5 


from the definition of the cumulative value : 
y cum, = f'"(a + bx)dx 


= a (*i.i 



it will be necessary to use the factorization of the difference of two 
squares, that is : 

xVi - X 2 = (Xj+I - Xj).(x i+ | + Xi) = Ax ( . (x i+ i + Xj) 


y, et y 

Property 10 


and then: 

y„n„ = aAXi + bAx, . Xi = Ax, (a + b. x,) 


y, = y, + b.^-x, )= y, + b^Y-y) = y. + b.^- = 


2.4 EXPONENTIAL MODEL 

Let y = f(x) 

Basic assumption of the mode! 

rir(y)=Constant=c in the interval (xi, Xj+i), with Axj = Xj+i - Xi 

Initial condition 

x* = Xi => y* = yi 
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Properties 

rir(y) ~ air(lny) means (hat the exponential model of y against x is equivalent to the linear model 
of Iny against x. So being, its properties can be deduced by backwards application of logarithm 
rules to the properties of the linear model of Iny against x. 



Figure 2.6 Graphical representation of the exponential model 



Figure 2.7 Graphical representation of the linear model of Iny against x 


1 . 

2 . 

3. 

4. 


General expression 

Value of yni at the end of the 
interval. Ax, 

Variation, Ay, , during the interval, 
Axi 

Central value, y„„ M , , in the interval 
Axi 


Exponential model of y 

Linear model Iny 

(y against x) 4 

_flny against x) 

y = y, •e cU ‘ Xi) ^_ 

Iny = lny,+c(x-x,) 

y..,=y,-e7 

Iny it |= lnyi+cAx, 

Ay, = y M -y, =y i (e , '“’ 

-l)calculated from 1 

c.A>, 

. . cAx 

y =y;-e * 

- toy,*..* =tay , +“ 

y«nmn, = (yi-yHt)', : 

- tay ranV = (lny,+lny, +1 )/2 

(y^ = geometric 


mean of the extremes y, 


and y i+ i) 
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Cumulative value, y (<m _ , during the 
interval, Ax, 

yen,, = Jy.dx=^ 

X, c 


Mean value, y, , during the interval, 

_ _ y t »», _ 1 Ay, 


Axi 

Ax, c Ax, 



e^'-l 

y i“ y >- A 

cAx, 

(replacing Ay, using 
Propriety 3) 


y- y.*.~y. 

In yi», - In y. 


y, = y cenmi, 


Other useful expressions 

7. Expressions of variation, Ay, 

8. Expression of amr (y) 

9. Expression ofrmr (y) 
in relation to y* 

10. Expression of rmr (y) 


11 . 

If Ay, < 0 alors 


Ay, = c. y^, 

Ay, = c.y.Ax, 

Ay _ 

amr(y) = ~r~ L = cy. 

Ax, 

rmr (y)^ l .,,o„„^ =c=tir (y) 

rmr(y) = amr (Iny) = ^-= c 

Ax, 


y decreases 
c<0 

amr(y)< 0 

rmr(y) < 0 and vice-versa 

y««, > o' 

y. >0 , 
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12. In the exponential model, the geometric mean of y ( and y,., is equal to the central value, 
yctmni, (Prop. 4) and approximately equal to the mean value, y, (Prop. 6), been the 
approximation better when Ax, is smaller. 


Demonstrations 


Cumulative value 
Property 5 


y 


cum, 


= [‘•'y-dx= [ , "y i -e ,< "-*‘>dx = 




Relation between 

and 


Property 6 - 
4 ,h expression 


From the approximation * Is- e h,2 with h = c.AXj 

and from property 6-2" d expression, will be: 

y; = yi 

Finally, by property 4-1“ expression, one can conclude that: 

* = y central, 
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CHAPTER 3 - COHORT 


3.1 COHORT - INTRODUCTION 

A cohort or annual class or a generation, is a group of individuals bom in the same spawning 
season. The following scheme illustrates the different phases of the life cycle of a cohort: 


PRE-EXPLOIT ABLE PHASE 


egg larvae Recruitment 

I 

juvenile 

I 


EXPLOITABLE PHASE 


Adults 

1st spawning | 2nd spawning 


-> 


\l/ ^ ^ 


new cohorts 


Figure 3.1 Cycle of life of a cohort 

Let us start, for example, with the egg phase. The phases that follow will be larvae, juvenile 
and adult. 


The number of individuals that arrive in the fishing area for the first time is called recruitment 
to the exploitable phase. These individuals grow, spawn (once or several times) and die. 

After the first spawning the individuals of the cohort are called adults and in general, they 
will spawn again every year, generating new cohorts. 

The phases of life of each cohort which precede the recruitment to the fishing area (egg, 
larvae, pre-recruits), are important phases of its life cycle but, during this time they are not 
usually subjected to exploitation. The variations in their abundances are mainly due to 
predation and environmental factors (winds, currents, temperature, salinity,...). In these non 
exploitable phases mortality is usually very high, particularly at the end of the larvae phase 
(Cushing, 1996). This results in a small percentage of survivors until the recruitment. Notice 
that this mortality is not directly caused by fishing. 

The recruitment of a cohort during the exploitable phase, may occur during several months in 
the following schematic ways : 


number of recruits 



1 year 1 year 1 year 

Figure 3.2 Types of annual recruitment to the exploitable phase 
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With some exceptions, the forms of recruitment can be simplified by considering that all the 
individuals are recruited at a certain instant, t r called age of recruitment to the exploitable 
phase. It was established that recruitments will occur on 1 January (beginning of the year in 
many countries). These two considerations do not usually change the results of the analyses, 
but simplify them and agree with the periods of time to which commercial statistics are 
referred. 

It should be mentioned that not all the individuals of the cohort spawn for the first time at the 
same age. The proportion of individuals which spawn increases with age, from 0 to 
100 percent. After the age at which 100 percent of the individuals spawned for the l sl time, all 
the individuals will be adult. The histogram or curve that represents these proportions is 
called maturity ogive. 

In certain cases, the mamrity ogive can also be simplified supposing that the P’ spawning 
occurs at the age t Mt designated as age of 1 st maturity. This simplification means that the 
individuals with an age inferior to W are considered juveniles and those with the same age or 
older, are considered adults. 

Figure 3.3 represents a maturity ogive with the shape of a histogram or curve : 



Figure 3.3 Maturity ogive 


3.2 EVOLUTION OF THE NUMBER OF A COHORT, 

IN AN INTERV AL OF TIME 

Consider the interval (tj, tj,|) with the size T, = t,,| - 1 , of the evolution of a cohort with time 
and Nt the number of survivors of the cohort at the instant t in the interval Ti (see 
Figure 3.4). 

The available information suggests that the mean rates of perccntual variation of N, can be 
considered approximately constant, that is, rmr (N.) = constant. 
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Basic assumption 

The relative instantaneous rate of variation of N, , in the interval T, is : 
rir (N t ) = constant negative - - Z, 



Ni = number of survivors 
at the beginning of the 
interval (tj, ti +1) 

Nj+i = number of survivors 
at the end of the interval 
(ti, ti +1) 


Figure 3.4 Evolution of N in the interval Ti 

The model of the evolution of Nt, in the interval Ti, is an exponential model (because rir(N, is 
constant). This model has the following properties: 

Properties 

1 . General expression. From the basic assumption 

riifNt) = -Zj 

with the initial condition that, for t = ti it will be N t = Ni then: 

N, = N, • e~ z<, ' ,,) 

2. Number of survivors, N,+i , at the end of interval Ti 

N nl = Nj -e ZT ' 

3. Number of deaths, Dj , during the interval Ti 

d, = n,-n w 

D, =N,(l-e' z ' ,T ‘) 

(notice that Di is positive but the variation AN, = Nj+i- Nj is negative) 
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4. Cumulative number of survivors, N cum , , during the interval Ti 


N =5l. 

Z 


N = N, 


1 — e - ‘ 


5. Approximate central value , N cmua i, , in the interval Tj 


Nccmnii 


= N. e 


-Z.T. 




6. Mean number, ^i, of survivors during the interval Tj 


— N 

JyJ. = 


N, = N, 


l-e^ T| 
1 Z.-T, 


Ni =- Dl 


Z; -T, 

N,=- N,_N “ 


N, =N, -e 


InN, -lnN h 


(Ricker) 


N, = N camli when is small (Zj.Tj < 1) 

Comments 

1 . The basic assumption is sometimes presented in terms of absolute instantaneous rates, that 
is: 

air (N t ) = - Z, .N t { air (N.) proportional to N ( } or 
air (InNt) = -Zj 

Z, = mortality total coefficient, assumed constant at the interval Tj 


Notice that: 

+Z, = rir of total mortality of N t 
-Zj - rir of variation of Ni 

2. Unit of Zj 
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From the definition, it can be deduced that Z, is expressed in units of [time]' * By 
agreement, the unit year 1 has been adopted, even when the interval of time is smaller or 
bigger than a year. 


The following expressions show, in a simplified way, the calculation of the unit of Zi, 
with the rules and usual symbols [..) of dimension in the determination of physical units . 


JlL 

[N,] 


|dN , ] 
[dt] 


= -Z. 


1 number r _ , . r _ t 

= +IZ, j then |Z, J = time 

number time 


3. Annual survival rate, S, 
When Tj =1 year, it will be : 


N , =N, = ^ 


N lt ,=N l e' 2 > 


and also 


Sj = Annual survival rate in the year i 

(or percentage of the initial number of individuals that survived at the end of the year). 


S, =e" Z| 

1 - S, = Annual mortality rate in the year i 

The percentage of the initial number of individuals that die during the year is, by 
definition, the relative mean rate rmr (N t ) of mortality of N, , over one year, in relation to 
the initial number, Ni 


1-S = -51 = 1-21^. 
Nj N, 

I — S, =l-e' z ' 

4. Absolute mean rate 

5. Relative mean rate 


N, ~ N„, 

N, 

amr(N,) = -Z^.N, 

rmr( N , ) = -Z, in relation to N , 
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6. Notice that Sj takes values between 0 and 1 , that is: 

0 < S, < 1 but Zj can be > 1 

7. If the limits of the interval T, were (t„ °°) then it would be: 

Tj = oo 

N i+ i=0 


Dj = Nj 

= y and 

“i 

N, =0 

3.3 CATCH, IN NUMBER, OVER AN INTERVAL OF TIME 

The causes of death of the individuals of the cohort due to fishing will be separated from all 
other causes of death. These other causes are grouped together as one cause designated as 
natural mortality. So, from the properties of the exponential model, the result will be 


rir(N t ) ll>Ia i = rir(N t ) natural + rirfN,) rBhlBg 


Supposing that, in the interval T„ the instantaneous rates of mortality due to natural causes 
and due to fishing arc constant and equal to M, and to Fj, respectively, then 

Z = F + M 

I 1 1 

Multiplying both factors of the previous equation (equality) by.N tul „ ( , then: 

Z.'N^ = Fj , N aal| +M| 

Zj-N cum is the number, Dj , of deaths due to total mortality, 

D,=Z,.N„„ 

In the same way F.-N^ will be the number of dead individuals due to fishing, that is, the 
Catch, Q, in number, and then: 

C* = F,.N cumi 

Notice too that, M <■, will be the number of dead individuals due to “natural” causes. 

The exploitation rate, E, , during the interval T, was defined by Beverton and Holt ( 1 956) as: 
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jj _ number capture _ C, 
number dead D, 

F.N F 

and then, E. = — = — or 

Z i- N a,m, Z, 

Ej =- S_ 

F i + M i 

The capture in number, Cj, in the interval T„ can be expressed in the following different ways: 
C, =F,.N cum) 


C, = F| . Nj.Tj 

Ci = Ej-D, 

F 

C, = — .D, 

Z, 

C = — 5 N,[l-e' (F|+M|>T ‘] 

Fj + M, 


Comments 

1 . Ricker (1975) defines the exploitation rate, Ei*, as the percentage of the initial number 
that is captured in the interval Ti, that is: E,* = C, / Nj. 


a) Ricker’s definition may be more natural, but mathematically Beverton & Holt’s 
definitions are more useful. 

b) It is easy to verify that E' = E,.(l - e~ z, T ‘ ) 

2. The exploitation rate, Ei , does not have any unit, it is an abstract number. 

3. The possible values of Ej are between 0 e 1, being 0 when there is no exploitation and 1 
when the capture C, is equal to the number of total deaths D,, that is, when M, = 0. 

3.4 INDIVIDUAL GROWTH 

In order to study the evolution of the biomass of a cohort, one can use the model of the 
evolution of a cohort, in number, and combine it with a model of the evolution of the mean 
weight of an individual of the cohort. In effect, the biomass B, is equal to N,.W, where W t is 
the individual mean weight at the instant t. 
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To define a model for the individual growth weight W t , there are then two possibilities : 
ALTERNATIVE I: 

A) To define a model for the mean individual growth in length, L, 

B) To define the relation Weight-Length. 

C) To combine A) with B) and obtain a mode for the mean individual growth in 
weight, Wt 

ALTERNATIVE 2: 

D) To define directly growth models for Wt and L t . 

ALTERNATIVE 1 

A) Model for the individual growth by length 

The models that arc used in fisheries biology are valid for the exploitable phase of the 
resource. The most well known is the von Bertalanffy Model (1938) adapted by Beverton 
and Holt (1957). The existing observations suggest that there is an asymptotical length, that 
is, there is a limit to which the individual length tends. 



Figure 3.5 von Bertalanffy Model 


t-age 

L, - individual mean length at the age t 

L» - asymptotical length 

t r - beginning of the exploitable phase 

So, Li presents an evolution where : 

air(Lt) is not constant (because growth is not linear) 

rir(L,) is not constant (because growth is not exponential) 

However, it can be observed that the variation of the quantity (L_ - L t ) (which we could call 
“what is left to grow"), presents a constant relative rate and can be described by an 
exponential model. So, we can adopt the: 
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Basic assumption 

rir (L» - L,) = - K. = negative constant during all the exploited life 

where K is the growth coefficient (attention: the growth coefficient K is not the velocity of 
growth but the relative velocity of what “is left to grow” !!). 

The properties of this model can be obtained directly from the general properties of the 
exponential model. The initial condition : 

t = t. => L, = L. 

where t a (and L a ), that corresponds to an instant within the exploitable phase, will be adopted. 

The properties of the model of individual growth by length by Beverton & Holt (1957) 
deduced from the exponential model of (L_ - L t ) are summarized as : 


Properties of the exponential model 
for(L„-L t ) 

von BertalanfTy Model 
for L t 

1. 

General expression 

(L.-L,)=(L.-L.)-e- K <--> 


L, = L,-(L_-L,)e' M '"- ) 

Parameters 


Parameters 

L,. = asymptotic length 
K. = growth coefficient 


L» = asymptotic length 
K = growth coefficient 

J a = initial condition 

k 


i a = initial condition 

k 

For t, = to and L, = 0 : 


For t, = to and L, = 0 : 

(L_ - L,)= L_ .e' K *'" ,,) 


L, = L_-L_e' K(, "" ) 



L, =L.|l-e- K,, -'" ) J 

2. Value at the end of the interval T» 

(L_-L,J = (L_-L,).e- KT ' 


Lj,i = L_ • (1 - e" KT ' ) + L, • e" KT ' 
Ford-Walford expression (1933-1 946) 
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3. Variation durin 

4 ( L --L i ) = (L 00 -L i+1 )-(L 00 -L i ) 

= (u.-L i )r kT i-(L„-L j ) 

g the interval Tj 

As A(L_ -L,)=-AL { 

It will be:: 

AL,=(L„-L,)(l-e- KT ') 

4. Cumulative value di 

, \ A(L_-L,) AL, 

\ - 

uring the interval Tj 

Ay 

L^, = L -Ti— ± 
from: 

A T 

(u-L i ) a _ i =L„-Ti-u. 11 “— 

5. Mean value duri 

(l — | ) AT 

ng the interval 7) 

L, =L_--^i- 
' K-T, 

Gulland e Holt expression ( 1 959) 
from: (L.-L,)=L_-L, 

L L = i/tum, _ AL, 

T< K-T, 

6. Central value of tl 

-K tccntral _, a 

(L.-LUm, =(L.-L.)e l * 

1 

and 

he interval TV 

I -1 _ tl -IV ' K, lco«nl "I. 1 

and 

^central ~ 


B) Relation Weight-Length 

It is common to use the power function to relate the individual weight to the total (or any 
other) length. Then : 

W, = a • L, b 

where constant a is designated as condition factor and constant b as allometric constant. This 
relation can be justified accepting that the percentage of growth in weight is proportional to 
the percentage of growth in length, otherwise, the individuals would become disproportionate. 
Thus, the basic assumption is : 

rir(W) = b -rir(L) 

where b is the constant of proportionality. 
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C) Combination of A) and B) and comments: 

1- From the combination of W = a • L b with L, = L_ • (1 - e" K1 '"'’ > ) we have 

W, = W. • (1 - e' IC(l ' , * , ) b with W_ = a • L b _ 

This relation of growth in weight is designated as the Richards equation (1959). 

When b=3 the equation is the von Bertalanffy growth equation (1938). 

2. From, W = a L b we have, by definition, W ctmrll ^ = a • , where is the value 

corresponding to 

3 Let w‘ be the weight corresponding to L, , that is, w’ = a • (l, f 
As. L, = then, w> a • = W cct „ k 

In practice, L cenm , and W CCTitrjii are preferred to the mean points 

4. The Richards and von Bertalanffy models are not the only models used for the evolution 
of W,. Other models which have also been used in stock assessment are : 
Gompertz model (1825) and Ricker model (1969). (see Alternative 2 - property 3) 

5. Historically, the von Bertalanffy model was developed from the basic assumption 

2 

tia(W,) = Cte,.W 5 -Cte 2 .W 

Where Ctei and Cte 2 were designated by von Bertalanffy as the anabolism and the 
catabolism constants, respectively. 

Adopting the relation W = a ■ L’ the basic assumption will become 
air(Li) = Ctei - Cte. . L 
(where Ctei e Cte 2 are other constants). 

The solution of this diferential equation gives the von Bertalanffy equation. 

ALTERNATIVE 2 

D) Model directly for W, and L, 

In alternative 1, a model was developed for growth in length, and then a model was built for 
growth in weight, using the relation between weight and length. 

The basic assumption adopted by alternative 1 used the characteristic (L_ - L, ) instead of 
weight in order to be able to have a characteristic with a rate rir constant, that is, an 
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exponential model. The relation W=a.L b was adopted to obtain the model of growth in 
weight. Notice that it can be said that L was considered as a function of W, that is, L= 
(W/a) l/b . It will then be possible to adopt, instead of that function of W, another function of 
the weight H(W,) , in order to be able to formulate directly the basic assumption: 

air[H(W_)~ H(W t )] = -K = constant 

with the initial condition 

t = t, => W, = W, 


Properties 

The properties of this model (once it is an exponential model) can be obtained directly from 
the general properties of the exponential model. It is particularly interesting to derive the 
general expression W, resulting from different choices of function H. 


1. General expression 


[H( W_) - H( W, )] = [H( W_) - H(W, )}e~ K "' B) 


or 


H(w,)= H(W_)-[h(W_)- H(w, )]e‘ M, “- ) 


2a. Richards equation in weight 

The result will be the general 
expression: 


Adopting the following function H (W,) = W, l/b 

I I I I 
W, b = W_ b - (W b - W, b ) • e ' K(l "“’ 


that is, the Richards equation; and when b=3, this is the equation of von Bertalanffy, so: 


2b. von Bertalanffy equation in 
weight, will be : 


■ i i i 

W, 3 * 5 = W l - (Wj - W t > ) • e' K,, " B) 


3. Gompertz equation in weight Adopting the function H(W,) = In W, 

The result will be the general expression: 

InW, = In W„ - (In W_ - In W, ) • e~ K, '~ a) ) 


4. 


The respective equations in length can be obtained by adopting other functions of 
H(W.) : 


4a. von Bertalanffy equation in 
length 


Adopting H(W,) = L, it will be: 
L, = U - (L.-L 1 ).e K(M * ) 


4b. Gompertz equation de Adopting H(W,) = In L, it will be: 

in length 

InL, = lnL„ - (lnL„ -InLJ- e ' K< " M 
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5. Simplified equations 

The individual growth equations, both in length and in weight, are simplified when 
one selects H(W,) = 0 for t a = t* 

So, the simplified general H(W, ) = H(W„ )i(l - e' K * t ' 1 "*) 

expression will be reduced to : 

5a. Simplified Richards equation, w = W • [l - e " lc ^' t, *] b 
in weight, will be : 

where /* was represented by to because H(W a ) = 0 in Richards model means that W> 
will also be zero. 


5b. 


Simplified Gompertz equation, 
in weight, will be: 


InW, = lnW_ 

In this case H(W t ) = 0 corresponds to W a = 1 


5c. Simplified Richards equation 
in length, will be: 


L, = L_ 

(with L„ = 0 for t a = to) 


5d. 


Simplified Gompertz equation, 
in length, will be: 


InL, = tnL„[l-e' K!, ' , ’ ) ] 
(with L„ = 1 for t, = t*) 


Comments 

1 . Gompertz equation , in weight, is similar to Gompertz equation, in length, but, in their 
simplified forms, t* represents different ages, because they will correspond, respectively, 
to W,= 1 and to L,= 1. 

2. Gompertz equation, in length, is similar to von Bertalanffy if L, is substituted by In L,. In 
practice, this fact allows the utilization of the same particular methods to estimate the 
parameters in both equations, using L in the von Bertalanffy expression and InL in the 
Gompertz expression. (See Section 7.4) 

3. It is important to notice, once again that, in practice, L cm iraii and W ccnlra ij arc used instead 
of the mean values, Li and Wi. 

4. Gompertz growth curve in length, has an inflection point, (t,„n, Li„n), with: 

tinfi = t„+( I /K).ln(ln(U7L a )) L ro n = L„ / c 


5. Gompertz growth curve in weight has an inflection point, (ti„n, W; 0 n) with: 

ti„n = t,+(l/K).ln(ln(W„yW a )) W infl = W. / e 

6. Richards growth curve in length does not have an inflection point but the growth curve in 
weight has an inflection point, (turn, W^n). 
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tinfl ti 


= t -— Ini 


vY 


i- 

lw. 


W,„ fl = (l-^-) b .W- 
b 


In the particular case of the von Bertalanffy equation it will be : 

W,„ n = (8/27). W„ and t mfl = to+(l/k).ln3 

7. Some authors refer Gompertz equation in other ways, for example, using the inflection 
point tj„n and the asymptotic weight W,„ instead of the parameters t a and W,. 

It will then be w, = w„.exp(-e _k( '"'"' ) )or L, = L_.exp(- 

Sometimes the length expression is presented in its general form: L t =a.exp(b.e cl ) 

The parameters of the length model will then be: L»= a; k = -c et t mn = (l/c).ln(-b) 

8. The growth in length presents an inflection in fish farming, where the study of growth 
covers very young ages and it is common to use the Gompertz equation. In fisheries, the 
tradition is to use the von Bertalanffy equation. 

9. A model that can sometimes be useful, is the Ricker model (1975). This model is valid 
for a certain interval of time T, and not necessarily for all the exploitable life of the fishery 
resource. In fact, the model is based on the basic assumption that the individual growth is 
exponential in the interval Ti. 

It will be, for example, L t = L i .e K ' 1,1,1 where K, can be different from one interval to the 
next. 

3.5 BIOMASS AND YIELD, DURING THE INTERVAL T, 


1. Biomass 

Theoretically, it could be said that the biomass at the instant t of the interval T, is given by: 


B,= N,. W, 


Thus, the cumulative biomass during the interval T, would be: 

Bcum, = fB.dt 
t. 

and the mean biomass in the interval T would be: 


B, = 


D 

‘-'cum, 

T, 
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In the same way, the mean weight of the cohort, W, , in the interval T, would be: 


— B„ 

W, = — 


The biomass can be obtained by dividing both terms of the fraction by Tj, as 
B, =N|.Wi 


2. Yield 

The yield, Yj , during the interval T, will be expressed as the product of the catch in numbers, 
times the individual mean weight : 

Yj = C, ■ Wi 


Comments 

In practice W, is considered approximately equal to W colOT , at the interval T . 

Other expressions of Yj will also be : 

Y = F N W 
= FiB cum( 

= F,N i .W,T 1 
= F.B.T, 

3.6 COHORT DURING THE EXPLOITABLE LIFE 

Consider the evolution of a cohort during the exploitable life, beginning at age t,, and intervals 
of time, Tj, covering all the exploitable phase (frequently the intervals are of 1 year...). 

Figure 3.6 illustrates the evolution of the number of survivors of the cohort, N, and the 
evolution of the catches in number, C, which are obtained during the successive intervals of 
time Tj. 
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Figure 3.6 Evolution of the number of survivors of the cohort, N, and the catches in 
number, C 

Figure 3.7 illustrates the evolution of the biomass of the cohort, B, and the evolution of the 
catches in weight, Y, which can be obtained during the successive intervals of time T,. 



Figure 3.7 Evolution of the biomass of the cohort, B, and the catches in weight, Y 


Values of the most important characteristics of the cohort, 
during all the exploitable phase 


Duration of the life of the cohort 

A. = £Tj 

Total number of deaths 

D = E D, (= R (recruitment) when all the 
individuals die) 

Cumulative number of survivors 

N C um = £ N CU m j 

Mean number of survivors 

N = N eura /X 
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Cumulative biomass of survivors 

B C um 21 B C umj 

Mean biomass of survivors 

B = B C um / ^ 

Catch in number 

C = £C, 

Catch in weight 

Y = EYi 

Mean weight of the cohort 

W_ = B_/N_ = B/N 

Mean weight of the catch 

W ( ^=Y/C 

Critical age 

Critical = age where B is maximum, when 
the cohort is not exploited. 


Comments 

1. At first sight it may seem that the values of the characteristics of a cohort during all the 
exploitable phase are of little interest, because, very rarely is fishing applied to an 
isolated cohort. At each moment, the survivors of several cohorts are simultaneously 
present and available for fishing. 

2. Despite this fact and for reasons which will be mentioned later, it is important to analyse 
the characteristics of a cohort during all its exploitable life. Knowledge of the evolution of 
a cohort, in number and in biomass, and particularly the critical age, is important for the 
success of the activities in fish farming. As B, = N, . W, , the critical age, t cl , llal i will be 
the age t in the interval T, where 

rir (W,) = - rir(N,)=M because, the critical age being the maximum biomass, the 
derivative of B will be equal to zero. 

3. Notice that N cum can be expressed in function of the recruitment as 

N cum =R.{...} 

where {...} represents a function of biological parameters and annual fishing mortality 
coefficients F, during the life of the cohort. B cum can also be expressed as: 

B CUII! = R .{...} 

where the functionj...} also includes growth parameters. 
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3.7 SIMPLIFICATION OF BEVERTON AND HOLT 

Beverton and Holt (1957) deduced algebraic expressions for the characteristics of a cohort 
during the exploited life, adopting the simple assumptions: 

1. The exploited phase of the cohort is initiated at age lc and is extended to the infinite. 

2. The natural mortality coefficient, M, is constant during all the exploitable phase. 



Figure 3.8 Life of a cohort during the exploitable phase 

3. The fishing mortality coefficient, F, is constant during all the exploited phase. 

4. Growth follows the von Bertalanffy equation with for L„ = 0 for t a = t 0 


Basic equations referring to the exploited phase 

c (ratio between length at age t c and 
the asymptotical length) 


c = — = 1 


2. Recruitment 

3. Cumulative number 

4. Catch in number 


U 

R c = R • 

_ R*. 

Ncum ^ 

C = ER, 


5. Cumulative biomass 


= R W. 


J ■, (l-c) [3 (l-c) ; (1-c)’ 


M + F M + F + K M + F + 2K M + F+3K 
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It can just be written 

6. Mean weight in the catch 

7. Catch in weight 

8. Mean age in the catch 

9. Mean length in the catch 

10. Critical age 


B curo ' = R c -W_.[...] 


w _ J3cum =Z .W. .[...] 

Ncum 

Y = C W = F . B cum = F . Rc . W„ .[...] 


t = 


_t_ 

Z 


L L- <L " Lc ^Z + K 



Comments 

1. The simplification of Beverton and Holt allows the calculation of any characteristic of the 
cohort during its life with algebric expressions, avoiding the addition of the values of the 
characteristic in the successive intervals Tj. This was useful for calculations in the 60-70's 
when computers were not available. It is also useful when the only available data is 
natural mortality, M, and growth parameters. 

2. At present, the simplified expressions are also useful to study the effects on the biomass, 
on yield and on the mean weight of the catch due to changes in the fishing mortality 
coefficient, F, and/or on the age of first capture, tc. These analyses are usually illustrated 
with figures. For example. Figure 3.9 exemplifies the analyses of the biomasses and of 
the catches in weight, obtained from a cohort during the exploited life, subjected to 
different fishing mortality coefficients, assuming a fixed age tc . 



Figure 3.9 Evolution of the biomass and of the catch in weight from a cohort subject 
to different fishing mortality coefficients and fixed t t (notice that the Figure 
illustrates only the analysis and does not take into consideration the scales of 
the axis). 
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3. Notice that the forms of the previous curves, Y and B cum against F, do not depend on the 
value of the recruitment and so, they are usually designated as curves of biomass and 
yield per recruit, B/R and Y/R, respectively. The calculations are usually made with 
R=1000. 

4. The mean weight, the mean age and the mean length of the catch do not depend on the 
value of the recruitment. The curves of the characteristics of a cohort during its life 
against the fishing level, F, or against the age of first catch, tc, deserve a careful study, 
for reasons which will be stressed in the chapter concerning the long-term projections of 
the stock . 

5. B cum was calculated as: 

B« = Jnt W, dt 

1C 

where 

-(M»F|<l-lc> , I. -K(l-u»|) 

N, = Rc-e and W, = W-P — e J 

The calculations can also be made using other values of the constant b, from the relation 
W-L, different from the isomorphic coefficient, b=3, using the incomplete mathematical 
function Beta (Jones, 1957). 

6. The means L and t can be calculated from the cumulative expression 


L 


cum 



and 



as 


F Fcum/ Ncum and t — t cum / 

These means are designated as weighted means, being the weighting factors, the number of 
survivals, N, , at each age t. 
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CHAPTER 4 - STOCK 


4.1 STOCK OVER A ONE YEAR PERIOD 

4.1.1 EVOLUTION OF THE AGE STRUCTURE OF THE STOCK 

Let us consider the evolution of a stock, with several cohorts, over the period of one year. 
Figure 4.1 shows the structure of the stock at the beginning of the year, the mean 
characteristics during the year and the age structure surviving at the end of the year. 


ages 


Begin During End 



Structure of the stock Survivals at the end 

at the beginning of the year a of year a 

Figure 4.1 Structure of the stock at the beginning and end of the year 


Let / = 0,1, 2,3, ... be the ages 

4.1.2 CHARACTERISTICS OF THE STOCK AT THE BEGINNING OF 
THE YEAR 

Let : 


N; - number of individuals at age i, at the beginning of the year, 
w, - individual weight at age i, at the beginning of the year. 

B, - total biomass of the individuals at age i, at the beginning of the year. 
N - total number of survivors of the stock at the beginning of the year. 

B - total biomass of the stock at the beginning of the same year. 
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Then : 

N = IN, 

i 

B = IB, =l(N i .W,) 

i i 

4.1.3 CHARACTERISTICS OF THE STOCK DURING THE YEAR 

Let : 

M„ Fj e Z, total mortality (natural and by fishing) coefficients, 
at age i during the year 

W individual mean weight at age i during the year 

N mean number of individuals during the year 

B mean biomass during the year 

C catch, in number, during the year 

Y catch, in weight, during the year 

W LjKh mean weight of the individuals caught during the year 

W (tock mean weight of the individuals of the stock during the 
year 


Then: 


Ncum = X Nci 


T 

c=sc 

Bcum = £ Bcum i 


T 

Y = IY, 

W 

catch £ 


w = — 

N 


(where T=1 year) 


(where T=1 year) 


with Ncmj, Q, B cllm , and Y, calculated for all the ages i according to the expressions given 
previously in Chapter 3 - Cohort. 
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4.1.4 CHARACTERISTICS OF THE STOCK AT THE END OF THE 
YEAR 

The number of individuals at beginning of age i+i will be: 

Ni*i = Ni-e'Z' 


Let : 

R - recruitment of the cohort in the following year 

Then the number and the biomass of the stock at the beginning of the following year will be : 
N = R + lNi*i and B = R.Wi + Sb»i 

i i 

where the product R.Wi is the biomass of the recruitment of the following year. 

Comments 

1. The end of the year coincides with the beginning of the following year. So, the number of 
survivors of age i at the end of the year will be Nj+i, with age i+l. 

2. The sum of the total number of survivors of the stock at the end of the year is not equal 
to the number of individuals of the stock at the beginning of the following year, because 
one has to count the recruits entering that year. 

3. The total number of deaths, D. during the year, would be D = £ D, 

i 

4. As the interval of time is 1 year, the cumulative values will be equal to the mean values, 
that is : 

=B, 

N = Zn, 

i 

B = Sb, 

5. The utilization of the same symbols N, B, D, etc., for the stock and for the cohort should 
not create any confusions. 


40 


Copyrighted material 



4.2 


FISHING PATTERN OVER A ONE YEAR PERIOD 


4.2.1 FISHING LEVEL AND EXPLOITATION PATTERN 

The direct action of fishing a stock can be represented by the coefficients of mortality by 
fishing Fj. These coefficients are associated with the quantity of effort, with the disponibility 
of the individuals of different sizes or ages, i, and with the fishing gears used by the vessels 
during the year. 

It is usual to separate the coefficients of mortality by fishing into two components: 

One is designated as the level of intensity of mortality by fishing, F , during the year, called 
fishing level, F. The level is associated with the quantity of fishing effort, (the number of 
vessels fishing), the number of days, hauls, fishing hours during the year, and with the 
efficiency or the fishing power of the vessels or gears. Another component, designated as 
exploitation pattern, Sj, is associated with the selective properties of the fishing gears relative 
to the sizes or ages of the individuals available to be captured, during that year. 

The combined set of the fishing level (a unique value for all ages) and the exploitation pattern 
(different values according to the size or the age), is designated as fishing pattern or fishing 
regime. The designation of fishing pattern may cause some confusion with the exploitation 
pattern, and the designation of fishing regime may be confused with what the economists and 
managers call fishing regime. 

The fishing pattern, Fj, of one age i during one year, is equal to the product of the fishing 
level of that year, F , times the exploitation pattern of each age, Sj. That is: 

F, = F. S j 

To analyse the effects of the coefficient of mortality by fishing, F , on the characteristics of 
the resource and on the catches, the exploitation pattern is generally taken as constant from 
one year to the next. Sometimes one analyses the effects of the changes of the exploitation 
pattern maintaining a fixed fishing level, but one can also analyse the combined effects of the 
two components. 

The fishing level, F , is often represented by F. 
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4.3 SHORT-TERM PROJECTIONS OF THE STOCK 

Knowing the structure of the stock at the beginning of one year, it is possible to estimate the 
characteristics of the stock during that year and project the structure of the stock for the 
beginning of the next year (with an exception to the recruitment of that year), for different 
values of the fishing level, F , (and for the exploitation pattern, s,). 

It is, of course, necessary to know the biological parameters of growth, maturity and the 
natural mortality coefficients, in each age during the year. 

Adopting one value for the recruitment of the following year, the projection could be repeated 
for one more year and so on. The inconvenience in projecting the stock for several years will 
be that those projections will depend more and more on the adopted annual recruitment 
values. That is why, in practice, the stock and the catches are projected for one, or at the most, 
two years. In Section 4.5, the estimation of recruitments will be discussed. 

Notice that to project a stock for the following year, it is necessary to have data from the 
previous year. So, the stock is firstly projected for the current year and then the catch and the 
biomass are projected for the following year. 

Let us suppose, for example, that in 1997 one wants to project the characteristics of the stock 
for 1998. As the available data will be, in the best hypothesis, referring to 1996, one will have 
to project the stock of 1996 for the beginning of 1997 and together with the recruitment of 
1997, project the stock until the end of 1997 and only after that, can one project to 1998, 
estimating previously the recruitment of 1998. 

4.4 LONG-TERM PROJECTIONS OF THE STOCK 

Let us consider the vector N with the components N|, N 2 , ..., Nj, ... representing the structure 
of the stock at the beginning of a year. Notice that Ni, N 2 , .... Ni, ... are the survivors of the 
different cohorts of the stock, at the beginning of the year. 

Let M, and F, be the natural and fishing mortality coefficients of age i, assumed to be constant 
in the future years. 

Figure 4.2 illustrates, with a theoretical stock, the projections of the numbers of survivors of 
the different cohorts at the beginning of 1980, for the years to come, from 1981-1986 (the 
values are expressed in million of individuals). 

Notice that the recruitments are missing during the years 1981 to 1986, as they have not yet 
occurred. So, it is clear that the respective survivors are also missing during those years. 

Let us suppose now that the recruitments of the same period of years were equal to those of 
1980, that is, 440 million of individuals. Figure 4.3 shows the projections in future years. It 
can be seen that the values of the age structure of the stock in 1986 are equal to the annual 
evolution of the cohort in 1980. 
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Year 

1980 1981 1982 1983 1984 1985 1986 


Age 

0 

1 

2 

3 

4 

5 

6 



Figure 4.2 Structure of the stock at the beginning of 1980 and projections of the 
number of survivors (million of individuals) of the different cohorts 
at the beginning of the years 1981-1986 


Year 


Age 

1980 

1981 

1982 

1983 

1984 

1985 

1986 

1987 

0 

|440 

440 

440 

440 

440 

440 

440 

440 

1 

995 

326 

326 

326 

326 

326 

326 

326 

2 

367 

680 

223 

223 

223 

223 

223 

223 

3 

68 

229 

425 

139 

139 

139 

139 

139 

4 

245 

41 

139 

258 

85 

85 

85 

85 

5 

345 

149 

25 

84 

156 

51 

51 

51 

6 

76 

209 

90 

15 

51 

95 

31 

31 


Figure 4.3 Long-term projections of the cohort of 1980 and structure of the stock 
assuming a constant annual recruitment of 440 million of individuals 

One practical conclusion seems to be that to obtain the structure of the stock in 1986 it would 
be enough to follow the evolution of the cohort of 1980 and then, it would not be necessary to 
have the complete structure of the stock at the beginning of the year. It would be enough to 
adopt one value for the recruitment ( R ) of a cohort (and, of course, assume that the 
biological parameters would be stable in the following years). An advice would be to make 
the calculations adopting a Recruitment of 1000 individuals (with computer software a 
recruitment equal to 1 is adopted). 

Except for the mean values - age, length, weight in the catch, etc, - which are independent of 
the adopted value of R, the characteristics of the stock in the long-term, under the previous 
conditions, are proportional to the recruitment. So, these projections are also designated as 
projections by recruit , by LT projections or even as equilibrium projections. 

The long-term projections do not only concern the survivors at the beginning of the year, but 
also the values of the mean abundances, N, during the year, and the catches in number, C, . 


43 


Copyrighted material 



One can also project the total biomasses, B or B, the spawning biomasses, SB or SB , as 
well as the catches in weight, Y, knowing the initial individual weights, Wj, the mean 
weights, W , and other parameters like those of maturity. 

It is important to verify that the cumulative values, Ncum; e Bcum,, of a cohort during 1 year 
are equal to the mean values, Nj e B, , during that year. The long-term projections can be 
calculated as: 


N S iock = Z Ncumj = Z Nj 

Bs,„ck = Z Bcum: = Z B ( ( SB stock =Z SB,) 
Wsiock = Bgtock / Naodt 
C = ZQ 
Y = ZYi 
W MKh = Y/C 


Several long-term projections can be made with different values of Fi, that is, with several 
fishing levels, F , and/or with several exploitation patterns, s,. 

As mentioned before, the analyses of the effects of the fishing pattern on the catches and 
stocks can be done with the two components (fishing level and exploitable pattern), separated 
or combined. 

Figure 4.4 (A-C) illustrates several types of yield per recruit curves, maintaining a fixed 
exploitation pattern. The curves are different for other exploitation patterns. 



Figure 4.4 Examples of types of curves of the Yield per Recruit (Y/R) against F, 
given the exploitation pattern: (A - with a maximum, B - flat-top, 

C - asymptotic) 
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Figure 4.5 (A-E) illustrates the relations between the most important characteristics of the 
stock and fishing level, maintaining a stable exploitation pattern: 



Figure 4.5 Relations between the long-term characteristics of the stock against the 
fishing level F, (A - mean number/R, B - catch in number/R, C - mean 
weight in the catch, D - mean biomass /R, E - yield/R) 

A more detailed analyses of these curves will be presented later on, in Chapter 5. 


Conclusion 


The long-term structure of the 



Evolution of the cohort 

stock by ages during 1 year 


during its life 


Comments 

1. One projects a cohort during its life in order to obtain the long-term projection of the 
stock for one year, assuming the annual recruitments to be constant. 

2. It is necessary to know M; for all ages of the cohort, as well as Wj , w , sj 
and the fishing level, F, that is assumed to be constant in the years that follow. 

3. Any recruitment size can be used. Adopt R = 1000 (or R = 1) with worksheets on 
your computer. 

4. The five most important characteristics of the stock arc N SIO di, B sloC k, C, Y and 
w ouch (see previous pages in this chapter for the respective expressions of 
calculation) 

5. A characteristic of the stock that is also important is the spawning biomass, SB. To 
calculate SB, it is necessary to know the maturity ogive (or histogram). 

6. Long-term projections are also designated as equilibrium situations. 

7. Long-term projections are useful to define the long-term management objectives. 
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Annual W gaK h , are independent of the recruitment size (such as Loich and t c „ch ). 

9. Economists transform the total yield, Y. into value YS, the mean weight of the catch, 
W cllth into price of the catches. ws cali:h , the catch per vessel (or the cpue ) into value 
of the production by vessel, US, and the fishing level, F into costs of exploitation, 

F$. The difference between the value of the catch and the cost of exploitation, YS-FS. 
is the profit, L$. Figure 4.6 illustrates an example of the LT relations between those 
characteristics used by the economists against the fishing level, F. 



Figure 4.6 Long-term relations between the value of the catch, the cost of exploit- 
ation and the profit against the fishing level 


4.5 STOCK-RECRUITMENT (S-R) RELATION 

The stock-recruitment relation, known by S-R relation, associates the size of the stock, for 
one year, with the recruitment which results from the stock spawning normally during that 
year. The recruitment can be the recruitment at the exploitable phase and the stock can be the 
spawning stock. This recruitment may occur one or more years after the spawning. 

The problem with the relation between the parental population and the new generation is not 
a special case of the fisheries resources, it is common to all the self renewable populations. 

It is important to determine, in each case, the stock and the recruitment to be used. In fact, 
that stock can be the total number of individuals (at the beginning of the year or the mean 
value during the year), the number of adult individuals of the stock, the number of adult 
females, etc. One can also adopt, not the abundances in number, but the corresponding 
biomasses. The decision will depend on the type of resource and on the available data. It is 
necessary to define if the recruitment is in weight or in number and if the recruitment is to the 
fishing area or to the exploited phase. 

In this manual, stock (S) will be considered as the spawning biomass and the recruitment (R) 
will be expressed in number. 

After the spawning, the individuals of the new generation will have to go through different 
phases of the biological cycle: eggs, larvae and juveniles, before they become recruits. These 
phases, which, in some species can take years, are not directly submitted to the fisheries 
exploitation. That is why fishing in those years does not directly affect the size of the new 
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recruitment. It is true that fishing acts on the size of the parental biomass, but that does not 
happen in the pre-recruit phases, and it is precisely for that interval of time that the relation 
S-R is applicable. 

There is, in these pre-recruit phases, a great mortality due to climate and environmental 
factors (winds, currents, temperatures, etc.) as well as due to biological factors (available 
food, predation and others). 

A great variety of factors (besides fishing) cause great fluctuations to the recruitment sizes, 
therefore, the relation S-R is a complex one. In conclusion, the relation S-R, between the 
stock and the resulting recruitment, can be considered independent of fishing. 


R 



“L 

0 


s 


Figure 4.7 Example of the dispersion of the points in the relation S-R 

Figure 4.7 shows the type of dispersion of values in the plot of recruitment (R) against 
parent stock (S). 

Despite the difficulties, some models have been proposed for the relation S-R. One of the 
reasons for this is that there must be a relation. One point of the relation (S = 0, R = 0) is 
even known already. 


4.5.1 BEVERTON AND HOLT MODEL (1957) 

Beverton and Holt, when analysing plaice fishery in British waters proposed one model, 
which can be re-written as: 



where a and k are constants. 
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s 


Figure 4.8 Beverton and Holt Relation S-R 
This model shows a curve where R is asymptotic when S— > 

4.5.2 RICKER MODEL (1954) 

Ricker, studying cod fishery in Canadian waters proposed another model that can be written 
as: 

R = a.S.e' 5 * 

where a and k are constants. 

R 
0 

0 s 



Figure 4.9 Ricker Relation S-R 

This model presents a maximum resulting recruitment at an intermediate value of the parental 
stock, S. 

4.5.3 OTHER MODELS 

Other models of S-R have been proposed, like the Deriso generalized model (1980), which 
can be re-written (Hilbom & Walters, 1992) as : 

R = aS.(l-c.-) 1/c 

k 

(the model is the Beverton and Holt for c = -1 and when c— » 0 it is the Ricker model) 
or Shepherd’s generalized model (1982), as: 
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(the model is the Beverton and Holt when c = 1 and when c = 2 the curve is an pproximation 
of the Ricker model) 

Notice that in every model presented, S=0 implies R=0, as expected and the slope of the 
tangent to the curve at that point, (0,0), is equal to the parameter a. Sometimes, the relations 
S-R are presented as R/S in function of S, so the Beverton and Holt equation would be : 


R, 

S 


— showing that the inverse of R/S is linear with S. 

1 + - 

k 


With the Ricker model the relation between R/S against S would be: 



that is, R/S is exponential negative with S (or ln(R/S) is linear with S). 

With the Deriso model it would be (R/S) c linear with S. 

Mathematically, these linear relations can be useful to estimate the parameters a and k, but 

statistically they cause some problems, because the variable S appears in the response 

variable and in the auxiliary variable. 

Comments 

1 . Remember that relations S-R are independent of the fishing level. 

2. Relations S-R may be introduced in the calculations of the stock projections. In that 
case, the projections will have to be made every year and they will require the structure 
of the stock at the beginning of the initial year. 

3. It must not be forgotten that there is a great dispersion of the observed points (S, R) 
around the model. 

4. To determine the limit of the Deriso equation when c — » 0 it is enough to remember that 
limit (l+A/n)“ = e A when n— > » and A is a constant. 

5. The Beverton and Holt model presents an asymptote, R —> ot.k when S — *=°. When S=k it 
will be R = a k/2. 


The Deriso model presents the maximum: 

Sm»= k / (1+c), Rmax = a.k / (l+e) |4<,,c) 
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The Shepherd model presents the maximum: 

S M x= k.(c-l)' /c , Rmax = (a/c).k.(e-l) M,/c) 

4.6 RELATION BETWEEN R AND B (R-S RELATION) 

Up to now the discussion have been centred on the relation S-R, that is the relation between 
the biomass S and the resulting recruitment R. There is another relation (which has already 
been referred to in Section 4.4 Long-term projections of the stock, particularly in the 
conclusion about the structure of the stock and the evolution of a cohort during all its life) that 
could be called the relation R-S, that is the relation between the recruitment R and the 
resulting cumulative biomass, B cum of a cohort during all its life for a given fishing pattern. 

The cumulative biomass (or spawning biomass) of a cohort during its life, B cun , , is, as it has 
already been seen, equal to: 

Bcum— R • {function of biological parameters and of the fishing pattern, F} 

As mentioned before, the cumulative biomass, B cum , of a cohort during its life is equal to the 
long-term mean biomass, B , of the stock during one year, so, one can refer to the B cum as B. 

it can then be said that in the long-term, if the biological parameters are considered constant, 
then for a certain fishing pattern, the mean biomass of a stock during one year is proportional 
to the recruitment. Figure 4.10 illustrates the proportionality. 



Figure 4.10 Illustration of the proportionality in the long-term between the mean 
Biomass B of a stock and the Recruitment (R), for a certain fishing level, 
F, assuming that the exploitation pattern and the biological parameters 
are constant 

Notice that in Figure 4.5-D, the relation B/R against F was shown. The two curves 4.5-D and 
that of Figure 4.10 are different representations of the same situation. While in Figure 4.5-D, 
a value of F corresponded to a value of B/R given by the curve, in Figure 4. 10 the value of F 
will correspond to a slope B/R of a given straight line. 

For other values of F, the straight line of Figure 4.10 will have different inclinations. 
Figure 4. 1 1 shows several lines according to the different values of F. 
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Figure 4.1 1 Illustration of four lines corresponding to different fishing levels, F 

Figure 4.12 shows the overlapping of the Ricker curve S-R (Figure 4.9) and the line in 
Figure 4.10, the axes of this last figure having been switched. To avoid confusions, it is 
convenient to refer the slope of the line corresponding to a value F in relation to the axis R. 

Long-term 



Figure 4.12 Overlapping of the curve S-R with the line R-S for a given F 

One can start from a value of biomass, and, through the relation S-R, determine the future 
recruitment, R. That R will give a resulting biomass, B , through the straight line. This process 
can be repeated until a situation of equilibrium is found. 

It will be possible to illustrate the combined analysis through the two relations for a certain 
fishing level. 

For example, let us select one value B|. The curve S-R allows the calculation of the resulting 
value Rj. Through the line (of a given F). the resulting value of Bj will correspond to that 
value of R;, and so on. Figure 4.13 shows that the process reaches the equilibrium point 
(Re, Be). This point would correspond to a value of the fishing level F, determined by the 
slope of the line in relation to the R axis. 
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LT 



Figure 4.13 Illustration of the process that, starting from a biomass Bi, theoretically 
leads to the equilibrium point (R E , B F ) 


Notice that the intersection of the two relations does not always lead to an equilibrium point. 
The existence of an equilibrium point depends on the angle between the straight line and the 
curve at the intersection point. 
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CHAPTER 5 - BIOLOGICAL REFERENCE POINTS 
AND REGULATION MEASURES 


5.1 BIOLOGICAL REFERENCE POINTS FOR THE MANAGEMENT 
AND CONSERVATION OF FISHERIES RESOURCES 

The long-term objectives for fisheries management should take into consideration scientific 
fishing research and population dynamics, as well as the climatic changes that may affect the 
stocks. 

In order to define these long-term objectives we have to consider the values of the fishing 
level, which allow bigger catches in weight, whilst also ensuring the conservation of the 
stocks. The extreme values of the biomass or the fishing level, which might seriously affect 
the self renovation of the stocks, also have to be considered. These fishing level values, of 
catch and biomass are designated as biological reference points (BRP). In this manual some 
of the different types of BRP will be considered (Caddy, & Mahon, 1995; FAO, 1996 and 
ICES, 1998). 

The Target Reference Points, TRP are BRP defined as the level of fishing mortality or of the 
biomass, which permit a long-term sustainable exploitation of the stocks, with the best 
possible catch. For this reason, these points are also designated as Reference Points for 
Management. They can be characterized as the fishing level Fo^, (or by the Biomass , B urgcI ). 

The most well known F UIgc , is Fo i but other values, like F max , F^d, and Fmsy will also be 
studied. 

For practical purposes of management, the TRP will be converted, directly or indirectly, into 
values of fishing effort, given as percentages of those verified in recent years. 

The Limit Reference Points, LRP are maximum values of fishing mortality or minimum 
values of the biomass, which must not be exceeded. Otherwise, it is considered that it might 
endanger the capacity of self-renewal of the stock. 

In the cases where fishing is already too intensive, the LRP can be important to correct the 
situation or to prevent it from getting worse. 

The LRP are limit values, mainly concerned with the conservation of marine stocks and they 
are therefore also referred to as reference points for conservation (this designation does not 
imply that the Fargo are not concerned with conservation). 

Several LRP have been suggested, which will generally be referred to as Fi, m or Bii m . In this 
manual the levels of biomass Bi<« s and MBAL will be referred to as Bi, ra and the fishing levels 
Floss and F ,ra,h as F i lm ■ 

The Precautionary Principle, proposed by FAO in the Conduct Code for Responsible 
Fisheries (FAO, 1995), declares that the limitations, uncertainties or lack of data for the 
assessment or for the estimation of parameters, cannot be justification for not applying 
regulation measures, especially when there is information that the stocks are over-exploited. 
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From this point of view, it is important to make clear which basic assumptions are necessary 
in order to estimate the consequences on the catches and on the abundance of the stocks. 

The uncertainties associated with the estimation of Fi, m , and B| lm , therefore lead us to 
determine new reference points, called Precautionary Reference Points, Fpa or Bpa. 

The assumptions and the consequences of adopting alternative hypotheses about the stock and 
fishing characteristics should always be presented to justify the estimated values of Fpa (or 
Bpa). 

The new limits (Fpa or Bpa) due to the application of the Precautionary Principle, will be 
more restrictive than the LRP's. The practical consequences of these new limits are the 
regulation measures designed to control the fishing effort which are more severe than in those 
cases where there is appropriate data. 

It can be said that this is the price to pay for not having the appropriate conditions to make 
available reliable data and information. 

The Precautionary Approach, suggests that the results of fisheries research should be adopted 
by management with regard to the formulation of the regulation measures and that these 
measures should also take into consideration the socio-economic and technical conditions of 
fishing (FAO, 1996). 

A final remark about all the Biological Reference Points mentioned above : 

The evaluation of the biological reference points has to be updated, taking into consideration 
the possible changes in the biological parameters or any other necessary correction of the 
exploitation pattern. This fact is important because the new biological reference points will be 
different from the previous ones. 
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5.2 BIOLOGICAL TARGET REFERENCE POINTS 

(^man I" 0.1 Fmed Ulld Emsy) 

5.2.1 F mas 
Definition 

1 Consider the long-term yield per recruit, Y/R, as a function of F, for a certain 
exploitation pattern. 

F ma . x is the point of the curve, Y/R against F, where Y/R is maximum. 

Figure 5.1 shows a curve of Y/R against F. 

MAX 



Figure 5.1 Y/R as function of F for a certain t c constant, 
showing F m „ and Y m „ 

In Chapter 4 it was mentioned that to estimate the long-term projections one 
could assume that the recruitment is constant and equal to 1 (R=l). In this way, 
the mathematical expressions are sometimes w ritten with Y instead of Y/R. 

2. Mathematically, at point F nu „, the derivative of Y/R against F is equal to zero, that 
is. 

For F = Fma,willbe ^777 = 0 
dF 


For F < Fma X will be — > 0 
dF 

For F > F m „ will be 7777 < 0 
dr 

Geometrically, the slope of the tangent to the curve is equal to zero for F = F ra „ , 
positive for F < F ra „ and negative for F > F nM1 . 


air(Y ) = 0 (Value of Y is maximum) 
air(Y ) > 0 (Y is increasing with F ) 
air(Y ) < 0 (Y is decreasing with F) 
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Comments 

1 . B max /R and Yn»j/R are the values at F nmx . 

It is also convenient to analyse the situation of B/R at the points F# F m „. 

F < F lna , corresponds to B > 

F > F„„ corresponds to B < B m „ 

Point Fmn does not depend on the value of the recruitment. 

2. For another relative pattern of exploitation there will be another F„ u . 

3. All the points of the curve Y/R against F, are long-term points or equilibrium 
points. 

4. When the level, F, is bigger than F^,* it is said that there is growth overfishing. 

It is convenient to present the two curves Y/R and B/R against F, in the same 
graph (usually with different scales). 


F 

Figure 5.2 Long-term curves of Y/R and B/R against F, given an 

exploitation pattern 

Fmax was adopted by the majority of the International Fisheries Commissions as a 
long-term objective of management (1950-1970). 

Even today F^ is used as a target-point having been proposed as a Limit 
Reference Point (LRP) in some cases. 

The flat-top and asymptotical curves do not allow the determination of an F mlx . 
The definition of F TO> does not consider the appropriate level of spawning 
biomass. 

only indicates the value of F which gives the maximum possible yield per 
recruit from a cohort during its life, for a given exploitation pattern. 

The analyses of these long-term curves, mainly of B, Y and W CM£ll against the 
fishing level, give information about the abundance of the resource (or catch per 
vessel), total yield of all the fleet and mean catch weight for different fishing 
levels. 
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5.2.2 F 0 ., 


1- Definition 

Consider the long-term yield per recruit, Y/R, as a function of the coefficient of fishing 
mortality, F. One value of air(Y/R), corresponds to each fishing level, F. The air(Y/R) is 
maximum when F = 0 and decreases, being zero when F = F„ m . 

The point F ol is the value of F where air(Y/R) is equal to 10 percent of air(Y/R) 


maximum. 


The figure 5.3 illustrates this situation. 


air(yM0%.air(y) 


air(y) maximum 


Figure 5.3 


Y/R showing the reference target point Fo.i 


2- For F = 0, the biomass per recruit, B/R will be B 0 /R, also designated as Virgin 
Biomass or Non-exploited Biomass. The airfy) at F=0 is also equal to Bo 


Y = F B implies — = B + F— 
dF dF 


Then, for 


>=M =1 

Uf l 


So, from the definition given in point 1 one can also say that Fo i is the value of F where 
air(Y) = 10 percent of the virgin biomass. 

3- Calculation of Fo.t 

Let the function V = Y -0.1.B„.F 

It can be proved that the function V is maximum when F = F„ , 
dV 

In fact, V is maximum when — = 0 . then: 
dF 


dV _ dY 
dF dF 


-0.1 B„=0 


Therefore, the value of F corresponding to the previous dY/dF is the value of Fo i. 
Fo.i can then be calculated by maximizing the function V = Y - 0. 1 B~.F 
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Bo can be calculated, for example, from the long-term relation of B against F, 
when F = 0. 

Graphically it will be: 

long-term 



Figure 5.4 Curve Y/R showing the maximum of the function V 

4. Why adopt air(Y/R) equal to 10 percent and not any other percentual value, for example, 
20 percent? 

Gulland and Boerema (1969) presented some arguments, including financial arguments. 
Some countries (like South Africa) adopt the value of 20 percent with a resulting 
reference point F 0 j that is more restrictive than F 0 |. 

5. Figure 5.5 illustrates the two biological reference points F™, and F 0 ( . 



Figure 5.5 Long-term variation of Y/R and B/R against F 
and points corresponding to F m „ and F 0 .i 

B 0 , and Y 0 , arc the values of B and Y corresponding to F 0 1 

F 0 I is always smaller than F nu , 

B 0 1 is always larger than B 

Y„ , is always smaller than Y^,, although, in practice, the difference is not large. 

The second sentence above indicates the advantages of B„ , over B raa ,. The last sentence 
shows that Y 0 , is not the largest possible catch, but is acceptable as a target point of 
management. The fact that B ol is larger than B suggests that the fishing level F 0 , is 
preferable to Fma X as TRP. 
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Notice that Fo.i can be calculated even when the curve is asymptotical or flat-top. 

Another value of Fo i will be obtained if the exploitation pattern changes. 

In the years 1960-70 Fo i started to be preferred to F, a „ as a target point by resource 
managers, and it was adopted in the 80’s, as a long-term objective by many International 
Fisheries Commissions and by the EEC. 

5.2.3 F raed 
1- Definition 

This target point considers the relation S-R between the stock and the resulting 
recruitment, in order to avoid the assumption of a constant recruitment. 

Let (the spawning or total) biomasses and the resulting recruitments for each year of a 
certain period of time be known. In this case, one can calculate the median value of the 
ratios between the annual spawning biomasses and the corresponding recruitments. 

Fined is defined as the F value corresponding to the median B/R ratio in the long-term 
B/R relation against F. 

Usually, Fmd is illustrated by considering the graph of the points corresponding to pairs 
of values of parental biomass (total or spawning), B, during that year and the resulting 
recruitment , R. Figure 5.6 shows this situation. 



Figure 5.6 Illustration of a median line 

The marked line is a line passing through the origin, which separates the total number of 
points in equal parts, that is, 50 percent of the points are in the upper part and 50 percent 
are in the lower part of the line. This line is designated as the median line, or 50 percent 
line, which can be explained as follows: in 50 percent of the years of the considered 
period the values of R were smaller than the values of R which were estimated by the 
median line (or, in 50 percent of the years of the referred period the values of R were 
bigger than the values of R estimated by the median line). 

As seen in Section 4.5 the slope (B/R) of each line marked in the graph, is associated 
with a certain value of the fishing level, F. The value of F associated with the median line 
is then, the median target point, F mc j. 
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It can be said that, given a certain level of parental biomass and knowing the B/R 
corresponding to F ine<] , then there is a 50 percent probability that the resulting recruitment 
will be less than (or greater than) the value indicated by the median line. 


Calculation of F mt<i 

In order to determine the value of Fmcd it is necessary to consider the long-term relation 
between the resulting biomass per recruit and the fishing level, F, (Section 4.4, Figure 4- 
D). 

The determination of Fmcd can be done mathematically or graphically. 

To make the mathematical calculation of Fmcd the ratio B/R has to be determined for each 
pair of values (B, R). Those values have to be ordered and then the median 
value, B/'Rmcdian can be calculated. 

In the long-term relation of the Biomass per Recruit against F, the value of Fmcd is the 
value of F that corresponds to the median value previously obtained. 

To make the graphical calculation, notice that in Figure 5.6 the slope of the median line 
against axis R is equal to B/R. This value will be the basis for the calculation of value 
Fmcd in the graph of B/R against F, in the long-term projection. Figure 5.7(A-B) 
illustrates the calculation. 

A) B) 




Figure 5.7 Illustration of the graphical calculation of Fmcd 

Other straight lines, corresponding to probabilities other than 50 percent, can be marked. 
Figure 5.8 illustrates the graphical calculation of Fio%. The line marked on Figure 5.8A, 
separates the points, in such a way that 10 percent of them remain below the line (or 
90 percent of the points stay above the line). So, this line has been designated as the 
10 percent line. 

Notice that the slope of the 10 percent line with axis R is larger than the slope of the 
median line, as shown in Figure 5.9B. In this way, Fio»/„ or F| ow , is smaller than F,„ cd . 
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A) 


B) 




Figure 5.8 Illustration of the graphical calculation of Fio% 


A) 


B) 


R 


median lines 


H.KI 




Figure 5.9 Graphical illustration of the slopes and corresponding F's 
for the median and 10 percent lines 


Comments 

1 . The target point F^ intends to ensure an acceptable level of biomass based on the 
empirical relation S-R. 

2. Other percentages can also be adopted, corresponding to straight lines which estimate 
different probabilities of recruitments which are less than those indicated by the median 
line. So, Fhigh would be the fishing level corresponding to the 90 percent line, for which 
the recruitments of 90 percent of the observed years would be less than those estimated 
by the line. 

3. F 9 o%, as will be seen in the following Section, is also considered as a Limit Point (LRP). 

4. Notice that the slope of the 10 percent line with axis R is larger than the slope of the 
median line and therefore, Fio% (or F| 0W ) is smaller than Fmcd (Figure 5.9B). 

5. Fmd was used in management, in recent years, particularly with iberic sardines. 

6. The biomass used can be the total biomass, B, but is frequently the spawning biomass, 
SP. 

7. If the median line does not pass through a marked point in the scatter plot (which 
happens when there is an even number of points) then one can use any straight line 
passing between the two central points, for instance, the mid-line. In any case, Fmd is 
always an approximate value. 
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5.2.4 F msv 

Definition 

F MS v is defined as being the value of F which produces the maximum yield in the long-term. It 
is necessary to select an S-R relation to estimate Fmsy- This point is different from F rn „. 

5.3 BIOLOGICAL LIMIT REFERENCE POINTS 
(B, OSS’ MBAL, F crash snd F| WS ) 

There are several proposals for F| lm and B| lm . For each stock, the adopted values of F [jm and 
Bi,m depend on the characteristics of the stock and on its exploitation. What is important is 
that the adopted LRP be a value that allows an exploitation level which avoids dangerous 
situations of stock renewal. 

Some of these points are derived from the observed values of Biomass and of Resulting 
Recruitment. Some examples of this type are B| oss and MBAL. These LRP are also usually 
classified by some authors as non-parametric, because their determination does not depend on 
any particular model of the S-R relation. 

Another category of LRP points, classified as parametric, is derived from S-R models. 

Foash will be mentioned. 

Let us also mention the category of LRP points involving observed values and values 
obtained by the application of S-R models, like, for example, F| 0SS . 

5.3.1 B lo „ 

Bios 5 is the smallest spawning biomass observed in the series of annual values of the spawning 
biomass (Lowest Observed Spawning Stock). 

5.3.2 MBAL 

More satisfactory is the LRP designated as Minimum Biological Acceptable Level, MBAL. 
In fact, this LRP is a spawning biomass level below which, observed spawning biomasses 
over a period of years, are considered unsatisfactory and the associated recruitments are 
smaller than the mean or median recruitment. 

5.3.3 F cras h 

The name itself shows that it is a limit that corresponds to a very high value of F, showing a 
great probability of collapse of the fishery. 

Fcrash is the fishing level F which will produce a long-term spawning biomass per recruit 
(S/R) equal to the inverse of the instantaneous rate of variation of R with the biomass, at the 
initial point (S = 0, R = 0). With the S-R models of Section 4.5 that value is the parameter 1/a 
of the models. 

In order to make the graphical determination of this LRP one can start by obtaining the slope 
of the angle that the tangent to the S-R curve makes with the R axis at the origin. Afterwards, 
and starting from the relation B/R against F, the value of F that corresponds to the value 
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B/R indicated by that slope is obtained. F craS h will then be the value of F corresponding to 
B/R equal to that slope, in the long-term relation B/R against F. 


F| 0SS is usually defined as the fishing level F which will produce a long-term spawning 
biomass per recruit (S/R) associated to B| 0i , . 

To determine this limit point, first obtain the value of R corresponding to Bi oss on the adjusted 
curve S-R. Then, calculate Bi 0 »/R and find the value of F, in the long-term relation B/R 
against F. 

Most of the Limit Points shown have been criticised for depending on the observed values or 
on the adjustment of the S-R relation. 

5.4 PRECAUTIONARY REFERENCE POINTS - Fpa, Bpa 

As previously mentioned, the Precautionary Principle recommends that the assessments 
should be done even when the basic data presents some gaps. This recommendation implies 
that, in this case, the determination of the Biological Reference Points will not be very 
precise. The uncertainities of the estimates should be calculated, and it is necessary to 
mention the assumptions and models which have been used. 

One suggestion to determine Fpa and Bpa might be to estimate F| lm or Bi, m and from these 
values, to apply the following empirical rules : 

Fpa = F| im .e' 16450 and Bpa = B, im .e* 1 6450 

The constant a is one measure associated with the uncertainity in the estimation of the fishing 
mortality level, F. The values obtained in several fisheries indicate that values of a are within 
the interval (0.2, 0.3) (ICES, 1997). In practice, it can be said that Fpa is between 0.47F|j m and 
0.61F Um ,and Bpa is between 1.39Bu m and 1 ,64Bi, m . 

It is important to make clear that the target points may also, in certain cases, be considered as 
limit or precautionary points depending on the combined analyses of the exploitation of the 
stock and of the biological reference points obtained. 
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5.5 FISHERIES REGULATION MEASURES 

The regulation measures aim to control the fishing level and the exploitation pattern applied 
to the stock for an adequate exploitation. 

The most common regulation measures to control fishing levels are: 

• Limitation of the number of fishing licences. 

• Limitation of the total fishing effort each year (limiting fishing days, number of trips, 
number of days at sea, etc.). 

• Limitation of Total Allowable Catch (TAC) 

TAC is a measure that directly controls the catch and, indirectly, the fishing level. It is 
convenient to combine the TAC with the allocation of quotas of this total TAC for each 
component of the fleet. In this way, the competition between vessels to fish the maximum 
possible catch, as quickly as possible, before the TAC is reached, can be avoided. 

The system of quotas allocated to each vessel is called Individual Quotas (IQ). 

The regulation measures to correct the exploitation pattern are usually called technical 
measures. Some of these measures are : 

• The minimum size (or weight) of the landed individuals. 

• The minimum mesh size of the fishing nets. 

• The prohibition of fishing in spawning. 

• The closed areas and periods for the protection of juveniles. 

The fishing management have the duty to promote legislation and the application of the 
regulation measures. (In the particular case of the EU, and for the stocks of the Economic 
Exclusive Zones (EEZ's) of the member states, the Commission decides on the measures to 
be taken). In any case, management needs the analyses on the state of the stock and its 
exploitation and on the effects of the recommended measures. That study must be done by the 
fisheries scientists of each country or region and their Fishing Research Institutes, who will 
have to calculate the projections of the stocks. The International Council for the Exploitation 
of the Sea (ICES) analyses the assessments and recommends regulation measures and the 
expected effects of the application of those measures to the Commission, as well as the 
consequences of their non-application. 

The short-term projections, as well as the regulation measures, only make sense if the 
long-term objectives of fisheries management are previously analysed and defined. 
Short-term projections of the stock and of its fishing must also be made by the scientists. 

Comments 

1. Management needs to define the fishing objectives, based on the long-term projections. 
Those objectives are valid for a period of years, even if they can be adjusted every year. 
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2. The regulation measures, on the contrary, have to be established every year, although 
some of them may be valid for more than one year. Some technical measures, like the 
minimum mesh sizes of the fishing nets or the minimum size of the landed individuals, 
are valid for several years. 

3. All the measures have advantages, difficulties and disadvantages regarding the purposes 
they intend to reach. 

The concession of fishing licences is a common practice almost everywhere, with a 
limited total number. 

TACs and quotas, because they control the catches, have caused misleading declarations 
about catches. 

The direct limitation of the total fishing effort (f), is based on the assumption that the 
measure causes a similar limitation on the fishing mortality coefficient (F). However, this 
relation may not be proportional. In the first place, it is difficult to measure the fishing 
effort of the different fishing gears and of all the involved fleets and it is also difficult to 
express it in units that respect the proportionality between F and f. Secondly, the 
capturability of several gears may increase (and consequently increase F) without 
increasing the fishing effort. Finally, the expected proportionality between F and f may 
not be true. In any case, what matters, is not to forget that there is a relation between F 
and f. 

4. The protection of the juveniles should be carried out during the whole year and will 
preferably control the fishing mortality throughout the year. The occasional measures, 
like areas and periods when fishing is forbidden for the protection of juveniles, require 
annual investigations in order to discover whether there are exclusive concentrations of 
juveniles, to assess the effects of that occasional interdiction, and to find out the 
consequences of the interdiction on other species, etc. The minimum size of the landed 
individuals, does not mean that smaller individuals are not caught, but only that they are 
not landed. The difference between the catch and the landing is the so-called rejections to 
the sea. It is clear that if the individuals are caught and rejected to the sea, the fishing 
mortality is larger than the one suggested by the landings. The minimum landing size of 
the individuals may have the effect of dissuading fishermen from catching small 
individuals. Currently, some countries are forcing the landing of all fish caught. 

5. The closed spawning areas and seasons, to save the spawning biomass and indirectly 
protect recruitment, is far from effective in the latter objective. In fact, large spawning 
biomasses correspond to a large number of eggs, but that does not necessarily imply 
bigger recruitments, as seen in Section 4.5. It is also not always true that forbidding 
fishing during the spawning, and not forbidding before (or after) the spawning, protects 
the spawning biomass. The only way to protect the spawning biomass will be to control 
fishing level during the whole year. Finally, it has to be said that, in any case, the 
interdiction of fishing in the spawning area and period, or on any other occasion, always 
represents a reduction of the fishing effort. This is not a major inconvenience and in 
some cases, may even be beneficial. 

6. It has to be stressed that no regulation measure will accomplish its objectives without 
observing two conditions : 

The understanding of the fishermen (broadly speaking) that the measure is good for the 
fishery. Hence, it is important to discuss the scientists' conclusions, their objectives, their 
reasons and the expected effects. 
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An efficient fiscalization in the ports and at sea! The 200-mile exclusive zone may be 
very vast and the fiscalization very expensive, but it is not necessary to fiscalize the 
whole area intensively. It is enough to control the areas of larger catches more intensively 
and the remaining areas less so. 

During the last few years, new ways of controlling access to fisheries resources and 
exploitation levels are being implemented. Some examples are the establishment of 
Individual Transferable Quotas systems (1TQ), co-management systems or even the 
system of regional or municipal management, where some management responsibilities 
are attributed to the resource users themselves. 

7. The 1TQ management system is based on the abusive assumption that only the 
economically efficient and profitable vessels deserve to be active in fishing. So, TAC's 
are divided into individual quotas, to be auctioned for the best offer. 

The co-management system delegates a great part of the responsibility of management to 
those who directly exploit the fishing resources - managers, fishermen and their 
professional organizations or unions. With this system neither are quotas sold in auctions, 
nor are the fishing licences lost. 

These are the most well known systems. 


The ITQ system presents the following inconveniences : permanent loss of the titles of 
quotas and of fishing licences; concentration of quotas in the hands of a small group of 
people (who may not even belong to the fishing sector or are even foreigners); and 
underestimation of the social, human and cultural aspects, in favour of economic 
efficiency criteria. 

On the contrary, the co-management system, is concerned with the social aspects of the 
people involved; it seeks their direct and conscientious co-participation with government 
authorities in the management responsibilities, including fiscalization. 
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CHAPTER 6 - PRODUCTION MODELS 


6.1 BASIC ASSUMPTION ABOUT THE EVOLUTION 
OF THE BIOMASS OF A NON EXPLOITED STOCK 

The production models (also called general production models, global models, sintetic models 
or Lotka-Volterra models) consider the stock globally, that is, they do not take into 
consideration the structure of the stock by age or size. 

The total biomass of a non exploited stock cannot grow beyond a certain limit. The value of 
that limit depends, for each resource, on the available space, on the feeding facilities, on the 
competition with other species, etc. In conclusion, it depends on the capacities of the 
ecosystem to maintain the stock. That size limit of the biomass will be designated by 
Carrying Capacity, k. 

The total biomass of a non exploited fishery resource has the tendency to increase with the 
time towards its carrying capacity, k, with a non constant absolute rate, air(B,). The rate, 
air(B t ), is small when the biomass is small, increases when the biomass grows and is again 
small when the biomass gets close to the carrying capacity. Changes, including reductions, 
can occur in the biomass due to fluctuations of the natural factors, but, in any case, the 
tendency will always be an increase towards its carrying capacity. 

The instantaneous rates air(B,) or rir(B,) are therefore not constant. 

In order to formulate the basic assumption of a model for the evolution of the non-cxploited 
biomass, one can adopt a function H of B, , as was done with the basic assumption of the 
individual growth, and define it: 

rir[H(k)-H(B,)] = -r with r constant 


r is the intrinsic growth rate of B,. The relative instantaneous rate, rir(B t ), of the non- 
exploited biomass can therefore be deduced as. 


rir(B t ) 


r[H(k) — H(B, )] 


R* 

dB, 



6.2 EXPLOITED STOCK 

When the stock is exploited, the rate of variation of the biomass due to all causes, that is, the 
total rir(Bi). can be separated into two components: natural rir(B,) due to all causes but not 
fishing and rir(B,) due to fishing: 

rir ( B , L. = rir (B, + rir(B, 
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In an interval of time T, and with a constant fishing level, it will be: 


rir(B, )»>,*, = F, = constant 

and: 

nr (B,L, =f(B,)^ nl 


The natural rate, rir(B c ), (which, according to the basic assumption of the natural evolution of 

the biomass, B, is supposed to be a function of the biomass, B,) is usually designated as f(Bi). 

Comments 

1. Historically, the production models were the first to be used on the analyses of the 
evolution of biological populations, Lotka-Volterra (1925-1928). 

2. Schaefer (1954) applied a production model to a fish stock subject to fishing. 

3. The carrying capacity, k , has been designated in fisheries biology, as B„ and also as B 0 . 
Currently the symbol k is preferred (notice that this symbol is different from the symbols, 
K, of the individual growth models and of the relation S-R). 

4. The basic assumption about natural rir(Bi) previously presented, can be mathematically 
formulated in different ways.. 

5. The production models can only be used in fisheries to analyse the effects of fishing level 
changes and not of changes in the exploitation pattern, because the models consider the 
biomass in a global way and do not take into consideration the age or size stock 
structure. 


6.3 VARIATION OF THE BIOMASS IN THE INTERVAL T i 

The “total”, “natural” and “by fishing” instantaneous rates can be approximated by the 
relative mean rates, rmr(B,). In fact, it can be said that rmr(B t ) s rir(B,) relative to the 
mean B,. (This relation is exact in the case of the exponential model) . 

The following general expression in terms of instantaneous rates 

rir (B,) wJ =f(B l ) mml -F, ridli „ | 


can then be approximated, replacing the rates by the respective mean rates in relation to the 
mean biomass during the interval Tj: 

rmr(B,) = f(Er)-F lnre , at , on ,„ B. 
or 

_L AB . - f fir) F 

7T x — * v ' / 'faiime 

Bj t «ii 


and 


AB, = f(B,).B,.T 1 -F,.B,.T i 
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The variation of the biomass due to all causes of mortality is then decomposed into the 
variation due to natural mortality and the variaton due to the fishing mortality: 

AB.mji = f(Bi).Bi.T i — YidAta, 

The value of the biomass, Bj+i, at the end of the interval T, is: 

B w = B, + B,.f(B i )lT i -Y j 

6.4 LONG-TERM PROJECTIONS (LT) 

(EQUILIBRIUM CONDITIONS) 

The situation of equilibrium at the interval T, implies that the biomass of the stock, at the end 
of the interval Tj (usually 1 year) is equal to the biomass at the beginning of the same interval, 
Bj*i=Bi.or the variation of the biomass is zero ABj=0. 

Bringing the instantaneous rates closer to the mean rates, when the stock is in equilibrium, 
during Tj, then ABj=0 and rmr(B,)=0. Thus, a equilibrium condition will be: 

f.(B“)= F, 

Then the equilibrium conditions, referred to with the subindex E are: 
f(B E )=F E 
y e = f„,b e .t 

6.5 BIOMASS AND FISHING LEVEL INDICES 

In practice the values of B, e F, are not always available and so, one has to look for 
quantities that are associated with the biomass, B, and the fishing level, F, preferably 
quantities (called indices) proportional to those parameters. 


Let U be an index of the mean biomass, B, then during the interval of time T we have: 


U = q.B 


and let f be the index of fishing mortality coefficient, F„ then during the interval of time T: 
f = const.F.T 


from Y = F.B.T and U-q.B 

will have ^ = F.(l /q).U.T 
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thus, to have Y = f .U is necessary to be 



A very common index of B is the catch per unit effort (cpue). The index of F will be the 
fishing effort in an appropriate unit, in order to be proportional to the fishing level. 

The constant of proportionality, q , is designated as the capturability or catchability 
coefficient and indicates the fraction of the biomass that is caught by unit of effort. 


6.6 BIOLOGICAL TARGET REFERENCE POINTS (TRP) 

Long-term (or equilibrium) biological reference points can also be defined for these models. 

Fmsy is the value of F that makes the long-term capture, Y, maximum 

Fmsy is different of Fm,,. In fact Fmsy maximizes the Catch in weight, while Fma*, maximizes 
the Catch in weight per Recruit. Notice that the value of Fm*., cannot be calculated with 
production models, because the age structure of the stock and the recruitment, R are 
considered implicit in the basic assumptions of the model. 

The biological reference points depend on the basic assumptions of the model , therefore the 
value of Fmsy of the structural models is different from the value of Fmsy of the production 
models because the relation S-R, as well as the natural mortality coefficient, M, are implicit in 
the production models. 

To compare results of the two types of models one has to take into consideration that each 
model is based on different basic assumptions. 

For the same reasons, Fo.i of the production models is a different concept to Fo.i of the 
structural models. 

Fo i, Bo i and Yo i of the productions models could be calculated directly from the basic 
assumptions but it is preferable to obtain those characteristics using the constant relations 
between the reference points 0.1 and MSY (Cadima, 1991). 


6.7 TYPES OF PRODUCTION MODELS 

The most common production models in fishery stocks assessment are the Schaefer model 
(1954), the Fox model (1970) and the Pella and Tomlinson model (1969), the latter is also 
designated as GENPROD (name of the computer program that the authors elaborated for the 
application of their model). Fox mentions that the elaboration of his model was based on an 
idea from Garrod ( 1969). 

Each one of these models corresponds to one particular function of H(Bi ) of the basic 
assumption. 
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6.7.1 SCHAEFER MODEL 

The function H(B,) of the basic assumption of this model is: 


Relative instantaneous rate, rir(B,), due to natural causes 

The general basic assumption of the Schaefer model is: 

r ' r .un.ni[k 1 ~ B' 1 ] = -r 

and then, the instantaneous rate of variation of the “natural" biomass can be mathematically 
deduced as: 


fir (®i )n«Mnt 



Equilibrium conditions 

The relative mean rate , rmr(Bt), in relation to B, will be: 


tmr(B,) = f(B E )=r. 

relative to B e 

and, as in equilibrium, ^(®e) = ^e, the equilibrium conditions can then be expressed as: 

B E =k.(l_-F E /r) 

Y E =F E .B E .T, 

Notice that Be is linear with Fe and for Fe =0 , Be = k = carrying capacity = virgin biomass. 

Graphically, the relation between Be and Fe shows a straight line with interception equal to k 
and slope equal to -k/r. 

Target point, F MSV 

The Schaefer equilibrium conditions during one year are: 

Y e = F e .Be B E = k.(l -F E /r) 

Y maximum will occur w hen dY/dF=0, then derivating the previous expression of Ye in order 
to F and making it equal to zero the target point Fmsy will be: 

Target point, Fmsy (Schaefer) 

Fmsy = r/2 BMSY = k/2 YMSY = rk/4 
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In fact, the derivative 


dY/dF = B + F (dB/dF) 

or dY/dF = k(l- F/r) + F (- k/r) = k - 2k.F/r 

and then, Fmsy = r/2 

the relations of the remaining characteristics are obtained by substituting this result in the 
equilibrium conditions. 

Target point, F 0 .i 

The ratio between Fo.i and Fmsy is constant and equal to 0.90, so: 

Target point, Fo.i (Schaefer) 


Fo.i/ Fmsy - 0.90 Bo.i/Bmsy _ 110 Yqi /Ymsy - 0.99 

In fact, as seen before, dY/dF = k - 2k.F/r and, as Fo.i corresponds to dY/dF = 0.1k, so: 


0.1k = k - 2kFo.|/r 


or 0.90 = 2 Fo.i/r 

or 0.90 =F 0 i/Fmsy 

Abundance indices, U, and fishing level indices, f 

As seen in Section 6.5, the indices U and f, are assumed to be proportional to B and F, so 
the equilibrium condition can be written as: 

U E = a + b-fn and Y E = f E . U E (a,b are constants). 

The target point, fiusY. is obtained by equating to zero the derivative of Yn in order to f E : 


Target point, f M sv (Schaefer) 


fMSY - -a/(2b) Umsy - a/2 Ymsy — a‘/(4b) 


In the production models, the ratios fo iTmsy e Uo.i/ Umsy are equal to the ratios Fo.j/Fmsy 
and Bo i / Bmsy- With Schaefer's model we will then have: 

Target point, fo.i (Schaefer) 


fo i/ f M sv = 0.90 U 0 1 / Umsy =1.10 Y 0 i /Y M sy = 0.99 
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From U = q. B, and FT=q.f , the previous expressions of Fmsy and f M sy , one can also obtain 
the relations between the parameters k and r and the coefficients a, b and q: 

k = a/q r = - aq /(bT) kr = - a : /(bT) 

When the value of the interval T is 1 year, T will not appear in these expressions. It is 
possible to calculate the parameters k and r, knowing the values of the capturability 
coefficient, q. Notice that the product k.r does not depend on q. 

6.7.2 FOX MODEL 

For the Fox production model the function H(B t ) will be: 

H(B,) = ln(B,) 

Relative instantaneous rate, rir(B,), due to natural causes 

For the Fox model, from the expression of the general basic assumption, we have: 


rirnamrai [ Ink -InB ] = - r 


and then, as previously referred to, the instantaneous rate of variation of the “natural” biomass 
can be mathematically deduced from that expression and written as: 

=r.ln(k/B.) 

Equilibrium conditions 

The equilibrium condition of the biomass can be expressed by: 
r.ln(k/B E )=F E 

Then, the equilibrium conditions will be: 


B E =k.e- F ‘" Y e = F e .Be.T, 

Notice that ln(Bt ) is linear with F E and that, for F E = 0, B E = k = virgin biomass or carrying 

capacity. The relation between !n( Bi ) and F E is linear, with interception equal to Ink, and 
slope = -1/r. 

Target point, F MSV 

Derivating Y e in order to F and equating the derivative to zero, F M sv, Bmsy and Y M sy will be: 
Fmsy = r BMSY = k/e YMSY = rk/e 
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Target point, F 0 .i 

In this model the ratio between F 0.1 and F MS y is constant and equal to 0.7815. So, it can be 
written: 


Fo,/ Fmsy = 0.7815 B„ , /B MSY = 1.2442 Y 0 , /Ymsy = 0.9724 

These results are obtained in a similar way to those for the Schaefer model. The equation to 
solve will be: 


e~ r * |,r .(l-F 01 /r) = 0.1 

which requires iterative methods to find the value of Fo |/r. The solution is Fo i/r = 0.7815 that 
is igual to Fo i/Fmsy- 


Abundance indices, V , and fishing level indices, f 

For the Fox model the equilibrium condition can be written as: 

Ue =e’ +l>ri or In U E = a + b.f 
and 

Y E =f E . U E 


(a,b are constants) 


Target point, f M sv 

The target point , f M sY , can be obtained by equating to zero the derivative of Ye in order to f E : 
f\isY ~ -lb Umsy = e a /e YMSY = -e a /be 

Target point, f 0 ., 

In the Fox model, the ratios foi/fsisY ct U u / U„ iY are equal to Fo.i/Fmsy and B„, / B« 1Y 
and then: 

fo i/ f MS Y = 0.7815 Uo., / Umsy = 1-2442 Y„ , /Y M sy = 0.9724 


From and from FT=q.f , the following can be deduced: 

k = e a /q. r = - q /(bT) kr = - e a /(bT) 

When the value of the interval T is one year, T will not appear in those expressions. The last 
expression allows the calculation of the product k.r. To calculate k and r separately it is 
necessary to know the value of the coefficient of capturability, q. 
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6.7.3 PELLA AND TOMLINSON MODEL (GENPROD) 

For this production model the function H(B,) will be: 

H(B,)=B' P 

Relative instantaneous rate, rir(B, ), due to natural causes 

The expression of the basic assumption of the GENPROD model, will be : 

- B‘ p )= -r 

therefore 

r * r ( B l)nrtnl = (r/ p)-[l - (k/ B)" 1 ’] 

Equilibrium conditions 

In equilibrium conditions, Fe will be: 

F e = fl B E ) = (r/p).l 1 - (k / B) p ] 

Then, the equilibrium conditions can be expressed as: 

B E =k.[l-p.F E /r]' ,p Y e =F e .B e .T, 


Notice that the relation between ( b e)' and Fe is linear with intercept equal to k and the slope 
equal to -pk/r, in conclusion, for Fe = 0, Be = k = carrying capacity=virgin biomass. 

Target point, F M sy 

Derivating Ye in order to F and equating to zero, we will have: 

Target point, F M sy ( Pella and Tomlinson ) 



i 


l 

1 

_ 


Fmsy ~ r - 

l 

.i+p. 

B msy = k| 

1 

. l + P_ 

p Y = rk 

1 MSY 1 ,K - 

\ 

Li+pJ 

ip ) 


Target point, F 0 .i 

The ratio between Fo i and Fmsy is constant for each value of p and can be obtained in a 
similar way to the previous cases. The equation to solve by iterative methods is: 

X= 1 -0.1.[l-p/(l+p).X] l/p where X=Foi/F M sY 

And also Bo i/ B M sy = [1+p - p./ (I+p).X] <m/w 


Yo.i/ Ymsy - [Fo.i/ Fmsy) [ B 0 1/ Bmsy] 
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The following Table summarizes the most important results: 



parameter p 

Fo.i/Fmsy 

Bo.i/Bmsy 

Yo.i/Ymsy 

Fox 

0.0* 

0.781521 

1.244182 

0.972355 


0.2 

0.819995 

1.193441 

0.978616 


0.4 

0.848355 

1.158613 

0.982915 


0.6 

0.869888 

1.133469 

0.985991 


0.8 

0.886657 

1.114599 

0.988268 

Schaefer 

1.0* 

0.900000 

1.100000 

0.990000 


1.2 

0.910816 

1.088420 

0.991350 


1.4 

0.919724 

1.079045 

0.992424 


1.6 

0.927165 

1.071323 

0.993293 


1.8 

0.933457 

1.064867 

0.994008 


2.0 

0.938835 

1.059401 

0.994602 


2.2 

0.94377 

1.054720 

0.995704 


2.4 

0.947516 

1.050674 

0.995531 


2.6 

0.951059 

1.047146 

0.995898 


2.8 

0.954188 

1.044045 

0.996216 


3.0 

0.956969 

1.041302 

0.996494 


Notice that (Fo.i/Fmsy) + (Bo i/Bmsy) = 2. From this result it can be said that when Fo.i is 
smaller than Fmsy by a certain percentage, the equivalent relation of the biomasses will be 
bigger by the same percentage. 

Abundance indices, U, and fishing level indices, f 

For the Pella and Tomlinson model, the equilibrium conditions can be written as: 

U = (a + b.f) l p or U p = (a + b.f)(a,b are constants). 

Y E =f E . U E 

The target point, f MS Y , can be obtained by equating to zero the derivative of Y E in order to f E : 

Target point, f M sv 


fvisv = -a/(b(l+I/p)) U MSY -[a/(l + p)] ,p Y M sY = -<p/b).(a/(l+p)) (,+, ' p) 

The ratios fo.i/fMSY and (Joy Umsy will be equal to the ratios F'o i/Fmsy and Bo.i/ Bmsy, 
respectively. These last ratios can be observed in the previous Table. 

The values of k, r e kr can also be obtained from U = q. B , and from F.T=q.f 
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k = a ,,p /q r/= - apq /(bT) kr = (-p/bT)a (m/p > 

When the size of the interval is one year, T will not appear in those expressions. The last 
expression in the previous table allows the calculation of the product k.r.The separate values 
of k and r can be calculated if the value of the coefficient of capturability, q is known. 

Comments 

1. The Pella and Tomlinson model has been criticized in its practical application because 
sometimes it produces better adjustments with non reliable values of the parameter p, 
resulting in extremely high values of F M sy- 

2. It is also important to notice that p, the additional parameter of that model, is written with 
different symbols depending on the authors. 

3. The values of the biological reference points relative to F, estimated by Schaefer's model, 
are more restrictive than the corresponding values estimated by the Fox or GENPROD 
production models. 

6.8 SHORT-TERM PROJECTIONS 
6.8.1 GENERAL METHODS 

Long-term projections have been estimated in fisheries since the 50's using these production 
models but in practice, it was only in the 90's that methods were developed for short-term 
projections. These methods are based in the Schaefer, Fox and Pella et Tomlinson expressions 
for the non-exploited biomass. 

By applying production models as referred to in Section 6.3, the variation of the biomass for 
1 year can be expressed, in a general way, as: 


AB iw *sf(Bi).B,.T 1 -Y lu _ i 


or B|.|=B, + B,.f( B, ).Ti - Y, 

where B* = biomass at the beginning of the year i 

B,+i= biomass at the end of the year i 

“i = mean biomass during the year i 
Y; = catch in weight during the year i 

f(B.) 

is the approximation of the mean rate of "natural" variation of the biomass, relative to 
during the year i. 

The expression of the variation of the biomass is the basis for most of the methods for 
short-term projections. Computer programs were prepared for the application of these 
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methods, which also determine long-term projections, biological reference points, etc. Some 
examples are CEDA and BYODIN, Rosenberg et al. (1990), and Punt and Hilbom (1996) 
respectively. 

Theoretically those methods suppose that B, et Yj are known for a period of years. The 
function f( B, ) can be that of Schaefer, Fox or Pella and Tomlinson. 

To determine the parameters r and k it would be necessary to adopt one of the expressions of 
ff Bj) and the value B) in the first interval of the period of years. 

In practice, the values of the annual mean biomasses are not available, only the associated 
quantities, usually assumed to be proportional to tire mean biomasses, that is indices Uj=q. B, 
so, the parameters to be estimated are r, k and q (see Chapter 7). 

The object function, of the least squares method, to be minimized is <I> = £ (U 0 b S - Umod) 2 that 
is, the sum of the squares of the residuals between the observed values and the estimated 
values (designated by error of the process). However, when the the relation U=q. B, is not 
determinant, but it is supposed to have an error, designated as observation error, then it is 
preferable to adopt the object function = £ (lnU 0 t, s - lnU mo d) 2 (Punt and Hilbom, 1996). 

6.8.2 PRAGER METHOD (1994) 

Prager (1994) adopted the Schaefer model and used the relative instantaneous rate of the 
variation of the biomass in the initial basic expression ( not the mean rate approximation) that 
is, 


rirfB,) = r[l - B t /k], B,-F,. B, 


He integrated this expression during the year i and obtained the relation between Bj+i and B, 

He also calculated the mean biomass, B, , integrating B, during the year i. Finally the catch in 
weight is calculated as: 

YrF, Bj. 

The estimation of the parameters can then be made using the least squares method. The 
computer program prepared for this estiamtion is called ASPIC (Prager, 1995). 


6.8.3 YOSH1MOTO AND CLARKE METHOD (1993) 

The short-term projections are derived from the basic assumption of the production models, 
. ri r ( B t ) kiU! — rir(B, ) BUuraJ (F, (fishing 

or, representing rir(B,) „ aIun ,i by f(B,): 

rir ( B .)u»i = f ( B ,)-(F,) fi «h i i„ 
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Integrating this expression during the interval of time T, and considering that: 
rir(B,) = air ln(B t ) and F t =F, =constant: 

li+l 

Jair(ln B.)dt = 

ti 

In (Bi+i) - In (Bi) = f(B,) cum, - F ( Ti = f(B,) T, - F ( Ti 

For the next interval Tj*i (which is the interval where one intends to project the stock and the 
catch): 

In (B 1+2 ) - In (Bj+i) = f(Bi +l ).T,+i - F W .T W 

Calculating the simple arithmetic mean of the two previous expressions and considering that : 

l/2[ln (B K1 ) + In (B, +2 )] - l/2[ln (B,) + In (B M )] = ln( B* it| ) - ln( B*i) 

where B*, is the geometric mean of B, and B,.i, and B*,+i is the geometric mean of Bt+i and 
B, +2> 

Therefore, the mean of the two expressions will be: 

ln( B*i+,) - ln( B*i) = (1/2) { «B,) Ti + f(B,.,) T i+I > - (1/2) {F, T, + F,., T i+ ,} 

The natural rir(Bt) of the Fox model is, as mentioned before f(Bi) = r (Ink - InB,), so the 
approximation, f(B,) can be written as: 

r[lnk- ln( B ,)] 

where B, is the geometric mean of B, and B,.|. 

Therefore, the previous expression relative to the geometric means, can be re-written as: 

In ( BVi)-ln( B*0 = 

= (1/2) { r(lnk- ln( B*0) T, + r(lnk- ln( §**,)) T i+ , } - (1/2) {F, T, + F jt , Th-i } 

To simplify, and as the intervals of T, are usually constant (and equal to one year), one can 
use T instead of Ti and Ti+i and the expression will be: 

ln( BV,) - In ( B*,) = (rT/2) { Ink- In( B*,) + Ink- ln( B* i+I ) }- (1/2) (F, T + Ft*, T( 

or reorganizing the terms of this expression, it will be: 

(1 + rT/2) . In ( B* 1+l ) = rT Ink + (1- rT/2) ln( B*i)) - (T/2) (F, + F M ) 
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Finally, the expression can be written as follows: 

In ( BVO = 2rT/(2+rT) . Ink + (2- rT)/(2+rT) . ln( B*,)) - T/(2+rT) . (F, + F w ) 

As seen in the long-term projections (or equilibrium), it is more common to have biomass 
indices, U, and fishing level indices, f, rather than B et F values. 

Using the indices U = q B and qf = FT 


the Yoshimoto and Clarke expression (1993) can be written as : 


lnUi«i = 


2rT 


,ln(qk) + 


2-rT 


.InUi — 




2 + rT 2 + rT 2+rT 

It is useful, in practice, to write this expression in the following way: 
In U W = b, + b 2 . In Uj+ b 3 . (f + f i+l ) 


where: 


b, 


2rT 
2 + rT 


ln(qk) 


b 2 


2-rT 
2 + rT 


From the coefficients bi, b; and b 3 one can estimate the parameters q, r and k (keep in mind 
that in the long-term projections it was not possible to obtain q separately) as : 

q = -4b3/(l+b2) 

rT = 2(l-b 2 )/(l+b 2 ) 

k = (( 1 +b 2 )/(-4b 3 )).e b,,(l ' bj) 


Comments 

1 . The fact of having developed, in this manual, the Y oshimoto and Clarke model for the 
short-term projections, does not mean a special preference for this model over other 
models for the short-term projections. 

2. Yoshimoto and Clarke designated their expression by the integrated expression of Fox, 
as it is based on the direct integration of the basic assumption. 
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3. Notice that f (BO and f (Bj+i) are, in general, different from f( BO and f( Bj+i). 
However, the means of f( ) may be considered equal to f(means of B) if another type of 
mean of B is used. 

• Definition of Bj’ through a fiinction 

Consider n values Bj and the simple arithmetic mean of f(B,), that is, 

1(B) ={l/n).X f(B0 

• Let B’ be a value such as f( B*) = f(B) . 

B* is designated as the mean of Bj through the function f. 

EXAMPLES 

4. If f(B) = ln(B) then ln( B*) = (l/n).Iln(B,) and B* is the mean of the values B, through 
the logarithm function, also designated as geometric mean of the values B, 

5. If f(B) = B 'then ( B*)' 1 = (l/n).X(B,)' and B* is designated as harmonic mean of the 
values Bj 

6. If f(B) = B' p then ( B*)' p = (l/n).X(B0" p and B* is designated as the mean of order (-p) 
of the values Bj . 

7. Another approach (Cadima & Pinho, 1995) of the integrated equation of Fox can be: 

In Uj+i = bi + b 2 . In U,+ b 3 . (fj + fi+i) 


where: 

b, 

bi= (l-e rT ) ln(qk) k= '~ b ; -e 1 ^ 

- 2 ■ b 3 • In b, 

b 2 = e rT r T= - lnb 2 

b 3 =-q(l-e rT )/(2rT) q = 2-lnb 2 — 

1 — b. 

This last approach of the integrated Fox model can be deduced from the basic assumption of 
the model, during the interval T,: 

nr (B, ),„ ul = r (ln k - In B, )„ mlF „ - (F, )«, 

Taking into account the properties of rates and assuming r, k and F, constant during Tj 
interval, the absolute instantaneous rate of [r.(lnK-lnB,)-Fj] will be: 

air[r.(ln k - lnB t )-Fj] = -r.air (lnBi)=-r.rir(Bi) 

So substituting rir (B,) by the Fox expression mentioned before, one can write: 

air[ r.(ln k - In Bi) - Fj] = -r. [r.(ln k - In B,) - Fj ] 
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or 


air[ r.(ln k - In B, )-F| ] 

[ r.(lnk-ln B,)-F, ] 

Finally, by the definition of rir the expression will be: 

rir[ r.(ln k - In B ( ) - F| ]= -r 

showing that [r.(ln k - In B,)-F,] follows an exponential model, during the interval T, with r 
constant. 

So, the final value of (r.ln k - r.ln B t - Fj), can be expressed as: 

(r.ln k - r.ln B i+ , - F,) = (r.ln k - r.ln Bi - Fi). e’ rTi 
or 

In Bj«i = (1- e' rTi ).ln k + e' rT, .ln B, - (1- e rTi ).F/r 
At the following interval, T,+i, the expression would be: 

In B ]+2 = (1- e rTi+l ).ln k + e‘ rT,+1 .ln B„, - (1- e' rT,+l ).F i+ ,/r 
then, the mean of the two previous expressions, considering T,=7Vi=T, will be: 

In Bj+i” = ( 1 - e rT ).ln k + e' rT In ((1- e rT ) / 2r).(F,+F 1+ i) 
where B,* = geometric mean of Bi and B,-i and Bj+i = geometric mean of Bi+i and B,- 2 . 

Using the indices Uj = q. Bj" and qf = F,T the expression will be: 

In Uh-i = (1- e' rT ).ln (qk) + e' rT . In U.-((l-e rT )/2rT) ,q(f,+f 1+ i) 
which is the initial expression of comment n° 7. 
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CHAPTER 7 - ESTIMATION OF PARAMETERS 

In the previous chapters, several models used in stock assessment were analysed, the 
respective parameters having been defined. In the corresponding exercises, it was not 
necessary to estimate the values of the parameters because they were given. In this chapter, 
several methods of estimating parameters will be analysed. In order to estimate the 
parameters, it is necessary to know’ the sampling theory and statistical inference. 

This manual will use one of the general methods most commonly used in the estimation of 
parameters - the least squares method. In many cases this method uses iterative processes, 
which require the adoption of initial values. Therefore, particular methods will also be 
presented, which obtain estimates close to the real values of the parameters. In many 
situations, these initial estimates also have a practical interest. These methods will be 
illustrated with the estimation of the growth parameters and the S-R stock-recruitment 
relation. 

The least squares method is presented under the forms of Simple linear Regression, multiple 
linear model and non linear models (method of Gauss-Newton). 

Subjects like residual analysis, sampling distribution of the estimators (asymptotic or empiric 
Bookstrap and jacknife), confidence limits and intervals, etc., are important. However, these 
matters would need a more extensive course. 


7.1 SIMPLE LINEAR REGRESSION - LEAST SQUARES METHOD 
Model 

Consider the following variables and parameters: 

Response or dependent variable = Y 

Auxiliary or independent variable = X 

Parameters = A,B 

The response variable is linear with the parameters 

Y = A+BX 


Objective 

The objective of the method is to estimate the parameters of the model, based on the observed 
pairs of values and applying a certain criterium function (the observed pairs of values are 
constituted by selected values of the auxiliary variable and by the corresponding observed 
values of the response variable), that is : 
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Observed values Xi and y, for each pair i, where i=l,2,...,i,...n 

Values to be estimated A and Band (Y|,Y2 Yj,...,Y„) for the n observed pairs of 

values 

A A . 

(Estimates values: A and B (or a and b) and (Y,, Y, Y, Y) 

Object function (or criterium function) 


i=l 


2 


Estimation method 

In the least squares method the estimators are the values of A and B which minimize the 
object function. Thus, one has to calculate the derivatives d<P/dA e 3®/dB, equate them to 
zero and solve the system of equations in A and B. 

The solution of the system can be presented as : , 

x=(l/n).Zx y=(l/n).Iy 

Sxx = I(x-x)(x-x) Sxy = I(x - x )(y - y ) 

b = Sxy/Sxx a = y - b. x 

Notice that the observed values y, for the same set of selected values of X, depend on the 
collected sample. For this reason, the problem of the simple linear regression is usually 
presented in the form : 

y = A + BX + e 

where e is a random variable with expected value equal to zero and variance equal to cr*. 

So, the expected value of y will be Y or A+BX and the variance of y will be equal to the 
variance of E. 

The terms deviation and residual will be used in the following ways: 

Deviation is the difference between y observed and y m „„ ( y ) i.e., deviation = (y- y ) 
while 

A A 

Residual is the difference between y observed and Y cslimaIM i ( Y ), i.e., residual = (y,- Y, ). 

To analyse the adjustment of the model to the observed data, it is necessary to consider the 
following characteristics: 


84 


Copyrighted material 



Sum of squares of the residuals: 


SQ^ =l(y-v ) 2 

This quantity indicates the residual variation of the observed values in relation to the 
estimated values of the response variable of the model, which can be considered as the 
variation of the observed values that is not explained by the model. 

Sum of squares of the deviations of the estimated values of the response variable of the 
model: 


sQ wol «,=l(v-y f 


This quantity indicates the variation of the estimated values of the response variable of the 
model in relation to its mean, that is the variation of the response variable explained by the 
model. 

Total sum of squares of the deviations of the observed values equal to: 

sq^i =X(y~y)’ 

This quantity indicates the total variation of the observed values in relation to the mean 
It is easy to verify the following relation: 

SQiotal “ SQmodcl SQttsidual 


or 


or 

where 


I _ SQ model t SQ residual 

SQlottl SQlota! 

1 =^ + (1 -r 2 ) 


r 2 (coefficient of determination) is the percentage of the total variation that is 
explained by the model and 

l-r : is the percentage of the total variation that is not explained by the model. 
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7.2 MULTIPLE LINEAR REGRESSION - LEAST SQUARES 
METHOD 

Model 

Consider (he following variables and parameters: 

Response or dependent variable = Y 

Auxiliary or independent variables = Xi, X 2 ,..., Xj,..., Xk 

Parameters = B|, B;,..., Bj,..., Bk 

The response variable is linear with the parameters 

Y = B,X|+B 2 X 2 +... + B k X k =X BjX, 


Objective 

The objective of the method is to estimate the parameters of the model, based on the observed 
n sets of values and by applying a certain criterium function (the observed sets of values are 
constituted by selected values of the auxiliary variable and by the corresponding observed 
values of the response variable), that is: 

Observed values xy xy Xj.*, ., x t .i and y-, for each set i, where i= 1 ,2,...,i,...n 
Values to be estimated Bi,B 2 ,...,Bj,...,Bket (Yi,Y 2 ,...,Y, Y n ) 


The estimated values can be represented by : 

Bi, B 2 Bj Bk (oubi,b 2 ,...,bj,...,bk)et Yi, Y: ,..., Y, ,..., Y „ 


Object function (or criterium function) 

<J> = Z ( y ,- Y ,) 2 

i=l 


Estimation method 

In the least squares method the estimators arc the values of Bj which minimize the object 
function. 

As with the simple linear model, the procedure of minimization requires equating the partial 
derivatives of <J> to zero in order to each parameter, B,. where j=l, 2, k. The system is 
preferably solved using matrix calculus. 
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Matrix version 

Matrix X (11 k) = Matrix of the n observed values of each of the k auxiliary 
variables 

Vector y (> „ = Vector of the n observed values of the response variable 

Vector Y (nJI = Vector of the values of the response variable given by the 
model (unknown) 

Vector B (k „ = Vector of the parameters 


Vector B or b (k „ = Vector of the estimators of the parameters 


Mode i 

Y ( n,l) = X (n ,k) . fifk.l) OU Y=X.B+£ 

Object function 


<l><i.i> = (y-Y) T .(y-Y) ou = (y-X.B) T .(y-X.B) 

To calculate the least squares estimators it will suffice to put the derivative dO/dB of <J> in 
order to vector B, equal to zero. d<t>/dB is a vector with components <W>/3B|, 3<t>/dB;, 
a<t>/3Bk. Thus: 

dd>/dB ( k,i) = -2.X T .(y-X.B) = 0 
or X T y - (X T .X). B = 0 

and b = B = (X T .X)'' . X T y 

The results can be written as: 


b ( k. l) = (X T Xy l .X T y 


Y(n.i) = X.b or Y(d-d = X (X T .X)''.X T y 


residuals,„i, = (y-Y) 


Comments 

In statistical analysis it is convenient to write the estimators and the sums of the squares using 
idempotent matrices. Then the idempotent matrices L, (I - L) and (I - M) with 
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L(iui) = X (X T . X)' 1 . X T , I = unity matrix and M (n . n ) = mean (n .i) matrix = 1/n [1] where [1] is a 
matrix with all its elements equal to one, are used. 

It is also important to consider the sampling distributions of the estimators assuming that the 
variables e, are independent and have a normal distribution. 

A summary of the main properties of the expected value and variance of the estimators is 
presented : 


E[ci+C’.u] = Ci+C2.E[u] V[ci+C2.u] = c;.V[u].c: T 


I - Random variable, e 
Expected value of E 
Variance of E 


£„ (independent) 

E[e] = 0. 

V[e] (n .n)=E[E.E T ]=I.CJ 2 


2 - Observed response variable y 
Expected value of y 
Variance of y 


y = Y+e 

E[y] = Y = X.B. 

v ty](n.n) = V[£]( 0 = I.o* 


3 - Estimator of B 

Expected value of B 
Variance of B 


B= (X T .X)'.X T .y 
E[B] = B 

V[B] (kk) =(X T X)-'.cr 


4 - Estimator of Y of the model 

A 

Expected value of Y 


Y = X. B = L.y 
E[Y] = Y. 


Variance of Y 


V[Y] = L.cr 


5 — Residual e 

Expected value of e 
Variance of e 


e = y- Y= (l-L).y 
E(e] = 0 
V[e] = (I-L).a 2 


6 - Sum of squares 

6.1 - Residual Sum of squares = SQ residual, 1 1 j = (y- Y ) T (y- Y ) = y T (I-L)y 

This quantity indicates the residual variation of the observed values in relation to the 
estimated values of the model, that is, the variation not explained by the model. 

6.2 - Sum of squares of the deviation of the model = SQ model ( i i ( = (Y-y) T (Y-y) = y T 
(L-M)y 

This quantity indicates the variation of the estimated response values of the model in relation 
to the mean, that is, the variation explained by the model. 
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6.3 - Total Sum of the squares of the deviations = SQ totalo » = (y-y) T (y-y ) = y T (I-M) y 

This quantity indicates the total variation of the observed values in relation to the mean. 

It is easy to verify the following relation: 

SQiolal — SQmodd "f" SQ rc5 j(juai Or 

[ _ ^QinoJtl t SQrenitol 

SQ,«| SQ„ U , 

or 1=R 2 + (1-R 2 ) 

where: 

R" is the percentage of the total variation that is explained by the model. In matrix terms it 
will be: 


R 2 = [y T (L - M)y].[ (y T (I - M)y]'‘ 

I-R 2 is the percentage of the total variation that is not explained 
by the model. 

The ranks of the matrices (I-L), (I-M) and (L-M) respectively equal to (n-k), (n-1) and (k-1), 
are the degrees of freedom associated with the respective sums of squares. 

7.3 NON-LINEAR MODEL - METHOD OF GAUSS-NEWTON - 
LEAST SQUARES METHOD 

Mode! 

Consider the following variables and parameters: 


Response or dependent variable 

= Y 

Auxiliary or independent variable 

= X 

Parameters 

= B|,B; 


The response variable is non-linear with the parameters 

Y = f(X;B) where B is a vector with the components B,,B: Bj Bk 

Objective 

The objective of the method is to estimate the parameters of the model, based on the n 
observed pairs of values and by applying a certain criterium function (the observed sets of 
values are constituted by selected values of the auxiliary variable and by the corresponding 
observed values of the response variable), that is: 
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Observed values Xi and yj for each pair i, where i=l,2,...,i,...n 

Values to be estimated and (Yi,Y2,...,Yi,...,Y„) form the n pairs of observed 

values. 

(Estimates = Bi, B2,..., Bj,..., Bi,orb|,b 2 ,...,bj,...,bk and Y i, Y 2,..., Y i,..., Y„) 

Object function or criterium function 

<i> = i(y 1 Y i ) 2 

1=1 

Estimation criterium 

The estimators will be the values of Bj for which the object function is minimum. 

(This criterium is called the least squares method). 

Matrix version 

It is convenient to present the problem using matrices. 

So: 


Vector X< 0 .n = Vector of the observed values of the auxiliary variable 

Vector y (t i) = Vector of the observed values of the response variable 

Vector Y(„j)= Vector of the values of the response variable given by the 
model 

Vector B(k,n = Vector of the parameters 

Vector b(k,n = Vector of the estimators of the parameters 


Model 


Y( n ,i) = f{X: B) 


Object function 

<t>(i..) = (y-Y) T (y-Y) 


In the case of the non linear model, it is not easy to solve the system of equations resulting 
from equating the derivative of the function <1> in order to the vector B, to zero. Estimation by 
the least squares method can, based on the Taylor series expansion of function Y, use iterative 
methods. 


Revision of the Taylor series expansion of a function 

Here is an example of the expansion of a function in the Taylor series in the case of a function 
with one variable. 
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The approximation of Taylor means to expand a function Y = ffc) around a selected point, Xo , 
in a power series of x : 

Y = f(x) = fi(xo) +(x-xo).f(xoVl ! + (x-xo) 2 r’(xoV2! +... + (x- xo) 1 f®( xo)/i!+... 
where f*°(xo) = ^derivatives of fix) in order to x, at the point xo- 


The expansion can be approximated to the desired power of x. When the expansion is 
approximated to the power 1 it is called a linear approximation, that is. 


Y = fl[xo) + (x-xo).f (xo) 


The Taylor expansion can be applied to functions with more than one variable. For example, 
for a function Y = fi[X|,X 2 ) of two variables, the linear expansion would be: 


Y = f(x„„,,x 
which may be written, in matrix notation, as 


. , ^ 5f(X| (0 ,,x 2(0 ,) Sf(x„ 0) ,yv 2 ( 0 , 

H0)» A 2(0)/ T ~ A l(0)/* j* T l A 2 a 2(0) /' 


) 


fix. 


8x, 


Y - Y,o)+A(0).(x-X(0)) 


where Y(o> is the value of the function at the point x<o> .with components X| ( oi and x 2 (0) .and 
A(oi is the matrix of derivatives whose elements are equal to the partial derivatives of f(xi,X2) 
in order to X[,X2 at the point (xi ( o), X2(0))- 

To estimate the parameters, the Taylor series expansion of function Y is made in order to the 
parameters B and not to the vector X. 

For example, the linear expansion of Y = ffx,B) in Bi. B2, .... B k , would be: 

Y = ffx;B) = f(x; B«») + (Bj-Bko)) <?f/<?Bi (x;B(oj) + + 

(B 2 -B ; , 0 ))<?f / (?B2 (x;B ( oi) + + + (B|j-Bk ( o,) dt /<?Bk (x;B ( oi) 

or, in matrix notation, it would be : 

Y(„,n = Y(0) (n, If + A(0) ( 0 ,k) • AB| 0 | <k,l) 

where 

A = matrix of order (n,k) of the partial derivatives of the matrix ffx;B) in order 10 the vector 
B at the point B,o> and 


AB ( oi = vector (B - B (0 )). 


Then, the object function will be: 


t> = (y-Y) T (y-Y) = (y-Y (0 ) - A, 0) .AB (0 , ) T (y-Y, 0 , - A(oi .AB(Of) 
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To obtain the minimum of this function it is more convenient to differentiate <t> in order to the 
vector AB than in relation to vector B and put it equal to zero. Thus: 


0 - -2(A(0)) T (y-Y(0)-A ( o) .AB«n) - -2A(0) T (y-Y(0))+ 2 A T Aoi -AB ( o) 

(0) 


or 


A t A(oi .AB(0) - A r (y-Y(0) ) 
(0) (0) 


Therefore: 

A B«,) = ( A (o)- A (o) )" A m (y-Y( 0) ) 

If AB ( o) is “equal to zero” then the estimate of B is equal to B<o). 

(In practice, when we say “equal to zero” in this process, we really mean smaller than the 
approximation vector one has to define beforehand). 

If AB( 0 ) is not “equal to zero” then the vector B ( o) will be replaced by : 

B(i)=B( 0 ) + AB ( o) 

And the process will be repeated, that is, there will be another iteration with B( 0 ) replaced by 
B(i> (and A,oi replaced by A<d ). The iterative process will go on until the convergence at the 
desired level of approximation is reached. 

Comments 

1 . It is not guaranteed that the process always converges. Sometimes it does not, some other 
times it is too slow (even for computers!) and some other times it converges to another 
limit! ! 

2. The above described method is the Gauss-Newton method which is the basis of many 
other methods. Some of those methods introduce modifications in order to obtain a faster 
convergence like the Marquardt method (1963), which is frequently used in fisheries 
research. Other methods use the second order Taylor expansion (Newton-Raphson 
method), looking for a better approximation. Some others, combine the two modifications. 

3. These methods need the calculation of the derivatives of the functions. Some computer 
programs require the introduction of the mathematical expressions of the derivatives, 
while others use sub-routines with numerical approximations of the derivatives. 

4. In fisheries research, there are methods to calculate the initial values of the parameters, for 
example in growth, mortality, selectivity or maturity analyses. 

5. It is important to point out that the convergence of the iterative methods is faster and more 
likely to approach the true limit when the initial value of the vector B«j) is close to the real 
value. 

7.4 ESTIMA HON OF GROWTH PARAMETERS 

The least squares method (non-linear regression) allows the estimation of the parameters K., 
L, and t„ of the individual growth equations. 
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The starting values of K, L» and to for the iterative process of estimation can be obtained by 
simple linear regression using the following methods : 

Ford- Walford (1933-1946) and Gulland and Holt (1959) Methods 

The Ford-Walford and Gulland and Holt expressions, which were presented in Section 3.4, 
are already in their linear form, allowing the estimation of K and L™ with methods of simple 
linear regression on observed Li and Tj, The Gulland and Holt expression allows the 
estimation of K and L. even when the intervals of time T, are not constant. In this case, it is 
convenient to re-write the expression as: 

AL/T, = K.L»-K. L 


Stamatopoulos and Caddy Method (1989) 

These authors also present a method to estimate K, L„ and t Q (or Lo) using the simple linear 
regression. In this case the von Bertalanffy equation should be expressed as a linear relation 
of L, against e' Kl . 

Consider n pairs of values t„ L, where t, is the age and Li the length of the individual i where 
i— 1,2, n. 

The von Bertalanffy equation , in its general form is (as previously seen): 

L» - Lt = (L»- L,). e' K(t '“ ) 


It can be written as: 


L,= U-(L»-L a ). e +Ku . e' Kl 

The equation has the simple linear form, y = a + bx, where: 

y = L, a = L» b - - (L.- L a ). e +Ku 


x = e 


-Kl 


If one takes L a = 0, then t a =to, but, if one considers t a = 0, then L a = L„. 

The parameters to estimate from a and b will be L«, t„ or L 0 . 

The authors propose adopting an initial value K<o), of K, and estimating a<o), b<o) and r 2 <oi by 
simple linear regression between y (= Li) and x(=e k ( 0 )). The procedure may be repealed for 

several values of K, that is, K (1) K^, One can then adopt the regression that results in the 

larger value of r 2 , to which Knux , a™, and b max correspond. From the values of a,™,, b max and 
K max one can obtain the values of the remaining parameters. 

One practical process towards finding K nulx can be: 
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(i) . To select two extreme values of K. which include the required value, for example K.= 0 

and K=2 (for practical difficulties, use K = 0.00001 instead of K = 0). 

(ii) . Calculate the 10 regressions for equally-spaced values of K between those tw r o values 

in regular intervals. 

(iii) . The corresponding 1 0 values of r 2 will allow one to select two new values of K. which 

determine another interval, smaller than the one in (i), containing another maximum 
value of r. 

(iv) . The steps (ii) and (iii) can be repeated until an interval of values of K with the desired 

approximation is obtained. Generally, the steps do not need many repetitions. 

7.5 ESTIMATION OF M - NATURAL MORTALITY COEFFICIENT 

Several methods were proposed to estimate M, and they are based on the association of M 
with other biological parameters of the resource. These methods can produce approximate 
results. 

7.5.1 RELATION OF M WITH THE LONGEVITY, r x 

Longevity: Maximum mean age h of the individuals in a non-exploited population. 

Duration of the exploitable life: t^ — t r = A (Figure 7. 1 ) 


N r = R 



Figure 7.1 Duration of the exploitable life 

Tanaka (1960) proposes "NATURAL" Survival Curves (Figure 7.2) to obtain the values of M 
from longevity. 
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A cohort practically vanishes when only a fraction, p, of the recruited individuals survives. In 
that case, N x = R ■ e' M * , and it can be written: 

p=-^- = e~ MA and so M = -(1/X).lnp 
R 

Different values of the survival fraction produce different survival curves of M in function 
of A. 


2.5 -r Survival curves of Tanaka 



Figure 7.2 Survival curves by Tanaka 

Any value of p can be chosen, for instance, p = 5%, (i.e. one in each twenty recruits survives 
until the age tx) as variable value of the survival curves. 

7.5.2 RELATION BETWEEN M AND GROWTH 
Bever ton and Holt Method (1959) 

Gulland (1969) mentions that Beverton and Holt verified that species with a larger mortality 
rate M also presented larger values of K. Looking for a simple relation between these two 
parameters, they concluded approximately that: 


CM 

VI 

VI 

for small pelagic fishes 

2 < — < 3 

for demersal fishes 


K 

Pauly Method (1980) 

Based on the following considerations: 

1 . Resources with a high mortality rate cannot have a very big maximum size; 

2. In warmer waters, the metabolism is accelerated, so the individuals can grow up to a 
larger size and reach the maximum size faster than in colder waters. 
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Based on data of 175 species, Pauly adjusted multiple linear regressions of transformed 
values of M against the corresponding transformed values of K. , L_ and temperature, T, and 
selected one that was considered to have a better adjustment, that is, the following empirical 
relation: 


In M = -0.0 1 52 - 0.0279 In L_ + 0.6543 In K + 0.463 In T° 


with the parameters expressed in the following units: 

M = year 1 

L_= cm of total length 
K = year 1 

T = surface temperature of the waters in 'C 

Pauly highlights the application of this expression to small pelagic fishes and crustaceans. 
The Pauly relation uses decimal logarithms to present the first coefficient different from the 
value -0.0152 which was given in the previous expression, written with natural logarithms. 


7.5.3 RELATION BETWEEN M AND REPRODUCTION 


Rikhter and Efanov Method (1976) 

These authors analysed the dependency between M and the age of first (or 50 percent) 
maturity. They used data from short, mean and long life species, and suggested the following 
relation of M with the, t™,, age of l“ maturity: 

(Units) 


M = 


1.521 

it v™ 

V l mat50% / 


- 0.155 


• muSCrtii 

M 


> year 
^ year" 


Gundersson Method (1980) 

Based on the assumption that the natural mortality rate should be related to the investment of 
the fish in reproduction, beyond the influence of other factors, Gundersson established several 
relations between M and those factors. 

He proposed, however, the following simple empirical relation, using the Gonadosomatic 
Index (GSI) (estimated for mature females in the spawning period) in order to calculate M: 


M = 4.64„GS1 - 0.37 
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7.5.4 KNOWING THE STOCK AGE STRUCTURE, AT BEGINNING 
AND END OF YEAR, AND CATCHES IN NUMBER, BY AGE, 
DURING THE YEAR 

The natural mortality coefficients M, , at age i can be calculated from the catch, C, , in 
numbers, and the survival numbers, N, and N,+i at the beginning and end of a year, by 
following the steps: 


calculate 

E = ''' 

’ N, - N i+I 

calculate 

Z, =lnN,-lnN itl 

calculate 

M i =Z,-(l-E i ) 


The several values of M obtained in each age could be combined to calculate a constant value, 
M, for all ages. 

Paloheimo Method (1961) 

Let us consider the supposition that F, is proportional to f, for several years i, that is 
F|=q~ for T, = 1 year, F, = q • f, , 

*i 

then: Z, = q ■ f, + M 

So, the linear regression between Z, and f, has a slope b = q and an intercept a = M . 


7.6 ESTIMATION OF Z - TOTAL MORTALITY COEFFICIENT 


There arc several methods of estimating the total mortality coefficient, Z, assumed to be 
constant during a certain interval of ages or years. 

It is convenient to group the methods, according to the basic data, into those using ages or 
those using lengths. 

7.6.1 METHODS USING AGE DATA 

The different methods are based on the general expression of the number of survivors of a 
cohort, at the instant t, submitted to the total mortality, Z, during an interval of time, that is: 

N, = N,.e' z( "° 

Z is supposed to be constant in the interval of time (t.,tb). 

Taking logarithms and re-arranging the terms, the expression will be: 

InN, = Cte - Z.t 

where Cte is a constant ( = In N,+Zt a ). 
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This expression shows that the logarithm of the number of survivors is linear with the age, 
being the slope equal to -Z. 

Any constant expression which does not affect the determination of Z will be referred to as 
Cte. 

1. If Z can be considered constant inside the interval (t a ,tb) and, having available abundance 
data, N, , or indices of abundance in number, U, in several ages, i, then, the application of the 
simple linear regression allows one to estimate the total mortality coefficient Z. 

In fact 

1 - e ‘ ZT ' 

Ni=N,. so Nj = Nj. Constant 

1 ZT, 

and, as 

N i =N I .e' z<, ‘ _, - ) 


then, by substitution: 

N, = Cte.e'* 1 ' (Ti = const = 1 year) 

and also 

InNi =Cte-Zt, 

The simple linear regression between In N, and ti allows the estimation of Z (notice that the 
constant, Cte is different from the previous one. In this case only the slope matters to 
estimate Z). 

2. If ages are not at constant intervals, the expression could be approximated and expressed 
in terms of the Ccntraii. For T, variable, it will be: 

"Nj = N, . e' ZTW 

and, as Nj = N, . e 

it will be Ni = Cte. e Zttcn,ra " 

and finally: In Nj = Cte - Z. t^ntmu 

3. When using indices Uj, the situation is similar because U, = q. Nj, with q constant, and 
then, also: 
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InUi =Cte-Zt. 


The simple linear regression between In U, and t, allows one to estimate Z. 

4. If the intervals are not constant, the expression should be modified to: 

In Uj~ Cte - Z. tc^trai, 

Simple linear regression can be applied to obtain Z, from catches, Ci , and ages, t, , supposing 
that Fj is constant. 

C, = F, N, Tj and so, InC. = Cte + In Nj when Tj is constant. So: 

InCj = Cte - Z. t, 

5. If the intervals are not constant, the expression should be modified to: 

lnCi/T; = Cte - Z. t^,™,, 

6. Let V| be the cumulative catch from ti until the end of the life, then: 

V,= IC k = IF k N Vcum , 

Where the sum goes from the last age until age i. 

As F k and Z k are supposed to be constant LN kcum = N,/Z and so: 

V, = FN/Z and InVi = C“ + InN, 

Therefore: 

In Vj = Cte - Z. ti 

7. Following Beverton and Holt (1956), Z can be expressed as : 





Then, it is possible to estimate Z from the mean age t 

This expression was derived, considering the interval (L, tb) as (t„, ~). 

7.6.2 METHODS USING LENGTH DATA 

When one has available data by length classes instead of by age, the methods previously 
referred to can still be applied. For that purpose, it is convenient to define the relative age. 

Using the von Bertalanffy equation one can obtain the age ! in function of the length, as: 

(the expression is written in the general form in relation to t a and not to to) 
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t = t, - (l/K.).ln[(L_- L, )/( Lx.- U)] 


or 


t = t.-—ln(l- 


L t -L, 
Loo L 


(This equation is referred to by some authors as the inverse von Bertalanffy equation). 

The difference t-ta is called relative age, t , . 

So: t* =-(l/K).ln[(L»- Lt )/( L»- L»)] or t* =-(l/K)ln[l-(L,-L s )/ ( L™- L.)] 

For t, = to , L, = 0 and: 

t" = - ~hi(l-— ) 

K L- 

t" is called a relative age because the absolute ages, t , are related to a constant age, t„. 

In this way, the duration of the interval Tj can either be calculated by the difference of the 
absolute ages or by the difference of the relative ages at the extremes of the interval: 

Tj = tj+i -t, = t*, +l -t', 


Also: 


t*wnlndi tcentrali t - Cte 

r= r+ cte 

So, the previous expressions still hold when the absolute ages are replaced by the relative 
ages: 


In Nj Cte - Z. t ccmrati 


In Uj = Cte-Z. t’cooo .11 

In Vj= Cte - Z. t*j _ 

In C/T = Cte - Z. t’ccmnn, 


Finally, the expression would also be : 



Beverton and Holt (1957) proved that : 
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Z = K-^ 


L-L, 


L must be calculated as the mean of the lengths weighted with abundances (or their 

indices) or with the catches in numbers. 

Comments 

1. The application of any of these methods must be preceeded by the graphical 
representation of the corresponding data, in order to verify if the assumptions of the 
methods are acceptable or not and also to determine the adequate interval, (t a , tb). 

2. These formulas are proved with the indications that were presented, but it is a good 
exercise to develop the demonstrations as they clarify the methods. 

3. It is useful to estimate a constant Z, even when it is not acceptable, because it gives a 
general orientation about the size of the values one can expect. 

4. The methods are sometimes referred to by the names of the authors. For example, the 
expression In V,= Cte - Z.t', is called the Jones and van Zalinge method (1981). 

5. The mean age as well as the mean length in the catch can be calculated from the 
following expressions: 

, _ Z ('central , i -Ci) with Cj = catch in number in the age class i 

ICi 

7“ _ Z(L c entntrCj) 

1C, 

where Q = catch in number in the length class i 

_ Z(* central! ■ C, ) 

ICi 


with Cj = catch in number in the age class. 

The relative age should be t* = - (l/K).ln[(L»- L, )/( L„- L a )] 
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Summary' of the Methods to Estimate the Total Mortality Coefficient, Z 


Assumption: Z is constant in the interval of ages, (t,, t h ) 


T Constant 

InNj =Cte-Z-tj 
In U, = Cte - Z ■ tj 
InC, = Cte - Z • tj 

in V, = Cte - Z ■ t. ("v = £ ^ ^ 

V ' k-ul 


Tj variable 

InN, = Cte-Z-t c , 
In U. = Cte - Z • t . 


central, » = t J L 

central, i 

V ^ . 


In — = Cte - Z • t„ 


In Vj = Cte - Z • t, 
1 


(Beverton and Holt equation of Z) 


Supposition: Z is constant in the interval of lengths, (L„ Lb) 

Relative age tj = — — • lnl 
b ' K 1 

r L--L,] 
^ L_ - La j 



In Ni = Cte - Z.t" 
InN, =Cte-Z-t" , 


InU, = Cte-Zt 


In — = Cte - Z • t’, 


T, = t’„ - 1 " (Gulland and Holt equation) 


In V, = Cte-Z-t, 


L -L 

Z = K 

L-L, 


(Jones and van Zalinge equation) 

(tj, = ■») (Beverton and Holt equation of Z) 
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7.7 ESTIMATION OF THE PARAMETERS OF THE STOCK- 
RECRUITMENT (S-R) RELATION 

The least squares method (non-linear model) can be used to estimate the parameters, a and k, 
of any of the S-R models. 

The initial values of the Beverton and Holt model (1957) can be obtained by re-writing the 
equation as: 


(R/Sy 1 or- = -+— S 
R a ak 


and estimating the simple linear regression between y (= S/R) and x (=S) which will give the 
estimations of 1/a and l/(ak). From these values, it will then be possible to estimate the 
parameters a and k. These values can be considered as the initial values in the application of 
the non-linear model. 

In the Ricker mode! (1954) the parameters can be obtained by re-writing the equation as: 

In — = In a — -S 
S k 

and applying the simple linear regression between y (= In R/S) and x (=S) to estimate lna and 
(-1/k). From these values, it will be possible to estimate the parameters (a and k) of the 
model, which can be considered as the initial values in the application of the non-linear 
model. 

It is useful to represent the graph of y against x in order to verify if the marked points are 
adjustable to a straight line before applying the linear regression in any of these models. 

In the models with the flexible parameter, c, like for example, the Deriso model (1980), the 
equation can be re-written as: 



For a given value of c the linear regression between y (= (R/S) c ) and x ( — S) allows the 
estimation of the parameters a and k. 

One can try several values of c to verify which one will have a better adjustment with the line 
y against x; for example, values of c between - 1 and 1 . 

The values thus obtained for a, k and c, can be considered as initial values in the application 
of the iterative method, to estimate the parameters a , k and c of the non-linear Deriso model. 
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7.8 


ESTIMATION OF THE MATRIX [F] AND OF THE MATRIX [N] - 
COHORT ANALYSIS - AC and LCA 


7.8.1 COHORT ANALYSIS BY AGE- (AC) 

The cohort analysis is a method to estimate the fishing mortality coefficients, Fj, and the 
number of survivors, Nj, at the beginning of each age, from the annual structures of the stock 
catches, in number, over a period of years. 

More specifically, consider a stock where the following is known: 

Data 

age, i, where i = 1,2,. ..,k 
year, j, where j = l,2,...,n 

Matrix of catches [C] with 

C,,j = Annual catch, in number, of the individuals with the age i and during the year j 
Matrix of natural mortality [MJ with 

My = natural mortality coefficient, at the age i and in the year j. 

Vector [T] where 

T( = Size of the age interval i (in general, Tj=T=l year) 

Objective 

To estimate 


and 


matrix [F] 
matrix [N], 


In the resolution of this problem, it is convenient to consider these estimations separately; one 
interval of age i (part 1); all the ages during the life of a cohort (part 2); and finally, all the 
ages and years (part 3). 

Part 1 (interval T,) 

Consider that the following characteristics of a cohort, in an interval Tj are known : 


Cj = Catch in number 


M, = Natural mortality coefficient 
Tj = Size of the interval 


Adopting a value of Fj, it is then possible to estimate the number of survivors at the 
beginning, N, , and at the end, Nj+i , of the interval. 
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In fact, from the expression: 


r = — E — .N .(l-e~ ,F '* M|,T ') 

F, + M, 

one can calculate Nj which is the only unknown variable in the expression. 

To calculate Nj+i one can use the expression N*,, = N 1 .e~ (F|4M '* T ' where the values Ni , F| and 
Mj were previously obtained. 


Part 2 (during the life) 

Suppose now that the catches Q of each age i, of a cohort during its life, the values of M, and 
the sizes of the interval T, are known. 


Adopting a certain value, F Huai , for the Fishing Mortality Coefficient in the last class of ages, 
it is possible, as mentioned in part 1, to estimate all the parameters (related to numbers) in 
that last age group. In this way, one will know the number of survivors at the beginning and 
end of the last age . 

The number at the beginning of that last class of ages, is also the number N^, at the end of the 
previous class, that is, Nami is the initial number of survivors of the class before last. 

Using the C, expression, resulting from the combination of the two expressions above : 


C, = 


F,+M, 


•N r«ui .(e* <F, * M||T ' -1) 


one can estimate F, in the previous class, which is the only unknown variable in the 
expression. The estimation may require iterative methods or trial and error methods. 

Finally, to estimate the number Ni of survivors at the beginning of the class i, the following 
expression can be used : 


Repeating this process for all previous classes, one will successively obtain the parameters in 
all ages, until the first age. 

In the case of a completely caught cohort, the number at the end of the last class is zero and 
the catch C has to be expressed as : 


C' — fin 

^"(F^ + M.) 


•N,, 
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Pope Method 

Pope (1972) presented a simple method to estimate the number of survivors at the beginning 
of each age of the cohort life, starting from the last age. 

It is enough to apply successively in a backward way, the expression: 

Ni= (Ni+i e + C,).e 

Pope indicates that the approximation is good when MT < 0.6 

Pope’s expression is obtained, supposing that the catch is made exactly at the central point of 
the interval Ti (Figure 7.3). 

Ml 


oi 

0 



Figure 7.3 Number of survivors during the interval T| = t| + i - 1| with the catch 
extracted at the central point of the interval 

Proceeding from the end to the beginning one calculates successively: 

N ’* — xr ~ -t-MTi/2 

= N i+1 e 
N’ = N” + Q 
N.-N’.e*" 172 

substituting N’ by N”+Cj, the expression will be: 

N l = (N” + C,).e MT/2 

Finally, substituting N” by Ni+i.e ‘ MTi/2 , it will be: 

N,=(^N lt ,e M ^+C^.e M ^ 


106 


Copyrighted material 



Part 3 (period of years) 


Let us suppose now that the Catch matrix [C], the natural mortality [M] matrix and the vector 
size of the intervals [T], are known for a period of years. 

Let us also assume that the values of F in the last age of all the years represented in the 
matrices and the values of F of all the ages of the last year were adopted. These values will be 
designated by F tcrm ,„ il | (Figure 7.4) 


Years 


Ages 2000 2001 2002 2003 

1 C C C C 

2 C C C C 

3 C C C C 

2h 2h 

1 i i 
III 

F terminal F terminal F terminal « 1 terminal 



Figure 7.4 Matrix of catch, [C], with F, in the last line and in the last 
column of the matrix C. The shadowed zones exemplify the catches 
of a cohort 

Notice that in this matrix the elements of the diagonal correspond to values of the same 

cohort, because one element of a certain age and a certain year will be followed, in the 

diagonal, by the element that is a year older. 

From parts 1 and 2 it will then be possible to estimate successively Fs and Ns for all the 

cohorts present in the catch matrix. 

Comments 

1. The values of M, 0 are considered constant and equal to M, when there is no information 
to adopt other values. 

2. When data is referred to ages, the values Tj will be equal to 1 year. 

3. The last age group of each year is, sometimes grouped agcs(+). The corresponding 

catches are composed of individuals caught during those years, with several ages. So, the 
cumulative values do not belong to the same cohorts, but are survivors of several 
previous cohorts with different recruitments and submitted to different fishing patterns. It 
would not be appropriate to use the catch of a group (+) and to apply cohort analysis. 
Despite this fact, the group (+) is important in order to calculate the annual totals of the 
catches in weight, Y, of total biomasses, B, and the spawning stock biomass. So, it is 

usual to start with the cohort analysis on the age immediately before the group (+) and 

use the group (+) only to calculate the annuals Y, B and (SP). The value of F in that 
group (+) in each year, can be estimated as being the same fishing mortality coefficient 
as the previous age or, in some cases, as being a reasonable value in relation to the values 
of F, in the year that is being considered. 

4. A difficulty in the technical application appears when the number of ages is small or 
when the years are few. In fact, in those cases, the cohorts have few age classes 
represented in the Matrix [C] and the estimations will be very dependent on the adopted 

Values of Ftcnninals* 


107 


Copyrighted material 







5. The cohort analysis (CA) has also been designated as: VPA (Virtual Population 
Analysis), Derzhavin method, Murphy method, Gulland method, Pope method, 
Sequential Analysis, etc. Sometimes, CA is referred to when the Pope formula and the 
VPA are used in other cases. Megrey (1989) presents a very complete revision about the 
cohort analyses. 

6. It is also possible to estimate the remaining parameters in an age related to numbers, 
that is, N cum „ Ni, D„ Z, and E,. When the information on initial individual or mean 
weights matrices [w] or [ w) are available, one can also calculate the matrices of annual 
catch in weight [Y], of biomasses at the beginning of the years, [B|, and of mean 
biomasses during the years [ B], If one has information on maturity ogives in each year, 
for example at the beginning of the year, spawning biomasses [SP] can also be 
calculated. Usually, only the total catches Y, the stock biomasses (total and spawning) at 
the beginning and the mean biomasses of the stock (total and spawning) in each year are 
estimated. 

7. The elements on the first line of the matrix [N] can be considered estimates of the 
recruitment to the fishery in each year. 

8. The fact that the F tcm ,j M is are adopted and that these values have influence on the 
resulting matrix [F] and matrix [N], forces the selection of values of F tcm ,, M is to be near 
the real ones. The agreement between the estimations of the parameters mentioned in the 
points 6. and 7. and other independent data or indices (for example, estimations by 
acoustic methods of recruitment or biomasses, estimations of abundance indices or 
epue's, of fishing efforts, etc) must be analysed. 

9. The hypothesis that the exploitation pattern is constant from year to year, means that the 
fishing level and the exploitation pattern can be separated, or F KP i = Fj x Sj . This 
hypothesis can be tested based on the matrix [ F ] obtained from the cohort analysis. 

It is usual to call this separation VPA-Separable (SVPA). 


We have 

i 

and 

£ f m = s ™, 

j 

and 

EFj.i = F.o. 




Then, if F^ = F r s, one can prove that F r s, = (F^. .s lol )/F ra . 

If the estimated values of F :J are the same as the previous Fsep,, = Fj.s, then the hypothesis 
is verified. This comparison can be carried out in two different ways, the simplest is to 
calculate the quotients F sep„ /Fy. If the hypothesis is true this quotient is equal to one. If 
the hypothesis is not verified it is always possible to consider other hypotheses with the 
annual vector [s] constant in some years only, mainly the last years. 

10. It is usual to consider an interval of ages, where it can be assumed that the individuals 
caught are “completely recruited”. In that case, the interval of ages corresponds to 
exploitation pattern constant (for the remaining ages, not completely recruited, the 
exploitation pattern should be smaller). For that interval of ages, the means of the values 
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of F.j in each year are then calculated. Those means, Fj, are considered as fishing levels 
in the respective years. The exploitation pattern in each cell, would then be the ratio F,j 
/ F, . The Si , for the period of years considered, can be taken as the mean of the relative 
pattern of exploitation calculated before. Alternatively, they can also be taken as referring 
to s,of an age chosen for reference. 

7.8.2 LENGTH COHORT ANALYSIS - (LCA) 

The technique of the cohort ansalysis, applied to the structure of the catches of a cohort 
during its life, can be made with non constant intervals of time, Tj,. This means that the length 
classes structure of the catches of a cohort during its life, can also be analysed. 

The methods of analysis of the cohort in those cases is called the LCA (Length Cohort 
Analysis). The same techniques; Pope method, iterative method, etc, of the CA for the ages, 
can be applied to the LCA analysis (the intervals T/s can be calculated from the relative 
ages). 

One way to apply the LCA to the length annual catch compositions, will be: to group the 
catches of length classes belonging to the same age interval in each year. The technique CA 
can then be applied directly to the resulting age composition of the catches by age of the 
matrix [C]. This technique is known as “slicing” the length compositions. To “slice”, one 
usually inverts the von Bertalanffy length growth equation and estimates the age t, for each 
length L, (sometimes using the relative ages t *) (Figure 7.5). It is possible that when grouping 
the length classes of the respective age interval, there are length classes composed by 
elements that belong to two consecutive age groups. In these cases, it will be necessary to 
“break” the catch of these extreme classes into two parts and distribute them to each of those 
ages. In the example of Figure 7.5, the catches of the length class (24-26] belong to age 0 and 
to age 1 . So, it is necessary to distribute that catch to the two ages. One simple method is to 
attribute to age 0 the fraction (1.00 - 0.98)/(1.06 - 0.98) = 0.25 of the annual catch of that 
length class and to age 1 the fraction (1.06 - 1.00)/(1.06 - 0.98) = 0.75. The method may not 
be the most appropriate one, because it is based on the assumption that, in the length classes, 
the distribution of the individuals by length is uniform. So, it is necessary to use the smallest 
possible interval of length classes, when applying this distribution technique. 

Another way to do the length cohort analysis is to use the catches in the length classes of the 
same age group. It is possible to follow the cohorts in the matrix [C], through the length 
classes belonging to a same age, in a certain year, with the length classes of the next age, in 
the following year, etc. In this way, the different cohorts existing in the matrix will be 
separated and the evolution of each one of them will be by length classes, not by age (see 
Figure 7.5). 
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2003 


Group Years 

Age Relative Length 2000 2001 2002 

Classes 



1.03 

20- 

41 

30 

17 

49 

0 

1.54 

22- 

400 

292 

166 

472 


1.98 

24- 

952 

699 

400 

1127 


2.06 

26- 

1766 

1317 

757 

2108 

1 

2.30 

28- 

2222 

1702 

985 

2688 


2.74 

30- 

2357 

1872 

1093 

2902 


2.88 

32- 

2175 

1091 

1067 

2739 


3.00 

34- 

1817 

948 

1416 

1445 


3.42 

36- 

1529 

812 

1270 

1250 

2 

3.64 

38- 

1251 

684 

980 

1053 


3.83 

40- 

1003 

560 

702 

710 


3.96 

42- 

787 

290 

310 

558 

3 

4.01 

44- 

595 

226 

179 

834 


4.25 

46- 

168 

70 

71 

112 


Cohort of the year 2000 

Figure 7.5 Example of a matrix [C] with the catches of the cohort shadowed, 
written in bold, recruited at year 2000, “sliced” by length classes, 

The LCA R. Jones method (1961), of analysing a length composition during the life of a 
cohort can then be applied. The different values of Tj are calculated as T, = t, t |*-t,’\ where t,* 
e tj+i* are the relative ages corresponding to the extremes of the length interval i. The vector 
[N] can also be obtained as the number of initial survivors in each length class of the cohort, 
and in each age class. 

Comments on cohort analyses 

1. Certain models, called integrated models, combine all the available information (catches, 
data collected on research cruises, effort and epue data, etc) with the matrix [C], and 
integrate in a unique model, in order to optimize the previously defined critcrium 
function. A model integrating CA and the hypothesis of constant exploitation pattern was 
developed and called SVPA, separable VPA, because the Fishing level and Exploitation 
pattern are “separable”. 

2. Fry (1949) considered the cumulative catches of a cohort by age during its life, from the 
end to the beginning, as an image of the number of survivors at the beginning of each 
age (which the author designated as “virtual population”): 

£ C k = V, = N, virtual 

ksfinal 

no 
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In the fishery that Fry studied, M was practically equal to zero. 

If M is different from zero it can also be said that the number Nj of survivors at the 
beginning of the interval i will be : 


where Dk represents the number of total deaths at the interval k. 


Adopting the initial values, Eyo) , for the exploitation rates, E , in all the classes, one can 
calculate the total deaths: 

Dk<0) = Ck/Ek<oy 


N j(0) can be calculated as the cumulative total deaths from the last class up to the i lh class, 
that is: 


N„„, = I I 

k=uh 


Then the expression will be: 


= ln(N lM(0| / N 1(0| ) 


and: 


p t = p 7 T 
r i(l) i l 'l(0) z 'i(l)' , i 


Comparing E 1( ij with Ej((», the new values of E will be: 


One can then estimate the values of E with the desired approximation by an iterative 
method, repeating the five calculations (of Dj, Nj, ZjT„ F,T, and E„) using E,(n instead of 
Ej(o). 


In the last class, the number, Num, can be taken as equal to the number of deaths. D| W , 
and in this case, Ni„, will be calculated as : 

N|a(l = = C) a3I /Ejaj, 
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3. Finally, the results of CA and of LCA give a perspective view of the stock in the 
previous years. That information is useful for the short and long-term projections. 
Usually, data concerning the catches is not available for the year in which the assessment 
is done and so it is necessary to project the catches and the biomasses to the beginning of 
the present year before calculating the short-term projection. 

4. When the relative ages are calculated, it is usual to adopt zero as the age t a 
corresponding to the value of L, , taken as the lower limit of the first length class 
represented in the catches. 
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CHAPTER 8 - EXERCISES 


8.1 MATHEMATICAL REVISION 

1. Calculate: 


A) 

10 4 

16 2 

8427° 

0.0 1 0J 

B) 

4 3 4 

4 5 

I 3 ")' 

5 2 +4 2 

2 2 x2 

C) 

log 1000 

log 0.01 

lo g(~T 
V. 10 3 ) 


D) 

Ine 

In- 

e 

Ine" 5 

e lnf 

E) 

0 A 

0 °° 

A 00 

0 


0 o A ~ “ A “ 0 


2. Verify that 

a) a = e' na b) a = lo' OB ° 

e x - 1 X _ , £ 

c) = 1 for -0.01 <x < +0.01 d) - = e 2 for -0.5<x<+0.5 

x x 

3. Solve the following expressions applying natural logarithms to both members of the 
equality: 

a) y-a x^ b) y = a ■ t+2 ' c ) c) y-a = b e c ^ x 

Note: a, b, e c are constants; e is the basis of natural logarithms (e = 2.7183..) ; x and y are 
variables. 

4. Determine the value of .v in the following expressions: 

a) e~ x = 5.2 b) 10* =5.5 c ) y-a = be c{x ~ h] 
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5. Calculate the derivatives of the following expressions: 


a) 

,v=13 

g) 

y=5 x 

m) 

y^4+2xf 

b) 

y=3-8x 

h) 


n) 

y=(x-6) 2 

c) 

>^ 5 

01 

>i=lnx 

o) 

y=a.(3-e b *) } 

d) 

y=x 2n 

j) 

>>=ln(5x+4) 

P) 

y= (4x+3).(<?'-4) 

e) 

y=S 

k) 

y= l/x 


El 


1) 

y= (2+4x)/(3-x) 



6. Calculate the indefinite integrals of the following functions: 


a) 

f(x) = 0 

0 

f( X ) = I 

x 

k) 

f(x)=e-° 3 ' 

b) 

f(x) = 5.34 

g) 

f ( x )=r— 

2-5-x 

1) 

f(x) = 3e 1 ** 1 

c) 

f(x)=x 6 

h) 

f(x) = - 1 + 2 x ,■ 
40 + x + x 2 

m) 

f(x)=x-e* 

d) 

f(x)= 1 + 3x 

i) 

f( x ) = e * 

n) 

f(x) = lnx 

e) 

f(x) = 4-x* 3 

j) 

f (x) = e 02 * 

o) 

f(x)= x ■ lnx 


7. Calculate the area under the function 

a) f(x) = 2 + 5x between x = 1 and x = 4 

b) f(x) = e 3 '* between x = 0 and x = 1 

c) ffx) = — - — between x = — and x = — 

5 + 2x 2 2 

d) f(x) = 1 + 3x between x = -2andx = 2 

8. Calculate the value of y cum uiaiivc with 

a) y- e between x = 0 an d x = 0.8 
f(x) — — — — 

b) 1 + 2x between x - 0 and x = 2 

c) ^ x ) _ ^' x between x = 0 an d x = 1 
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9. Calculate the Mean Value of y with 


a) y = 3 e 7s between* = 0 and x = 1 

b) y = 4 ■ (l -e" 0,2 * ) between *=1 and * = 3 

c) y = 2 - x between * = 0 and * = 1.2 

1 0. Calculate the integral of 

a) f(x) = 2 5 x with the initial condition x = 1=> F(x) = 4whereF(x) = (f(x)dx 

b) f(x) = — with the initial condition F(l) = 2 

x 

c) — = 0.2 • y with the initial condition * = 0 => y = 1 0 
dx 

dy 3 

d) — = with the initial condition x = 0 => y = 0 

dx l + 3x 

8.2 RATES (2.2) 

Consider the function, y = 40 -35.e _0 ' x at the interval (0,10) 

1. Calculate: 

a) The values of y for * = 0,1,2, 3, 4,5,6, 7,8,9,1 0 ; 

b) Represent graphically the function y at the interval (0,10) of*; 

c) The variation. Ay , corresponding to the interval (l,2) of*; 

d) The absolute mean rate of variation of v, amr(y), at the intervals (1,7), (2,5), (5,6) and 
(8,9) of*; 

e) The absolute instantaneous rate of variation of_y, air(y), at the points x=3 and x=4; 

f) Calculate the relative mean rate of variation of y, r.m.r.fy), at the interval (8,9) in 
relation to the value of y corresponding to the initial point, to the final point and to the 
central point of that interval; 

g) Calculate a relative instantaneous rate of variation of y, r.i.r.(y) at the central point of 
the interval (8,9). 
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2. Calculate the air(y) of the following functions: 

a) y = 1 + 1 Ox 

b) y = x 3 - 2x + 3 

c) y = e x 

d) y = lnx 

3. Calculate the rir(y) of the following functions: 

a) y = 4 + x 

b) y = e K 

c) y = 6-e : ‘ 

d) y = a • x with a = constant 

4. Calculate the air of the air(y) of y = 3x ' - 4x - 1 2 

5. Given the function, y = 3 e~ “ verify that rir(y)=air(ln y) 

8.3 SIMPLE LINEAR MODEL (2.3) 

Consider a model that relates the characteristic y with time l, where the basic assumption is: 
air(y)= -3, for 0 < t < «*> 

Adopt the initial conditionl: for t = 0 , y=30 

1 . How would you designate this model? 

Write the general expression for the value of the characteristic y at the instant t; 

2. Calculate the value ofy when t = 0,1,2, 3, 4,5, 6 and draw the graph ofy against I. 

3. Considering the interval of time, A t, from t = 2 to t = 4 

a) Calculate the variation of y during the mentioned interval At; 

b) Calculate the central value ofy in the interval At; 

c) Calculate the cumulative value ofy in that interval, y^ 

d) Calculate the mean value, y , of y, in the interval At; 
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e) Calculate the simple arithmetic mean ofy in the interval A/; 

f) Verify that the arithmetic mean ofy is equal to the mean value, y , and equal to the 
central value, y ccnM |, ofy in that interval. 

g) Verify that in the linear model, the amr(y) = air(y) = constant. To do that, calculate, 
for the above mentioned interval. At, the amr(y) and the air(y) and compare the 
results. 

Repeat exercise 3. considering the interval from / = 0 to / = 10. 


8.4 EXPONENTIAL MODEL (2.4) 

Consider a model that relates the characteristic y with time t, through the following basic 
assumption: 

rir(y) = -0.4 for 0 < f < <» 

Adopt the initial condition : for t = 0, y =100 

1 . Write the general expression for the value of the characteristic y at the instant t; 

a) Calculate the value of y at the instants t =1,2, 3, 4, 5, 6. 

b) Represent, graphically, the values of y calculated above, against the corresponding 
values of t. 

c) Represent, graphically, the values of lny against the given values of t. 

2. Considering the interval of time At = (3,6) 

a) Calculate the variation of y, Ay, during the interval At. 

b) Calculate y ccnua i in the interval At. 

c) Calculate the value of y cum in the interval At. 

d) Calculate y in the interval At. 

e) Show that the geometric mean of the values of y for / = 3 y / = 6 is equal to y C cnirai 
and approximately equal to yin that interval. 

f) Show that, in that interval, rmr(y) reWm ^ = rir(y) = -0.4 

3. Consider the interval of time from t = 0 to t = 1 0 . Repeat the calculations of questions 2 
item a), c) and d) for this interval. 
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8.5 COHORT - EVOLUTION IN NUMBER 


(3.2) 


GROUP I 


From the stock of megrim, Lepidorhombus whiffiagonis, from Divisions VIIIc and IXa of the 
ICES, the Assessment Working Group of ICES (ICES, 1997a) estimated that the fish recruit 
to the exploitable phase, at the beginning of age 1 year, and that in 1996 the instantaneous 
rate of total mortality during the exploitable phase was 0.7 year 1 . 

Consider 1000 individuals of a cohort of megrim, recruited to the exploited phase, which 
starts at the beginning of age I and finishes at the end of age 7 years. 


a) What is the value of the rir of the variation of N, in this interval ? 

b) What is the value of the rir of the mortality of N, in this interval ? 

c) Calculate the annual rate of survival during the interval. 

d) Calculate the annual rate of mortality during the interval. 

a) Calculate the number of survivors at the beginning of each age of that interval. 

b) Calculate the number of survivors at the end of 7 years of age. 

c) Draw the graph of the number of survivors in each age of that interval. 

d) Calculate the number of deaths in each age of the interval. 

e) Calculate the number of deaths during all the exploitable phase. 

f) Determine the percentage of the initial number of 3 year-olds that survive until the 
beginning of their 6th year. 

g) Determine the percentage of the initial number of 3 year-olds that die before the 
beginning of their 6th year. 

h) Calculate the mean number of survivors during each age of the given interval. 

i) Calculate the cumulative number of survivors during ages 3 to 5. 

j) Calculate the mean number of survivors between the beginnings of ages 3 and 6. 
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GROUP II 


The Working Group of ICES which evaluated the stock of ibcric sardine, Sardina pilchardus , 
estimated the mortality rates in each age of 1995 (ICES, 1997b). presented in the following 
table : 


Age Group 

1) 

1 

2 

3 

4 

5 

6 

Annual Rate of 
Mortality 

0.36 

0.43 

0.54 

0.63 

0.66 

0,68 



Suppose that these rates correspond to a cohort. 

1 . As previously seen, there may be several types of rates (ex: amr, rmr, air and rir and the 
relative rates were referred to several values of the characteristics). What type of rate is 
the annual rate of mortality? 

2. Calculate the survival rate in each age class. 

3. Calculate the total mortality coefficient for each age class. 

4. Calculate the survival rate between the beginning of age 1 and the end of age 4. 

5. Calculate the annual mean rate of survival in the same interval of ages. 

GROUP III 

Consider a cohort of a certain species for which the number of survivors at the beginning of 
age 2 years is 4325, while the number of survivors at the end of age 2 years is 2040. 

1 . Calculate the mean number of 2 year-old individuals and the number of individuals at the 
age of 2.5 years. Compare the results. 

2. If the annual rate of mortality of this cohort during the ages 3 and 4 is 70 and 60 
percent respectively, calculate the percentage of the initial number of individuals at age 3 that 
will survive until the end of age 4. 

Give the relation between the survival rate during the period that covers the ages 3 and 4 
years and the annual survival rates of ages 3 and 4 years. 

8.6 COHORT - CATCH IN NUMBER (3.3) 

GROUP I 

According to the Assessment Working Group of ICES (ICES, 1996a) the relative 
instantaneous rates of total and natural mortality for the age of 3 years, of the stock of blue 
whiting, Micromesistius poutassou, in 1995, were estimated as being, respectively, 
Z 3 =0.4 year ' and M , =0.2 year' 1 . In that year, the number of survivors at the beginning of 
age 3 year was 2600 million individuals. 
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1 . Calculate, for the cohort of 1992 and for the age of 3 years : 

a) The total annual rates of survival and of mortality. 

b) The relative instantaneous rate of fishing mortality. 

c) The exploitation rate. 

d) The number of deaths during the age. 

e) The mean number of survivors during the age. 

f) The total catch in number of 3 years old individuals. 

g) The number of survivors at the end of the age. 

GROUP II 

The 4 year-old age group of the stock (Div. ICES Vlle-h) of whiting, Merlangius merlangus 
merlangus, is simultaneously exploited by the crustaceans trawl fleet and demersal fish trawl 
fleet. 

The Working Group of ICES that evaluates this stock estimated (ICES, 1996a) that, in 1995, 
the total instantaneous mortality rate of that age (4 years) was Z = 1.3 year '. Suppose, for this 
exercise, that the instantaneous rate of fishing mortality caused by the crustaceans trawl fleet 
was Fc = 0.5 year' 1 , while the corresponding value for the demersal fish trawl fleet was Ff = 

0. 6 year' 1 . The Group also considered the natural instantaneous mortality rate, M = 0.2 year' 1 . 

1 . In 1 995 1 7.66 million individuals recruited at 4 years of age. 

a) Calculate the total number of deaths during that age. 

b) Calculate the mean number of survivors during the age. 

c) Calculate the exploitation rate of each fleet. 

d) Calculate the total exploitation rate. 

e) Calculate the catch in number by each fleet and the total catch in number. 

f) Calculate the number of survivors at the end of the age. (The solutions to the questions 
can be done in any order). 
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GROUP III 


For this exercise, suppose that in 1990, the mean number of survivors of the cohort of a stock 
of anchovy, Engrautis encrasicholus, during the period of 2 years of age was calculated as 
about 50 million individuals. During 1992, 70 million individuals were caught, from which, 
40 percent were caught by the national fleet, and it is estimated that 80 million died of natural 
causes. 

1 . For this age and this cohort: 

a) Calculate the total exploitation rate and the exploitation rate of the national and 
foreign fleets. 

b) Calculate the total, natural and fishing mortality coefficients. 

c) Calculate the fishing instantaneous mortality rates caused by the national fleet and by 
the foreign fleet. 

d) Calculate the number of survivors at the beginning and at the end of the age. 
(The solutions to the questions can be done in any order). 

GROUP IV 

According to the Assessment Working Group of ICES (ICES, 1997a) the fishing mortality 
coefficients applied to the 1976 cohort of the stock of common sole, Solea vulgaris, of the 
Celtic Sea were estimated in each age from 2 to 8 years (following table). The natural 
mortality coefficient for this stock is considered constant and equal to 0.1 year 1 . It was 
estimated that, at the beginning of age 6 years there were 1112 million survivors of this 
cohort. 


Age 

2 

3 

4 

5 

6 

7 

8 

Fi 

0.07 

0.22 

0.33 

0.41 

0.45 

0.41 



1. Calculate the numbers of survivors of this cohort at the beginning of each of the above 
mentioned ages. 

2. Calculate the number of deaths in each age referred in the table. 

3. Calculate the exploitation rates of this cohort in each age. 

4. Calculate the mean numbers of survivors during each of the above ages. 

5. Calculate the catches, in numbers, extracted from this cohort during each of the ages 
mentioned above, using two different methods. 


8.7 INDIVIDUAL GROWTH IN LENGTH AND WEIGHT (3.4) 

GROUP I 

The parameters of the von Bertalanfly length growth equation of the stock of European 
anglerfish (Div. VIIIc and IXa of ICES), Lophius budegassa, were estimated as (Duarte el a/., 
1997): 
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Asymptotic length = 101.69 cm 
Coefficient of growth in length = 0.08 year' 1 
Theoretical age, when the length is zero = -0.2 year 

1. Calculate the theoretical length corresponding to the age 3.84 years. 

2. Calculate the length at the beginning of the ages 1 to 12 years. 

3. Calculate, for each of the above mentioned ages, the central length. 

4. Represent, graphically, the Bertalanffy curve of growth in length for this stock. 

GROUP II 

Using the growth parameters given in Group I: 

1. Calculate the length that corresponds to each age interval between 1 and 12 years as 
being the simple arithmetic mean of the length at the beginning and at the end of each 
class. 

2. Calculate the mean length in each age for the same interval from 1 to 12 years 
accordingly to the von Bertalanffy model. 

3. Compare the lengths obtained in 1) with those obtained in 2) and with the central values 
in each age of the interval, calculated in Group 1-3. 

GROUP III 

The data presented in the following table represents the mean length (cm) by age (years) 
obtained from direct age reading of individuals of the stock of European anglerfish, Lophius 
budegassa, (Div. VUlc and IXa). 


t 

L,(cm) 

T 

Li(cm) 

1 9.2 

7 

44.4 

2 

16.5 

8 

49.0 

3 

22.9 

9 

52.3 

4 

28.8 

10 

55.0 

5 

34.7 

11 

60.8 

6 

38.6 

12 

63.4 


Based on this data the parameters of the growth equation were estimated according to the 
Gompertz model as being : 

Gompertz: L_ = 73.7 cm; k = 0.22 year 1 t*=- 2.76 year 

(t* is the age corresponding to L=1 cm) 


122 


Copyrighted material 




1 . Represent, graphically, the observed values. 

2. Calculate, for the interval 1-12 years, the values of the length at the beginning of each age, 
according to the Bcrtalanffy growth model and draw the corresponding growth curve. 

3. Calculate, for the interval 1-12 years, the values of the length at the beginning of each age, 
according to the Gompertz growth model and draw the corresponding growth curve. 
Determine the inflection point of the curve. 

4. Say which growth model you consider more appropriate for this case and justify your 
answer. 


GROUP IV 

The data presented on the following table concern the stock of European anglerfish Lophius 
budegassa (Div. VIIIc and IXa). 

Table of individual weights by length class, taken from the samples of European anglerfish, 
Lophius budegassa collected by IEO and 1P1MAR in 1994. 


HUS3B 

cm 

n 


mm 

1) 

20- 

129 

3 

50- 

1685 


22- 

163 

2 

52- 

1896 


24- 

219 

4 

54- 

2107 

24 

26- 

265 

14 

56- 

2345 

41 

28- 

320 

8 

58- 

2569 

41 

30- 

397 

10 

60- 

2848 

32 

32- 

486 

9 

62- 

3126 

35 

34- 

545 

57 

64- 

3407 

28 

36- 

664 

60 

66- 

3700 

19 

38- 

773 

61 

68- 

4056 

23 

40- 

890 

58 

70- 

4411 

17 

42- 

1027 

64 

72- 

4764 

13 

44- 

1122 

56 

74- 

5203 

8 

46- 

1334 

50 

76- 

5587 

4 

48- 

1503 

37 

78- 

5982 

3 


The mean of the observed weights and the number (n) of individuals, are indicated for each 
length class. 

Based on this data, the parameters of the weight-length relation were estimated for this stock 
as : 

a = 0.021 
b = 2.88 

1 . Calculate the theoretical weight for each length class. 

2. On a graph, plot the observed and the theoretical weights against the length classes. 
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3. Suppose you want to use the weight-length relation, with b=3 (the constant estimated for 
this relation being a=0.013). Calculate, for this case, the theoretical weights for each 
length class. Compare these values with the theoretical weights estimated in 1). 

4. Using the results obtained up to now, write the Bertalanffy growth equation, in weight, 
for this stock. 


8.8 COHORT DURING ALL LIFE - BIOMASS AND CATCH IN (3.6) 
WEIGHT 


GROUP I 

The recruitment to the exploitable phase of horse-mackerel, Trachurus trachurus , distributed 
in the Ibero-Atlantic waters (Div. VIlIc and IXa) occurs at age 1. In order to make the 
calculations, let us consider the exploitable phase between the ages 1 and 10 years and the 
recruitment of a cohort equal to 1 000 individuals. 

The parameters of the von Bertalanffy equation were estimated by the Working Group of 
ICES (ICES, 1998a), using the mean lengths, at age of the catch, as being: 

L_ = 34.46 cm 
K = 0.225 year 1 
to = -1.66 year 

The weight-length relation was also estimated, using the mean weights by age adopted by the 
WG for the long-term projections (ICES, 1998a), as being: 

W(g) = 0.01 1 L(cm) 190 

The mortality of this stock is characterized as : 

• A constant natural mortality coefficient during all the exploitable phase : M = 0. 1 5 
year 1 . 

• Fishing mortality coefficients in 1996 (ICES, 1998a) are different with age as : 


Age (year) 

1 

2 


4 

5 

6 

7 

8 

9 

10 

Fi (year 1 ) 

0.24 

0.26 



0.06 

0.09 

0.12 



0.24 


1. Organize the calculations in a spreadsheet, in order to calculate, for each age of the life 
of the cohort (from age 1 to 10): 

a) The number at the beginning of the age 

b) The individual weight at the beginning of the age 

c) The biomass at the beginning of the age 

d) The number of deaths during the age 
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e) The mean number of survivors during the age 

f) The mean individual weight during the age 

g) The mean biomass during the age 

h) The Catch in number during the age 

i) The Catch in weight during the age 

2. Determine: 

a) The cumulative number of survivors during all the exploitable phase 

b) The cumulative biomass of the cohort during all the exploitable life 

c) The total catch, in numbers, from the cohort during all the exploitable life; 

d) The total catch, in weight, from the cohort during all the exploitable life; 

e) The mean weight of the individuals caught, during all the exploitable life; 

1) The mean weight of the individuals of the cohort, during all the exploitable life. 

GROUP 11 

Present the histograms of : 

1 . Mean numbers of survivors of the cohort in each age, during all the exploitable life. 

2. Mean biomasses of the cohort at each age, during all the exploitable life. 

3. Catches in number by age of the cohort, during all the exploitable life. 

4. Catches in weight by age of the cohort, during all the exploitable life. 

GROUP III 

Suppose now, that one intends to analyse the case in which all the fishing mortality 
coefficients were 30 percent bigger than those indicated in the table of Group I. 

1. Calculate the characteristics of the Group 1-2 that could be obtained from the same 
cohort during all its life, in this alternative situation. 

2. Compare the values of the characteristics obtained under these conditions with those 
obtained in Group 1-2, calculating the percentage of variation of those characteristics in 
relation to the corresponding values of the previous situation. 
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8.9 COHORT DURING ITS LIFE - SIMPLIFICATION OF (3.7) 

BEVERTON AND HOLT MODEL 

GROUP I 

The recruitment to the exploitable phase of horse-mackerel, Trachurus trachurus , in Ibero- 
Atlantic waters (Div. VIIIc and IXa) occurs at age 1 . 

The recruitment at the exploitable phase was simplified by adopting the age tc=2 year. 

The parameters of the von Bertalanffy equation estimated for this stock are the following : 

L» = 34.46 cm 
K - 0.225 year 1 
to = -1.66 year 

The weight-length relation: W(g) = 0.01 1 L(cm) 2 90 
The mortality of this stock is characterized as : 

• A constant natural mortality coefficient during all the exploitable phase 
M = 0.15 year' 1 ; 

• Fishing mortality coefficient, F = 0. 14 year 1 , constant during all the exploited phase. 

1 . Calculate, using the simplification of Beverton and Holt : 

a) The recruitment R c to the exploited phase. 

b) The number of deaths during all the exploited life. 

c) The cumulative number of survivors during all the exploited life. 

d) The cumulative biomass during all the exploited life. 

e) The catch in number during all the exploited life. 

f) The catch in weight during all the exploited life. 

g) The mean weight of the individuals of the cohort during all the exploited life. 

h) The mean weight of the individuals caught during all the exploited life. 

GROUP II 

The data presented in Section 8.8 - Group 1 show a great variety in the values of F. However, 
in Section 8.9, the simplification of an F constant was adopted, F = 0.14 year' 1 . The purpose, 
now, is to compare the results from Section 8.8 with those obtained in this exercise, using the 
simplification of Beverton and Holt. 
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1 . 


a) Determinate the cumulative number, the cumulative biomass, the total catch in 
number, the total catch in weight and the mean weight in the catch. 

b) Compare the results with those obtained in Section 8.8 (Group 1-2). 


a) Using a 30 percent bigger value of F calculate the cumulative number, the 
cumulative biomass, the total catch in number, the total catch in weight and the mean 
weight in the catch and compare with the corresponding values obtained in the 
Group I-l.a). 

b) Compare the percentages of variation of the characteristics obtained in 2.a) with 
those obtained in Section 8.8 (Group III). 


8.10 STOCK - SHORT-TERM PROJECTION 


(4.3) 


The mortality parameters and the exploitation pattern of the Iberic Stock of hake, Merluccius 
merluccius, from Divisions VIIIc and IXa of ICES were estimated by the Assessment 
Working Group of ICES (ICES, 1998b) as : 

• Natural Mortality Coefficient = 0.20 year 1 

• Fishing Mortality Level in 1996 = 0.24 year' 1 


Exploitation pattern in 1996 


age 

0 

1 

2 

3 

4 

5 

6 

7 

8 

S| 

0.00 

0.09 

0.29 

1.31 

1.25 

1.12 

1.32 

1.55 

1.55 


The growth parameters were estimated as being : 


Growth parameters of von Bertalanffy 
(CE, 1994): 

Weight length relation, 

W,(g)= a.L(cm) b (Cardador, 1988): 

= 100 cm 

a = 0.004 

K = 0.08 year' 1 

b = 3.2 

/ o =-1.4 year 



At the beginning of 1996 the considered stock had the following age structure, with i 
representing the age and N, the number of survivors at the beginning of the age i, expressed in 
millions of individuals : 


age 

0 

1 

2 

3 

4 

5 

6 

7 

8 

*r, 

83 

27 

41 

30 

22 

11 

6 

3 

5 


127 


Copyrighted material 




























GROUP I 

1. Estimate the individual weights at the beginning of each age, the total number of 
individuals and the total biomass of the stock at the beginning of the year. 

2. Estimate the individual mean weight, the mean number of survivors, the mean biomass, 
the total catch in number, the total catch in weight and the mean weight in the catch and in 
the stock, during 1996. 

3. Supposing there was a recruitment in 1997 equal to 100 million individuals, calculate the 
following, for the beginning of 1997: 

a) The age structure of the stock, in number 

b) The age structure of the stock, in biomass 

c) The total number of individuals of the stock 

d) The total biomass of the stock 

(notice that the stock at the beginning of 1997 is equivalent to the stock at the end of 1996, 
except for the 1997 recruitment). 


GROUP II 

The scientists drew attention to the fact that the fishing level in 1996 was very high and 
should, therefore, be reduced. They suggested a reduction of about 40 percent so that one 
could get adequate catches in weight and biomasses in the future. 

As an alternative to improve the exploitation of this stock, they also suggested increasing the 
mesh size of the fishing nets. 

1 . Fishing managers asked the scientists to evaluate the mean biomass of the stock, the catch 
in number, the catch in weight and the mean weight in the catch during 1997 in case of : 

a) Maintaining the 1996 fishing pattern (status quo situation) 

b) Reducing the 1996 fishing level by 40 percent. 

Making the necessary calculations in order to estimate a) and b), present the results and 
comment on the changes in the catch in weight and mean biomass, relative to 1996. 


128 


Copyrighted material 



2. Fishing managers also asked the scientists to evaluate the catch in weight and the mean 
biomass that would result from maintaining the fishing level of 1996 in 1997, but using a 
new mesh size of the fishing net. To solve this question, the scientists proposed the 
following new exploitation pattern: 


age 

0 

1 

2 

3 

4 

5 

6 

7 

8 

Si 

0.00 

0.00 

0.03 

1.50 

1.65 

1.75 

1.80 

1.80 

1.80 


Make the calculations, present the results and make your comments on the changes in the 
catch in weight and the mean biomass relative to 1996. 

3. Evaluate the effects on the catch in weight and on the mean biomass, resulting from the 
simultaneous adoption (in 1997) of the reduction by 40 percent of the fishing level in 
1996 and the introduction of the new exploitation pattern. Present the results and make 
your comments . 


8.1 1 STOCK - LONG-TERM PROJECTION 


(4.4) 


The mortality parameters of the Iberic Stock (Div. Vlllc e IXa) of Sardine, Sardina 
pilchardus, were estimated by the Assessment Working Group of ICES (ICES, 1997b) as : 

• Natural Mortality Coefficient = 0.33 year ' 1 

• Fishing Mortality Level in 1996 = 0.56 year 


Assume that the exploitable phase occurs from the beginning of age zero until the end of age 
six. 


Relative pattern of exploitation, si, during 1996 


s„ 

Sl 

S2 

S3 

s< 

*5 

S6 

0.21 

0.41 

0.79 

1.18 

1.34 

1.43 

1.68 
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The parameters of the individual growth and the weight length relation of this stock were 
estimated (Pestana, 1989) as : 


Von BertalanfTy 
growth parameters 


L™ = 22.3 cm 
K = 0.40 year' 1 

to = - 1 .6 year 


Weight length relation 

W = a.L b 

with L in cm and W in g 

a = 0.0044 
b = 3.185 


GROUP I 

1. Calculate the evolution in number and the biomass of a cohort during its life at the 
beginning of each age, supposing that growth, natural and fishing mortality parameters are 
the given values. 

2. The recruitment at age 0 in 1996 was estimated as being 4300 million individuals. 
Calculate the cumulative number, the cumulative biomass, the catch in number and in 
weight, during all the life of the cohort. 


GROUP II 

It was estimated that, at the beginning of 1996, the considered stock had the following age 
structure, in number, representing i the age and Nj the number of survivors at the beginning of 
age i, expressed in millions of individuals : 


1 


1 

2 

3 

4 

5 

6 

N, 

4300 

279 

591 

233 

561 

384 

180 


To make a long-term projection of the stock, all the mortality and growth parameters will be 
considered stable during the following years. Consider, also, that the recruitment in those 
future years will be equal to the recruitment of 1996. 

1 . Based on these assumptions, project the numbers of survivors at the beginning of each 
year and age, until 2006. 

2. Compare the structures of the stock in the years 2003 and 2006. 

3. Compare the evolution of the 1996 cohort with the structure of the stock at the beginning 
of 2003. 
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8.12 STOCK - RECRUITMENT RELATION 


(4.5) 


During the meeting of the assessment Working Group (ICES, 1998b) of the Iberic stock 
(Div. VIIIc and IXa) of hake, Merluccius meriuccius, the following stock parameters were 
estimated for the period 1982-1996: 


Year 

N (age 0) 
(million) 

Spawning Biomass 
(thousand tons) 

1982 

125 

59.8 

1983 

107 

61.4 

1984 

136 

58.8 

1985 

97 

44.1 

1986 

104 

26.4 

1987 

97 

24.2 

1988 

84 

22.8 

1989 

56 

18.9 

1990 

59 

19.4 

1991 

69 

20.5 

1992 

86 

21.5 

1993 

70 

21.0 

1994 

63 

16.5 

1995 

32 

15.2 

1996 

83 

18.0 


The annual recruitment at the exploitable phase is considered as being the number of 
individuals with age 0. 

1. Draw the dispersion graph of the resulting recruitments, against the parental spawning 
biomass. 

2. The parameters of Shepherd, Ricker, Beverton & Holt and Deriso S-R relations were 
estimated and are shown in the following table, as well as the respective determination 
coefficients, r : 


Parameters 

Shepherd 

Ricker 

Beverton & Holt 

Deriso 

cc(R/g) 

3.50 

4.43 

4.91 

4.40 

k ( 1 000 tons) 

64.94 

78.13 

45.39 

106.27 

C 

3.52 



0.896 

r 

0.71 

0.68 

0.66 

0.75 
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a) Calculate the expected recruitments in each year of the period 1982-1996 using the 
four S-R models. 

b) The determination coefficient, r, can be used as an indicator of a good or bad 
adjustment of the model to the observed data, depending on the number of 
observations (r 2 can be interpreted as a percentage of variation of the observed points 
that is explained by the model. Values close to 1 indicate a good adjustment and 
values close to zero indicate a bad adjustment). Using this indicator, make your 
comments about the adjustments of each model. 


8.13 F m „ (5.2.1) 

Consider the stock of cod fish, Cadus morhua, in the Irish sea (Div. Vila). The following 
mortality and biological parameters were estimated by the Working Group of ICES (ICES, 
1998c): 

Natural Mortality Coefficient: M = 0.20 year' 1 

Fishing Level in 1996: F% = 0.58 year' 1 

Mean weight (kg) in the catch and in the stock : 


Age 

0 

1 

2 

3 

4 

5 

6 

7 

Wi 

0.001 

0.883 

1.778 

3.597 

5.695 

7.904 

8.502 

9.200 


GROUP I 


Assume that the stock was constituted by the age groups 0 to 7 years and that it was exploited 
with the following exploitation pattern : 



0 

i 

2 

3 

4 

5 

6 

7 

s i 



0.96 



1.00 

1.00 

1.00 


1. Calculate, for 1000 recruits, the long-term annual catch in weight and the annual mean 
biomass, corresponding to the fishing level of 1996. 

2. Adopting the factor Ff ac ior between 0 and 2.5 year' 1 with intervals of 0.1 year' 1 draw the 
curve of the annual catch in weight against Ff aclor . Represent, in the same graph, the curve 
of the mean biomass against Ff ac , or . 

3. Calculate the biological reference point Fma*. 
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GROUP II 


1 . Answer the previous questions, but considering the following exploitation pattern. 


Relative pattern of exploitation adopted by the Working Group: 


Age 

0 

1 

2 

3 

4 

5 

6 

7 


0.00 

0.20 

0.96 

1.30 

1.12 

0.67 

0.58 

0.58 


a) Compare the value of F max of this exercise with the one obtained in the exercise in 
Group I. 

b) Calculate the biological reference point Fm,* knowing that the Working Group considers 
the group of age 7 as a group of cumulative ages (7+). The mean weights in the catch and 
in the stock presented at the beginning of the text are maintained, except for the last 
group of ages which will now be the age group 7+ with the mean weight equal to 10.873 
kg. The exploitation pattern is also maintained with the value 0.58 for the group 7+. 

c) Compare the two calculated values of F m „, considering the last age group as 7 or as 7+. 


8.14 F 0 ., 


(5.2.2) 


Consider the Iberic stock (Div. VlIIc e IXa) of four-spot megrim, Lepidorhombus boscii. 
The following parameters were estimated by the Working Group of ICES (ICES, 1998b): 
Natural Mortality Coefficient: M = 0.20 year' 1 

Fishing mortality in 1996: F 96 = 0.36 year' 1 

Exploitation pattern: 


Age 

1 

2 

3 

4 

5 

6 

7+ 

s i 

0.06 

0.43 

0.89 

1.65 

1.66 

1.22 

1.22 


Mean weight (kg) in the catch and in the stock: 


Age 

1 

2 

3 

4 

5 

6 

7+ 

Wj 

0.037 

0.067 

0.086 

0.109 

0.144 

0.188 

0.244 
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1. Adopting the value 1000 for the recruitment to the fishing area, calculate the long-term 
annual catch in weight, the annual mean biomass, and the annual mean weight in the 
catch, corresponding to the fishing mortality level of 1996. 

2. Adopting the factor Ff, clof between 0 and 2.5 year' 1 with intervals of 0.1 year 1 draw' the 
curve of the annual catch in weight against Ff uctor . Represent, in the same graph, the curve 
of the mean biomass against Fr actot . 

3. Calculate the biological reference point Fo i 

4. Calculate the biological reference point F IMX 

5. In the graph of the curves of the annual catch in weight and of the mean biomass against 
F factor* mark the Biological reference points Fo.i and F max which were previously 
calculated. Make your comments. 

6. Calculate the long-term mean biomass, the catch in weight and the mean weight in the 
catch for Fo.i. Compare those characteristics with the values obtained in question 1 and 
make your comments. 


8.15 F raed and F msv (5.2.3) and (5.2.4) 

During the meeting of the Assessment Working Group (ICES, 1998b) on the Iberic stock 
(Div. Vlllc and IXa) of hake, Merluccius merluccius, the following population parameters 
were estimated : 

Natural mortality coefficient: M = 0.2 year. 

Fishing mortality in 1996: F 96 = 0.24 year’ 1 


Mean weight in the catch (g): 



For the long-term projections, the Working Group adopted the mean exploitation pattern for 
the period 1994-1996, as shown in the next table : 


Age (year) 

0 

1 

2 

3 

4 

5 

6 

7 

8+ 

s i 

0.001 

0.11 

0.398 

1.3 

1.261 

1.019 

1.473 

1.874 

1.874 
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1. Calculate the long-term annual mean biomass, the spawning biomass, the annual catch in 
weight and the mean weight in the catch. 

2. Draw the curve of the annual catch in weight and the mean biomass against F, for values 
of Ffacior between 0 and 2.5 year' 1 , with intervals of 0.1 year' 1 . Calculate the TRPs F rrux 
and Fo i- 

3. Calculate the biological reference point F roc d , knowing the recruitments (million 
individuals at age 0) and the spawning biomasses (thousand tons) between 1982 and 1996, 
estimated by the Working Group, presented in the following table : 


Year 

N (age 0) 
(million) 

Spawning biomass 
(1000 tons) 

1982 

125 

59.8 

1983 

107 

61.4 

1984 

136 

58.8 

1985 

97 

44.1 

1986 

104 

26.4 

1987 

97 

24.2 

1988 

84 

22.8 

1989 

56 

18.9 

1990 

59 

19.4 

1991 

69 

20.5 

1992 

86 

21.5 

1993 

70 

21.0 

1994 

63 

16.5 

1995 

32 

15.2 

1996 

83 

18.0 


4. Calculate the long-term mean biomass, the spawning biomass, the catch in weight and 
the mean weight in the catch for F m «i. Compare with the values obtained in question 1 
and make your comments. 

5. Adopting the Ricker stock-recruitment relation, estimated by the Working Group (a = 
4.43 R/ kg and 1C = 78.13 thousand tons) , calculate F M sy, B M sy and Y M sy- and 
compare the different F-target estimated. 
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8.16 MBALANDB, 


(5.3.4 & 5.3.5) 


loss 


During the meeting of the Assessment Working Group (ICES, 1998b) on the Iberic stock 
(Div. VIIIc and IXa) of hake, Mertuccius merluccius, the recruitment (million individuals at 
age 0) and the spawning biomass (thousand tons) was estimated for the period 1982-1996. 
The values obtained are presented in the following table : 


Year 

N (age 0) 
(million) 

Spawning biomass 
(1000 tons) 

1982 

125 

59.8 

1983 

107 

61.4 

1984 

136 

58.8 

1985 

97 

44.1 

1986 

104 

26.4 

1987 

97 

24.2 

1988 

84 

22.8 

1989 

56 

18.9 

1990 

59 

19.4 

1991 

69 

20.5 

1992 

86 

21.5 

1993 

70 

21.0 

1994 

63 

16.5 

1995 

32 

15.2 

1996 

83 

18.0 


1. Using the spawning biomasses and the resulting recruitments, calculate the Biological 
Reference Limit Points, MBAL and B|„ SS . 
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8.17 F| OSS and F cr , sh 


(5.3.5 & 5.3.6) 


GROUP I 

During the meeting of the Assessment Working Group (ICES, 1998b) on the Iberic stock 
(Div. VIlIc and IXa) of hake, Merluccius merluccius, the recruitment (million individuals at 
age 0) and the spawning biomass (thousand tons) was estimated for the period 1982-1996. 
The values obtained are shown in the following table : 


Year 

N (age 0) 
(million) 

Spawning biomass 
(thousand tons) 

1982 

125 

59.8 

1983 

107 

61.4 

1984 

136 

58.8 

1985 

97 

44.1 

1986 

104 

26.4 

1987 

97 

24.2 

1988 

84 

22.8 

1989 

56 

18.9 

1990 

59 

19.4 

1991 

69 

20.5 

1992 

86 

21.5 

1993 

70 

21.0 

1994 

63 

16.5 

1995 

32 

15.2 

1996 

83 

18.0 


The Shepherd S-R relation was adjusted to the pairs of values in the table (r=0.71), and the 
relation parameters are the following : 

a =3.5 Kg' 
k = 64.94 thousand tons 
c = 3.52 

1. Draw the dispersion graph of the resulting recruitments, against the parental spawning 
biomasses. 

2. Calculate the expected recruitments in each year of the period 1982-1996 according to the 
Shepherd S-R model and in the previous dispersion graph, draw the respective curve. 
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GROUP II 


1 . Calculate the annual catch in weight and the spawning biomass per recruit for the stock of 
hake, using the mortality and biological parameters estimated by the Working Group for 
the long-term projections (given in Section 8.15), namely : 

Natural mortality coefficient : M = 0.2 year' 1 

Fishing mortality in 1996: F96 = 0.24 year' 1 

Mean weight in the catch (g): 


Age (year) 

0 

1 

2 

3 

4 

5 

6 

. 

8+ 

Wj 

4 

37 





706 


1508 


Maturity ogive ( percent): 


Age (year) 

D 

1 

2 

3 

4 

5 

6 

7 

8+ 

% maturei 

1 0 

0 

D 

6 

20 

49 

76 

91 

100 


Mean relative pattern of exploitation of the period 1994-1996: 


Age (year) 

D 

1 

2 

3 

4 

5 

D 

D 

8+ 

S j 

0.001 

0.11 

0.398 

1.3 

1.261 

1.019 

1.473 

1.874 

1.874 


GROUP III 

1. Using the results of Groups 1 and II, calculate biological reference Limit-Points, F] oss and 
Fcrash- 
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8.18 PRODUCTION MODELS (EQUILIBRIUM) -SCHAEFER (6.7.1) 


The following table presents the annual total catch Y, (t) and mean biomass B, (t) for the 
fishery of the Iberic stock (Div. VlIIc and IXa) of sardine, Sardina pilchardus, between 1977 
andl996, used by the Working Group of ICES (ICES, 1998a). 


Year 

Y (t) 

B(t) 

1977 

125750 

750289 

1978 

139990 

759192 

1979 

153441 

763313 

1980 

191682 

804765 

1981 

214133 

842091 

1982 

204504 

802573 

1983 

181139 

713376 

1984 

202686 

794856 

1985 

204107 

810539 

1986 

180606 

679808 

1987 

168825 

547179 

1988 

158540 

481295 

1989 

137126 

431719 

1990 

139157 

368099 

1991 

127756 

316365 

1992 

126054 

453161 

1993 

138795 

539096 

1994 

132800 

416842 

1995 

121384 

368158 

1996 

111431 

246037 


GROUP I 

1 . Calculate F, corresponding to each year i. 

2. Calculate the biomasses, B,. at the beginning of each year. (Use the procedure proposed 
by Schaefer, that is, the biomass at the beginning of a year is approximated by the 
arithmetic mean of the mean biomasses of the previous and the following year). 

3. Calculate the equilibrium catches. Ye, which would correspond to the observed values 
ofF. 
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4. Calculate the equilibrium mean biomasses, Be, which would correspond to the observed 
values of F. 

5. Draw the graph of Be against F;. 

GROUP II 

The Schaefer model was adjusted and the following parameters were estimated : 

k= 1562851 t 
r = 0.426 year' 1 

1. Calculate the equilibrium biomasses and the equilibrium catches corresponding to the 
fishing levels observed in each year, using the Schaefer model. 

2. Draw, in the graph of Group I point 5, the equilibrium biomasses calculated with the 
Schaefer model. 

3. Calculate Fmsv , Bmsy and Ymsy 

4. Calculate Fo i , Bo i and Y 0 .|. 


8.19 PRODUCTION MODELS (EQUILIBRIUM) (6.7.1 & 6.7.2) 


The following table presents the annual total catches and the corresponding total fishing 
efforts of a shrimp stock in the Arabian Sea during the period 1969 to 1978 (Sparre and 
Venema, 1992). 


Year 

Catch, Y 
(0 

Total effort 
(1000 days) 

1969 

546.7 

1.224 

1970 

812.4 

2.202 

1971 

2493.3 

6.684 

1972 

4358.6 

12.418 

1973 

6891.5 

16.019 

1974 

6532.0 

21.552 

1975 

4737.1 

24.570 

1976 

5567.4 

29.441 

1977 

5687.7 

28.575 

1978 

5984.0 

30.172 


GROUP I 

1 . Draw a graph of the annual abundance index against the corresponding fishing mortality 
index. 

2. The Fox model was adjusted to the data. The following parameters were obtained : 

a = 6.150 b = - 0.028 with a determination coefficient, r 2 = 0.78. 
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a) Calculate and draw the curves of the equilibrium conditions of the abundance index 
and of the total catch against fishing effort. 

b) Determine the target-points MSY and 0. 1 . 

c) Determine the parameters (kq), (r/q) and (kr). 

GROUP 11 

Knowing that the adjustment of the Schaefer model to the same set of data, resulted in the 
following values of the parameters : 

a = 444.454 
b = - 8.055 
(r 2 = 0.77) 

1 . Repeat the calculations of the previous Group. 

8.20 PRODUCTION MODELS - SHORT-TERM PROJECTION (6.8) 

The following table presents the annual total catches (in tons) and the corresponding catches 
by fishing unit effort (kg/ fishing day of the fleet PESCRUL) of the stock of Deepwater rose 
shrimp, Parapenaeus longirostris, of the Algarve during the period 1983 to 1994 (Mattos 
Silva, 1995). 


Year 

Y(t) 

cpue 

(kg/day) 

1983 

538 

235 

1984 

638 

131 

1985 

431 

63 

1986 

99 

22 

1987 

37 

8 

1988 

62 

21 

1989 

437 

77 

1990 

146 

28 

1991 

126 

26 

1992 

53 

25 

1993 

91 

41 

1994 

232 

66 


Using the following parameters of the integrated Fox model and Yoshimoto & Clarke (1993): 
k= 1580 t, q=0.39 thousand days' 1 and r=0.55 year' 1 : 
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1. Project the cpue and the catch in weight for the year 1995, supposing that the fishing 
effort will be maintained equal to that of 1994 (situation status quo). 

2. Determine the target reference points, Ymsy. Bmsy and Fmsy and the indices Umsy an d 

fMSY- 

3. Determine the target reference points, Yo.i, Bo i and Fn i and the indices Uo i and fo |. 

4. Determine the percentages of the carrying capacity corresponding to the target points 
Fmsy and Fo.i 

5. Based on the results obtained in the previous questions, make your comments on the state 
of the stock and of its exploitation. 

6. Suppose that one intends to reduce the fishing effort in 1995 by about 20 percent relative 
to the effort of 1994. Project the catch in weight for 1995 and present the variations 
resulting from that effort reduction on the catch in weight and on the biomass. 


8.21 SIMPLE LINEAR REGRESSION - ESTIMATION OF THE (7.2) 
PARAMETERS OF THE W-L RELATION AND GROWTH 
PARAMETERS (FORD-WALFORD, GULLAND AND HOLT 
AND STAMATOPOULOS AND CADDY) 


GROUP 1 

Consider the following 10 pairs of values of x and y: 


*i 

2 

6 

7 

8 

11 

15 

16 

18 

19 

21 

yt 

13 

40 

52 

56 

78 

105 

111 

130 

132 

149 


1 . Estimate the constants A and B of the straight line. 

2. Estimate the values of Y corresponding to the given values of x. 

3. Calculate the determination coefficient r z . 

4. Draw a graph with the observed values and with the estimated line. Observe the 
adjustment and say if you consider the linear model adequate. 
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GROUP II 


The data presented on the following table represent the individual weights by length class of 
European anglerfish, Lophius budegassa samples from the Iberic coast in 1994. 


Li (cm) 

Wt(g) 

Li (cm) 

W«(g) 

20- 

129 

50- 

1685 

22- 

163 

52- 

1896 

24- 

219 

54- 

2107 

26- 

265 

56- 

2345 

28- 

320 

58- 

2569 

30- 

397 

60- 

2848 

32- 

486 

62- 

3126 

34- 

545 

64- 

3407 

36- 

664 

66- 

3700 

38- 

773 

68- 

4056 

40- 

890 

70- 

4411 

42- 

1027 

72- 

4764 

44- 

1122 

74- 

5203 

46- 

1334 

76- 

5587 

48- 

1503 

78- 

5982 


I. Using the simple linear regression model, estimate the parameters of the weight-length 
relation for this stock, considering that W CM n, a i=W mta „. (Notice that the In W ccm „! is linear 
to the In L ccn aai). 


GROUP III 


The data presented in the following table represent the mean length (cm) of the individuals at 
the beginning of the age (years), obtained from direct age reading of the individuals of the 
stock of European anglerfish, Lophius budegassa, (ICES Div. VlIIc and IXa), Section 8.7 


T 

L, 

t 

L, 

1 

9.2 

7 

44.4 

2 

16.5 

8 

49.0 

3 

22.9 

9 

52.3 

4 

28.8 

10 

55.0 

5 

34.7 

11 

60.8 

6 

38.6 

12 

63.4 
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From this data and using the simple linear regression model, estimate the growth parameters 
K and L» with the expressions : 

1 . Ford-Walford (1933-1 946) 

2. Gulland and Holt (1959) 

3. Stamatopoulos and Caddy (1989). (The expressions for 1), 2) and 3), were studied in the 
Chapter concerning Individual Growth). 

4. Comment on the results obtained in the previous questions with the values of the 
parameters given in Section 8.7. 


8.22 MULTIPLE LINEAR MODEL - REVISION OF MATRICES - (7.3) 
ESTIMATION OF THE PARAMETERS OF FOX 
INTEGRATED MODEL (IFOX) 


Revision of matrices 


group i 


Consider the matrices A and B : 



A = 

230 1 

B = 

1103 



114 1 


1325 



0422 


2 160 



1033 


22 10 

1. 

Using a spreadsheet, calculate: A + B, A * B, Det(A), Det(B), A' 1 2 3 e B 

2. 

Show that (A.B) 1 

= B’'.A 1 



3. 

Show that (A.B) t 

II 

to 

H 

GROUP II 


Let the Matrices: 




M(44)— 

(1/4) 1 1 1 1 

0(4,4) _ 

0000 

1(4,4) ~ 10 0 


1 1 1 1 


0000 

0 1 0 


1111 


0000 

00 1 


till 


0000 

000 


1 . Verify that the matrix null 0 is idempotent. 

2. Verify that the matrix identity f is idempotent. 

3. Verify that the matrix M is idempotent. 
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4. What are the traces of M and I? 

5. Calculate the ranks, r , of M and of I. 

6. What is the value of the determinant of M and I? 

GROUP III 

1. Verify that the product Mx, where x is the vector given by x T =(3 4 8 1), is a vector 
with all elements equal to the arithmetic mean, x, of the 4 elements of vector x. 

2. Verify that (I-M)x is the vector of the deviation. 

3. Verify that the sum of the squares of x s , £(x 2 ) can be written as : x T . x 

4. Verify that the sum of the squares of the deviations, £(x, - x) 2 , can also be written in 

a matricial form, as : x r (I-M)x 

GROUP IV 

1 . Consider the vector x = 2 + 0 where 0 is an unknown parameter. 

30 

5-0 

a) Write the derivative — of the vector x 

dO 

b) Calculate x T x 

d t 

c) Calculate — (x x) 

dO 

d) Show that — (x T x) =2(— ) r x 

dO dO 

2. Consider the vector x = 2 + 40, - 502 where 0i and 02 are two unknown constants. 

1 + 0 | + 02 

0i : +402 

a) Write the derivative matrix (take 0i and 02 as variables) 

aO 

b) Calculate x r x 

c) Transpose (x T x) 
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d dx 

d) Show that the transposed matrix — (x T x) = 2 ( — ) r x 

dO dO 


GROUP V 

Consider the following system of 2 equations with 2 unknowns; 

5 = 2 A + 3 B 
4= A- 2 B 

1 . Show that the equation system can be written in matrix form as, 

Y(2,l) = X(2,2) 6(2,1) 

where Y is the vector of the independent terms (5 e 4) of the system, 

0 is the vector of the unknowns A and B 

and X is the matrix of the coefficients of the unknowns 

2. Verify that the solution of the system can be given as 0 = (X r X)'X T Y 

3. Show that X is a square, non singular matrix, and then that the solution of the system can 
be 0=X'Y 

Estimation of the parameters of the Yoshimoto and Clarke model 

( 1993 ) 

4. Estimate the parameters k, q and r, of the Fox integrated model (IFOX) and of 
Yoshimoto & Clarke (1993) using the following data : 


Year 

Y 

(t) 

CPUE 

(kg/day) 

1983 

538 

235 

1984 

638 

131 

1985 

431 

63 

1986 

99 

22 

1987 

37 

8 

1988 

62 

21 

1989 

437 

77 

1990 

146 

28 

1991 

126 

26 

1992 

53 

25 

1993 

91 

41 

1994 

232 

66 
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which represent the total annual catches (in tons) and the respective catches by fishing effort 
unit (kg/fishing day of the fleet PESCRUL) of the stock of Deepwater rose shrimp, 
Parapenaeus longirostris of the Algarve during the period 1983 to 1994 (Mattos Silva, 1995). 

Comment on the obtained results comparing them with those presented in Section 8.20. 


8.23 NON LINEAR REGRESSION - ESTIMATION OF THE (7.4) 

GROWTH PARAMETERS AND OF THE S- R RELATION 
(GAUSS-NEWTON METHOD) 


The data in the following table represent the mean length (cm) of the individuals at the 
beginning of the age (years), obtained from direct age reading of individuals of the stock of 
anglerfish, Lophius budegassa, (Div. VlIIc andlXa), Section 8.7. 


t 

L, 

t 

L, 

1 

9.2 

7 

44.4 

2 

16.5 

8 


3 

22.9 

9 

52.3 

4 

28.8 



5 

34.7 

11 


6 

38.6 

12 

63.4 


GROUP 1 


1 . Represent graphically the values of L, against t. 

2. Estimate the growth parameters K, L„ and to for the Bertalanffy growth model using a 
non-linear regression model. 

3. Estimate the values of L t corresponding to given values of t and mark on a graph the 
observed values and the estimate curve. Comment on the adjustment. 

4. Compare the results obtained in the previous questions with the values of the parameters 
given in Section 8.7. 


GROUP 11 

Using the data on spawning biomass and on recruitments of hake, presented in Section 8.12, 
estimate the parameters of the Beverton and Holt, Ricker, Deriso and Shepherd S-R models. 

1 . Compare the obtained values with those presented in Section 8.12. 
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8.24 ESTIMATION OF M 


(7.6) 


GROUP I 

1 . Estimate the 5 percent Tanaka survival curve for natural mortality coefficients between 
0.0 and 5.0 year 1 and longevities between 0 and 30 years. 

2. Calculate the values of M corresponding to the longevity 1, 2, 3, 10, 15 and 30 years. 

GROUP II 

The parameters of von Bertaianffy equation for the Iberic stock (Div. VIIIc and IXa of ICES) 
of horse-mackerel, Trachurus trachurus, are the following : 

= 34.46 cm, TL (carapace length) 

K. = 0.225 year' 1 

1. Estimate the value of M for this stock, knowing that from 1985-95 in the distribution area 
of this stock in the Iberic Peninsula, the mean temperature of the sea water at the surface, 
was : T = 13°C 


GROUP III 

Consider a certain fishing stock for which the first maturity mean age was estimated as 2.3 
year. 

1 . Obtain an approximate estimate of M for this stock. 

GROUP IV 

The reproductive biology of the stock of Atlantic mackerel, Scomber scombrus was studied 
and it was estimated that the mean gonadsomatic index (gonad weight/total weight) of the 
mature females in the spawning period was 0.13. 

1. Estimate an approximate value for the natural mortality coefficient M for this stock, 
assuming that it is constant for all the ages and years. 


GROUP V 

The scientists responsible for the evaluation of the stock of a certain fishing resource, 
implemented an acoustic cruise every year in January, to estimate the abundance of the stock 
by age classes. Fishery statistics are also needed to estimate the catch by age during the year. 

The following table presents the estimations of the abundance of the stock, by age classes, 
obtained on the cruises in 1993 and 1994, as well as the structure of the catches during 1993. 
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Age 

Number of survivals, N, in million 

Total catch in 1993 


January 1993 

January 1994 

(million) 

2 

243 

353 

11 

3 

99 

189 

15 

4 

86 

67 

20 

5 

37 

52 

9 

6 

13 

22 

3 

7 

6 

8 

1 


Although this information is available, the scientists responsible for this stock have 
difficulties applying the assessment models to the resource because they do not have an 
estimation of M. 

1. So, estimate the natural mortality coefficient for this resource in 1993, and help the 
group of scientists responsible for the resource. 


GROUP VI 

The following table presents the data on fishing effort, in million trawl hours, and the total 
mortality coefficent, Z, for a certain fishery for the period 1987 to 1995. 


Years 

Effort 
(10* hours) 

Z (year ') 

1987 

2.08 

1.97 

1988 

2.80 

2.05 

1989 

3.5 

1.82 

1990 

3.6 

2.32 

1991 

3.8 

2.58 

1992 

- 

- 

1993 

- 

- 

1994 

9.94 

3.74 

1995 

6.06 

3.74 


1. Determine M (natural mortality coefficient), assumed to be constant for the period 1987- 
1995. 

2. Determine the catchability coefficient q. 
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1. Estimate the total mortality coefficient, Z, of the stock in each of the periods. 

2. Consider that the mean fishing effort on the North Sea plaice during the two periods was 5 
million hours of trawl in 1929-1938 and 3.1 million hours of trawl in 1950-1958. Estimate 
for each period: 

a) the natural mortality coefficient, M; 

b) the catchability coefficient, q; 

c) and the fishing mortality coefficient, F. 


GROUP III 

The following table presents the annual composition of the catches by age from 1988 to 1994, 
in millions of individuals, for a certain resource : 
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CATCHES (million individuals) 


Age 

1988 

1989 

1990 

1991 

1992 

1993 

1994 

0 

599 

239 

424 

664 

685 

478 

330 

1 




1004 

418 


288 

2 




532 

335 


323 

3 



159 

269 

203 



4 


mm 

75 

32 

69 

86 


5 

6 

9 

13 

18 

8 

25 

31 

6 

1 

8 

3 

5 

5 

3 

8 

7 

6 

0 

1 

0 

1 

1 

1 


1 . Calculate the mean annual composition during 1988-1 994, 

2. Estimate Z, based on that mean composition. 

3. Estimate Z, based on the mean age of the mean composition of the catch. 

4. Estimate Z for each year of the given period. 

5. Compare the annual Zs with the values of Z obtained in questions 2 and 3. 

GROUP IV 


The following table shows the length composition, in equilibrium, of a certain resource, with 
L x = 100 cm and K = 0.2 year 1 2 3 4 5 . 


Length class (cm) 

35- 












95- 

Catch (Ci) in million 

7 








33 



17 

8 


1 . Calculate the relative ages corresponding to the lower limit of each length class. 

2. Determine the age interval corresponding to each length class. 

3. From which class can one consider Z constant? 

4. Determine Z using : 

a) The catches in each class. 

b) The cumulative catches. 

c) The mean length in the catch. 

5. Compare the values of Z obtained by the different methods of question 4. 
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GROUP V 


The length compositions of the catches for three different periods of time are known for a 
certain fishing resource. 


Period 

Length 
classes (cm) 

45- 








1 


>95 

1960-69 


256 









H 

12 



268 

226 

180 

141 

105 

76 

50 



6 

1 

1980-89 


212 

161 

116 

79 

52 

31 

17 

8 

3 

1 

0 


Consider the 45 cm length class as the first class completely recruited . 

Adopt K. = 0.3 year ' and Z, M = 100 cm as the von Bertalanffy growth parameters for this 
resource. 

1. Estimate the values of the total mortality coefficient, Z, for each period and comment on 
the results. 


8.26 AGE COHORT ANALYSIS (CA) 

GROUP I 


(7.9.1) 


1 . Consider a stock and an interval of time i, (/,• , r i+) ) . Knowing that for this interval of 
time: 

Mi = 0.4 year' 1 
Tj = 2.3 year 

Ci = 230 million individuals 

a) Adopt the value 0.5 year 1 for the fishing mortality coefficient for the interval and 
calculate the numbers of survivors at the beginning and end of the interval. 

2. Consider the interval of time i, (r, , t, + | ) . Knowing that in this interval of time: 

Mj = 0.6 year 1 
Tj = 0.9 year 

Q = 98 million individuals 
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Calculate the value of the fishing mortality coefficient, Fj , for the interval , taking the number 
of survivors, N, , at the beginning of the interval, i, to be 1 72 million individuals. 

3. Consider the interval of time (t„ t i+i). Knowing that in that interval of time : 

Mj= 0.5 year' 1 
Tj = I year 

Ci = 42 million individuals 

a) Calculate the value of the fishing mortality coefficient for the interval, knowing 
that the number of survivors at the end of the year was N,h= 85 million 
individuals. Calculate the value of F, using the Pope formula. 

GROUP II 

The data in the following table represent the catches in millions, of a cohort of hake, 
Merluccius merluccius, in the Iberic Peninsula waters. 


Age (years) 

0 

1 

2 

3 

4 

5 

6 

7 

8 

Ci (million) 

712 

3941 

8191 

10311 

5515 

4149 

3081 

1185 

549 


Adopt a value of 0.2 year' 1 for the natural mortality coefficient, constant for all the ages. 

1. Suppose that the value of the fishing mortality coefficient at the last age (8 years) was 
1.0 year Calculate, by an iterative method and by the Pope method, for each age of the 
cohort : 

a) The value of the fishing mortality coefficient. 

b) The number of survivors at the beginning of the age. 

c) Compare the results obtained by the two methods. 

d) Represent, on a graph the values of F, estimated against the age, and say what the 
recruitment of this cohort is at the exploited phase. 

GROUP III 

1 . Aiming to analyse the influence of the chosen F Krmma t , repeat the calculations of question 
1 of Group II, using one of the previous methods, with 0.3 and 1.5 year' 1 for the value of 

F terminal • 

a) Draw a graph with the estimated values of F, and N, against the age. 

b) Comment on the differences between the graphs for the different values of F tcml mai. 

2. Aiming to analyse the influence of the choice of M, repeat the calculations of question 1 
of Group II, using one of the previous methods, for values of M of 0. 1 and 0.4 year 1 . 
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a) Represent, on a graph, the estimated values of F, and N; against the age. 

b) Comment on the differences between the graphs for the different values of M. 


CROUP IV 


The annual catches by age class, of a certain resource, for the years of 1985 to 1994, are 
presented in the following table. 



Years 

Age 

(years) 

1985 

1986 

1987 

1988 

1989 

1990 

1991 



1994 

0 

67 

88 





63 

38 

52 

90 

1 

532 

1908 







541 

704 

2 







2826 

1911 

1322 

741 

3 




m 


1972 



1382 

890 

4 



3309 



1170 

883 

2102 

896 

540 

5 

97 

439 

733 

1730 

558 

827 

479 

356 

807 

316 

6 

16 



333 

656 

207 

291 

166 

117 

243 

7 

25 

8 

73 

136 



73 

101 

54 

35 

8 

5 

7 


33 

52 

47 

85 

25 

33 

16 


The modus operandi of the fishing fleet was constant during the period, but the number of 
vessels increased significantly. It is considered that, at present, the resource is intensively 
exploited. 

Besides the information on the fishery, the estimates of the growth parameters of this resource 
and of the natural mortality coefficient are also available: 

L» = 38.5 cm a = 0.02 1 of the relation W(g)-L(cm) 

K. = 0.25 year' 1 b = 2.784 of the relation W(g)-L(cm) 

t„ = -0.51 year M = 0.3 year' 1 

1 . Estimate the fishing mortality coefficient and the number of survivors at the beginning of 
the year for each age class and each year. Use the Pope Cohort Analyses method. 

a) Start by selecting F tcmlma i = 0.5 year' 1 for the last age of every year and for all the ages of 
the last year. 

b) After analysing matrix F obtained in a), select new values for Ftcmmui and repeat the 
application of Pope’s method. 
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2. Besides the information given in the previous question, it is also known that the spawning 
takes place in a restricted period, around the beginning of the year. Research cruises using 
acoustic methods took place during the spawning period, in order to estimate the 
spawning biomass (kg/hour of trawl). The results obtained are shown in the following 
table : 


Years 

1985 

1986 

1987 

1988 

1989 

mill 





Spawning 

biomass 

Index 


1613 










The biological information collected during those cruises was also used to estimate the 
maturity ogive of the stock at the spawning period: 


Age (years) 

D 

1 

2 

» 

4 

5 

6 

7 

8 

% Matures 

0 

1 

20 

50 







a) Calculate the spawning biomass in the spawning period of each year from 1985 to 1994 
using the results of the Cohort Analyses obtained in question 1 .b. 

b) Use the information of the acoustic cruises to tune the Cohort Analyses. 

c) Comment on the tuning results. 


8.27 LENGTH COHORT ANALYSIS (LCA) (7.9.2) 

GROUP I 

The following table presents the annual catch length composition, of a cohort of a resource 
with L» = 1 30 cm and K = 0. 1 year' 1 . 
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Length 

classes(cm) 

Catch, C| 
(million) 

6- 

1823 

12- 

14463 

18- 

25227 

24- 



3889 

36- 

2959 

42- 

1871 

48- 

653 

54- 

322 


228 

66- 

181 

72- 

96 

78- 

16 

X4- 

0 


The natural mortality coefficient was estimated as being M = 0.3 year' 1 . 

1. Using the Pope method, and adopting E=0.5 as being the exploitation rate in the (78-) 
length class of the catch, estimate the number of survivors at the beginning of each length 
class, the fishing mortality coefficient F and the exploitation rate E in each class. 

2. Calculate the mean number of survivors of the cohort. 

GROUP H 

The following tables 1 and 2 present the basic information on a hypothetical stock during the 

years 1985 to 1994. 

1. Apply the slicing technique to the Catch matrix and comment on the validity of applying 
cohort analyses by ages. 

2. Estimate the matrices TF] and [N] by length classes and years. 

3. Calculate the matrix [Fsep] and comment on the hypothesis that the exploitation pattern 
can be considered to be constant during those years. 


Table 1. Growth parameters of the von-Bertalanffy curve, l,„ and K Natural Mortality 
Coefficient, M and constants a and b of the wcight/length relation 


Growth 

Natural Mortality 

Weight/length relation 
w ( = a . (Lj) b 

L» (cm) 

42 

M (year 1 ) 

0.8 

a 

0.0023 

K (year 1 ) 

0.5 



h 

3 
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Table 2. Catch matrix in thousands of individuals, by length classes and years 
in the period 1985-94 


Age 

(sliced) 

Length 

classes 

(cm) 

Years 

1985 

1986 

1987 

1988 

1989 

1990 

1991 

1992 

1993 

1994 

0 


35 

m 


17 


69 

34 

61 

46 

29 


BB 




167 


662 

327 

593 

442 

276 


22- 

805 





1575 

777 



657 


23- 

1500 

1766 

1317 

757 

2108 

2923 

1436 



1220 


24- 

1901 

2222 

1702 

985 

2688 

3678 

1795 

3175 

2485 

1535 


25- 




mm 


3900 




1627 


26- 




am 


3600 


2925 

2482 

1502 

1 

27- 



1228 



2932 

3376 

1695 

1785 

2376 


28- 



1091 

1276 

1250 

2467 

2801 



1999 


29- 




1125 

1053 

2018 

2258 

1071 

1265 

1636 


30- 

1127 

1003 

812 

980 

873 

1619 

1782 

818 

1031 

1312 


31- 


787 

684 



1269 

1372 

607 

823 



32- 






959 

1017 

432 

635 


2 

33- 




389 

834 

511 






34- 




310 

618 

361 




365 


35- 

403 

267 

170 

240 

439 

242 


mm 

110 

245 


36- 

262 

168 


178 

294 

152 

207 

191 

72 

154 

3 

37- 

165 

168 

66 

71 

175 

214 

93 


75 

46 


38- 

86 

84 

40 

45 

96 

107 

44 

55 

39 

23 


Consider L a = 20 cm and t„=0 

(Extracted from : Cadima, E. & Palma, C.,1997. Cohort Analysis from annual length catch 
compositions. Working document presented to the Working Group of the Demersal Stocks 
Assessment of the South Shelf, held in Copenhagen from 1-10 September, 1997.) 
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8.28 EXAMINATION - WRITTEN TEST 


TRAINING COURSE ON FISH STOCK ASSESSMENT 
INSTITUTO DE INVESTIGAQAO DAS PESCAS E DO MAR 
(LISBON, 3 NOV. -10 DEC. 1997) 

QUESTION 1 

Consider a certain stock with the following parameters: 


Natural Mortality Coefficient: M = 0.20 year' 1 

Fishing mortality in 1996: F% = 1.08 year' 1 

Exploitation pattern : 



Age (year) 

1 

2 

3 

4 

5 

6 

7+ 

% matures, 

34 

90 

100 

100 

100 

100 

100 
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The estimated recruitments and the spawning biomasses between 1986 and 1996 are 
presented in the following table : 


Table 1 


Year 

N (age 1) 
(thousands) 

Spawning biomass 
(t) 

1986 

8751 

1957 

1987 

8305 

1591 

1988 

7123 

1956 

1989 

7596 

2073 

1990 

6013 

2287 

1991 

5054 

1506 

1992 

9713 

1400 

1993 

5520 

1275 

1994 

6000 

980 

1995 

7329 

675 

1996 

6840 

917 


1. Calculate the biological reference point Fo i. Indicate the value of the virgin biomass that 
you estimated and calculate the percentageBo i/B^u,.. 

2. Calculate the biological reference point F m , 

3. Calculate the biological reference point F^d 

4. Estimate the parameters of the S-R Ricker model and indicate what the value of the 
spawning biomass/recruit is, corresponding to F clajh . 

5. Comment on the present state of the stock and its exploitation. 
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QUESTION 2 

Consider a stock where the length compositions of the catches during 1986-1995 are known 
(Table 2). 

The following parameters were estimated from this stock : 

Natural mortality coefficient = L2 year 1 
Asymptotic length = 39.8 cm 
Growth coefficient= QJS year ' 


Table 2 - Catch Matrix (thousand individuals) 


Length 
classes (cm) 

Year 

1986 

1987 

1988 

1989 

1990 

1991 

1992 

1993 

1994 

1995 

14- 

m 

2 

8 

2 

2 

18 

13 

12 

15 

8 

15- 

211 

12 

15 

18 

O 

35 

26 

23 

22 

16 

16- 

22 

25 

23 

26 

12 

52 

32 

H 

42 

24 

17- 

42 

41 

12 

43 

31 

84 

63 

56 

62 

32 

18- 

22 

80 

22 

83 

60 

164 

122 

ms 

134 

26 

19- 

224 

194 

175 

202 

146 

398 

297 

263 

325 

185 

20- 

261 

226 

203 

234 

169 

461 

343 

305 

376 

214 

21- 

420 

363 

326 

376 

271 

736 

547 

488 

601 

342 

22- 

335 

525 

603 

506 

511 

365 

983 

571 

657 

741 

23- 

345 

540 

618 

516 

520 

370 

995 

581 

666 

752 

24- 

422 

661 

751 

622 

625 

444 

1189 

699 

799 

902 

25- 

442 

693 

781 

642 

643 

454 

1212 

718 

818 

923 

26- 

415 

650 

726 

592 

590 

415 

1102 

659 

747 

843 

27- 

388 

607 

672 

542 

537 

376 

995 

601 

677 

765 

28- 

360 

564 

617 

493 

486 

339 

890 

543 

609 

688 

29- 

332 

520 

563 

444 

435 

301 

788 

487 

543 

613 

30- 

304 

475 

508 

396 

386 

265 

689 

431 

478 

539 

31- 

275 

430 

453 

348 

337 

230 

594 

377 

414 

468 

32- 

272 

246 

509 

494 

330 

314 

212 

428 

351 

353 

33- 

239 

216 

439 

419 

277 

261 

174 

359 

292 

294 

34- 

206 

186 

369 

345 

226 

211 

132 

293 

235 

237 

35- 

171 

LiS 

300 

273 

176 

162 

106 

228 

181 

182 

36- 

06 

121 

230 

202 

129 

III 

25 

167 

130 

111 

37- 

82 

22 

Ufi 

196 

142 

21 

80 

40 

UL2 

83 

38- 

42 

60 

66 

103 

25 

45 

38 

21 

58 

41 

39- 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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1. Estimate the values of the fishing mortality coefficient in each length class of the cohort 
of 1987 (to simplify, adopt F, er minai = 0.5 year' 1 ). 

2. Estimate the corresponding values of the number of survivors at the beginning of each 
class. 

3. Say what the recruitment of this cohort is. 


QUESTION 3 


The Fox Production model was adjusted to a certain stock and the following parameters were 
obtained : 


K = 300 thousand tons 
r = 0.50 year 1 

1 . Determine the biological reference points Fmsy and Fo |. 

2. Knowing that in recent years the biomass of this stock is about 30 percent of the virgin 
biomass, comment on the present state of the stock and its exploitation, based on the 
adopted model. 
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This manual starts with an introduction to tha mathematical models applied in fish stock 
assessment The basic assumptions about a model and the concepts of variation rates of 
a characteristic in relation to time are presented. The concept of cohort and models for 
the individual growth of the cohort are developed. In the chapter concerning the study 
of the stock, the fishing pattern and its components are defined, the most used models 
for the stock-recruitment relation are presented, as well as the short- and long-term 
projections of a stock. With regard to fisheries resources management the discussion 
is focused on the biological reference points. Finally, the general methods of estimating 
parameters with special rele va nce to the cohort analysis by age and length are descr i bed. 
The exercises from the last course held in the Institute de Investigate das Pescas 
a do Mar, Portugal are pre s ented by the author and the scientist Manuals Azevedo. 
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